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Abstract 

Light Front Field Theory Calculation of Deuteron Properties 

by Jason Randolph Cooke 

Chair of Supervisory Committee: 

Professor Gerald A. Miller 
Department of Physics 


Light front dynamics is a promising approach for solving bound state problems in 
nuclear physics. It is also ideal for calculating the deuteron form factors at high 
momentum transfers. However, in light-front dynamics rotational invariance is not 
manifest, which results in a splitting in the binding energies of states with different 
magnetic quantum numbers and a breaking of the angular condition for the matrix 
elements of the deuteron current operator. The objective of this work is to inves¬ 
tigate the symptoms of broken rotational invariance for deuterons calculated using 
one-meson-exchange (OME) and two-meson-exchange (TME) potentials for various 
models. 

We first consider the Wick-Cutkosky model. The binding energies of states with 
different m values are split when calculated with the OME potential, but this splitting 
is reduced for all states when the TME potential is included. In addition, we find 
that appropriate OME+TME potentials give almost identical results to the ladder 
and crossed ladder Bethe-Salpeter equation. 

Next, we derive light-front nucleon-nucleon OME and TME potentials from an 
effective nuclear Lagrangian. We consider first the potentials generated by the ex- 




change of pions only. There is a large splitting in the binding energies between the 
rn — 0 and m — 1 denteron states when the one-pion-exchange (OPE) potential 
is used. Including the chiral two-pion-exchange (TPE) potential in the calculation 
reduces this splitting. We then proceed to use the potentials generated by all the 
mesons, and find that the deuteron mass splittings are smaller for both the OME and 
OME+TME potentials than in the pion-only model. 

The deuteron wave functions are used to calculate the electromagnetic and axial 
current matrix elements and form factors. The matrix elements of the electromag¬ 
netic current operator have better transformation properties under rotations when 
we use the OME potential instead of the OME+TME potential to calculate the wave 
functions. The axial current matrix elements have about the same transformation 
properties regardless of whether the OME or OME+TME potential is used. Finally, 
at momentum transfers greater than about 2 GeV 2 , the breaking of rotational invari¬ 
ance causes less uncertainty in the A and B form factors than do the uncertainties in 
the nucleon form factors. 
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Chapter 1 

INTRODUCTION 

“If you always do what interests you, at least one person is pleased.” 

—Katharine Hepburn 

Recent experiments at Thomas Jefferson National Accelerator Facility have mea¬ 
sured the A(Q 2 ) structure function of the deuteron for momentum transfers up to 6 
(GeV/c) 2 [jl]], and measurements for B{Q 2 ) are planned. At such large momentum 
transfers, a relativistic description of the deuteron is required. Even at lower momen¬ 
tum transfers, a relativistic description is important to understand the details of the 
form factors. In addition, incorporating relativity is important for the deuteron wave 
function to transform correctly under boosts to large momentum, which is important 
for calculating form factors. 

One approach that gives a relativistic description of the deuteron is light-front dy¬ 
namics. The subject of this work is to investigate the consequences of combining light 
front dynamics with various nuclear models to calculate bound state wave functions, 
and using them to calculate the deuteron form factors. 

1.1 Theoretical Descriptions of the Deuteron 

In principle, the two-particle deuteron wave function is obtained by solving the full 
quantum held theory that describes the interacting nucleons for the lowest energy 
bound state. In practice, many approximations must be made before even a numeric 
solution is possible. 


2 


There are many approximations that can be used to obtain relativistic two-particle 
wave functions. One approach is to use the Feynman-Schwinger representation (FSR) 
of the two-particle Green’s function 0, which effectively includes all two-particle to 
two-particle ladder and crossed ladder diagrams. However, it requires a path integral 
to be done numerically, and intensive computation is required to obtain an accurate 
answers. While the FSR has been applied to scalar theories, such as the Wick- 
Cutkosky model, it has not been used with more realistic models of the deuteron 
involving spin. 

Another approach comes from using the Bethe-Salpeter equation (BSE) M • 
Solving the BSE is equivalent to solving the FSR with a truncated kernel. This 
simplification allows the BSE to be solved much more quickly than the FSR, however, 
the truncation of the BSE kernel causes the BSE results to differ from the FSR 
results. It is well known that any finite truncation of the kernel yields bound-state 
wave functions with problems, such as the incorrect one-body limit [0. In addition, 
other approximations can be made to the propagators used in the BSE to obtain a 
quasi-potential equation [B-|TTf which addresses some of the problems with the BSE. 

Constructing potentials from a relativistic Hamiltonian theory provides yet an¬ 
other approach. It can be shown that the potentials derived in this manner are phys¬ 
ically equivalent to the Bethe-Salpeter equation. There are several different forms of 
Hamiltonian dynamics [12], each associated with a hypersurface on which the commu¬ 
tation relations for the fields are defined. Equal-time dynamics is the most commonly 
known form of Hamiltonian dynamics, but we will consider light-front dynamics in 
this work, where the orientation of the light-front is fixed. 

There is another light-front dynamics approach that is in use, explicitly covariant 
light-front dynamics [T3|,|T4|] . In that approach, manifest covariance is kept by using a 
null-plane whose orientation is variable. Although that approach has some benefits, 
we are interested in numerical calculations, for which the usual light-front dynamics 
is best suited. 
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1.2 Benefits of Light-front Dynamics 


The utility of the light-front dynamics was first discussed by Dirac Light-front 
dynamics makes use of the light-front coordinate system, where a four-vector x u is 
expressed as x u = (x + ,x~,x 1 ,x 2 ), with x ± = x° ±x 3 . This is simply a change of 
variables, but an especially convenient one. Using this coordinate system and defining 
the commutation relations at equal light-front time (x + = t^p), we obtain a light-front 
Hamiltonian |[T5|-[T7|. We use Hamiltonian in the light-front Schrodinger equation to 
solve for bound states. There are many desirable features of the light-front dynamics 
and the use of light-front coordinates. 

First of all, high-energy experiments are naturally described using light-front coor¬ 
dinates. The wave front of a beam of high-energy particles traveling in the (negative) 
three-direction is defined by a surface where x + is (approximately) constant. Such 
a beam can probe the wave function of a target described in terms of light-front 
variables |15|, 181: the Bjorken x variable used to describe high-energy experiments is 
simply the ratio of the plus momentum of the struck constituent particle to the total 
plus momentum (p + ) of the bound state. 

Secondly, the vacuum for a theory with massive particles can be very simple on the 
light front. This is because all massive particles and anti-particles have positive plus 
momentum, and the total plus momentum is a conserved quantity. Thus, the naive 
vacuum (with p + = 0) is empty, and diagrams that couple to this vacuum are zero. 
This greatly reduces the number of non-trivial light-front time-ordered diagrams. 

Thirdly, the generators of boosts in the one, two, and plus directions are kinematic, 
meaning they are independent of the interaction. Thus, even when the Hamiltonian is 
truncated, the wave functions will transform correctly under boosts. Thus, light-front 
dynamics is useful for describing form factors at high momentum transfers. 

Finally, it is easy to perform relativistic calculations using light-front dynamics. 
This is partly due to the simplicity of the vacuum, and partly due the the fact that, 
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with light-front variables, center-of-momentum variables can be cleanly separated 
from the relative momentum variables. This allows us to write relativistic equations 
which have the simple form of a non-relativistic Schrodinger equation. 


1.3 Breaking of Rotational Invariance 


One serious drawback of light-front dynamics is that rotational invariance is not 
manifest in any light-front Hamiltonian^. This is a result of selecting a particular 
direction in space for the orientation of the light-front. 

An untruncated light-front Hamiltonian will commute with the total relative an¬ 
gular momentum operator, since the total momentum commutes with the relative 
momentum. Thus, eigenstates of the full Hamiltonian will also be eigenstates of the 
angular momentum. However, as mentioned earlier, a Fock-space truncation of the 
light-front Hamiltonian results in the momentum operator four-vector losing covari¬ 
ance under rotations. Hence J 2 and the truncated Hamiltonian do not commute and 
this implies that the eigenstates of the truncated Hamiltonian will not be eigenstates 
of the angular momentum. 

How will this violation of rotational invariance affect physical observables? One 
way to observe this violation is to note that on the light front, rotational invariance 
about the z-axis is maintained. This allows us to classify states as eigenstates of J 3 
with eigenvalues m. We compare the energies of states with the same angular momen¬ 
tum quantum number j but different m values. If the Hamiltonian were rotationally 
invariant, the energies should be the same; the breaking of rotational invariance causes 


the energies to be different [19 


We expect that the higher Fock-space components of the full Hamiltonian will 
be small if the coupling constant is small enough. Thus, truncation at successively 
higher orders in the Fock-space expansion should reduce the violation of rotational 


lr This problem is not unique to light-front dynamics. In any form of Hamiltonian dynamics, at 
least three of the Poincare group generators are dynamic, and thus complicated. 
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invariance of the truncated Hamiltonian. In particular, by retaining enough terms in 
the perturbation expansion of the Hamiltonian, the violation of rotational invariance 
can be reduced to an arbitrarily small amount, provided only that a perturbation 
expansion is valid. Thus, the only real question is: How many terms are required? 


1.4 Outline 


In chapter |2|, the massive Wick-Cutkosky model is used to investigate the degree 
to which rotational invariance is broken for deeply bound states. We use light-front 
dynamics to obtain one-meson-exchange (OME) and two-meson-exchange (TME) po¬ 
tentials, and then we use those potentials to calculate the bound states. We develop 
several methods for quantifying the extent to which rotational invariance is broken 
by the bound states. The binding energies we calculate are compared with those 
calculated using other methods in the literature, such as the Bethe-Salpeter equa¬ 
tion J2^,|2T|] and the Feynman-Schwinger representation of the Greens function p[,|22 


A brief discussion about the approximation we use for the Wick-Cutkosky poten¬ 
tials is in order. It is well known that the vacuum of the full Wick-Cutkosky model 
is unstable [^3j] due to the cubic coupling which provides the interaction. However, 
when the bound-state calculation is restricted to the two-particle sector, the quenched 
approximation is used, and the self-energy and vertex-renormalization diagrams are 
neglected, then the theory has a well defined ground state. In this chapter, we com¬ 
pare the results of our light-front Hamiltonian calculation to the Bethe-Salpeter and 
Feynman-Schwinger representation calculations, both of which use the same approx¬ 
imations. The use of these simplifying approximations allows us to highlight the 
differences between the various approaches. The inclusion of the self-energy diagrams 


and counterterms for the light-front Hamiltonian |2J| and for the Feynman-Schwinger 
representation ||25|-p7f will not be discussed here. 

In chapter II. we introduce a model Lagrangian for nuclear physics which includes 
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chiral symmetry JTB| . The methods introduced in chapter are generalized for use 
with this nuclear model. The Hamiltonian is derived and used to calculate the OME 
and TME potentials for a new light-front nucleon-nucleon potential. We also for¬ 
mulate a new light-front pion-only model. The potentials are used to calculate the 
binding energy of the deuteron. We have some freedom in how to choose the TME 
potentials, and we consider several different choices. 

We are interested mainly in the rotational properties of the potentials and wave 
functions, not in achieving detailed agreement with experimental data. This is be¬ 
cause detailed agreement between calculations in light-front dynamics and experi¬ 
mental data cannot be achieved until the problems of rotational invariance have been 
addressed. From this standpoint, we may simplify our calculations by omitting ef¬ 
fects that are formally rotationally invariant, such at crossed TME potentials, but are 
required for precise calculations. 

The wave functions obtained in chapter |3] are used in chapter f| to calculate the 
electromagnetic and axial form factors of the deuteron. Although rotational invari¬ 
ance demands that there be only three independent matrix elements of the deuteron 
current, the light-front calculation of the deuteron current results in four independent 
matrix elements. This is a result of the lack of manifest rotational invariance on 
the light front. There are several prescriptions for choosing which deuteron current 
matrix element should be eliminated, and in principle this choice will affect the form 
factors. We attempt to find currents that transform correctly (or well enough) under 
rotations so that the choice of “bad” component does not matter too much. 
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Chapter 2 

WICK-CUTKOSKY MODEL 


The Wick-Cutkosky model is an ideal starting point for any investigation of the 
deuteron. The model, introduced with an analytic solution in 1954 by Wick m 


and Cutkosky |29|, is concerned with the interaction of two heavy scalar particles 
mediated by a massless scalar particle. The model can be extended by allowing the 
exchange particle to have a mass. The massive Wick-Cutkosky model cannot be 
solved analytically, but it is amenable to numeric calculations |14,^, 3Q-j35||. 

A large amount of the material presented in this chapter is based on previously 
published work by the author in Refs. 


2.1 Formalism 


We consider an isospin doublet of two uncharged scalars 0 = (0i,0 2 ) with mass M 
(which we will refer to as nucleons ), that couple to a third, uncharged scalar y with 
mass /i (which we will refer to as a meson ) by a 0 2 y interaction. This is denoted as 


a (j) 2 x theory, which is an extension of the Wick-Cutkosky model 


to include 


massive mesons. This model has been used on the light front to study scattering 
states jn| as well as bound states |3B| . The Lagrangian is 


c = \ (^0^0 - M 2 <j> 2 ) + i (0 /t ycCy - // 2 y 2 ) + #y0 2 y, (2.1) 

where g is a dimensionless coupling constant and 0 2 = 0 2 + (j)\. This Lagrangian 
will be used in two formalisms, the Bethe-Salpeter equation and the Hamiltonian 
approach. 






















2.1.1 Light-front Hamiltonian 


To obtain the light-front Hamiltonian from the Lagrangian in Eq. (|2.1|) , we follow 
the approach of by Miller [|H| and many others (see the review in Ref. [[HJ) to write 
the light-front Hamiltonian ( P ~) as the sum of a free, non-interacting part and a 
term containing the interactions. We use the conventions given in Appendix [A]. The 
operators can be expressed in terms of Fock space operators since for this theory in 
light-front dynamics, the physical vacuum is the Fock space vacuum, and thus the 
Hilbert space is simply the Fock space. The Hamiltonian is obtained by using the 
energy-momentum tensor in 

^ J dx~d 2 x± T +fi (x + = 0,x~,x_i_). (2.2) 

The usual relations determine T +M , with 

r „, = + y (2.3) 

“ 9(d|A) 

in which the degrees of freedom (the fields 0 and x) are labeled by 0 r . 

It is worthwhile to consider the limit in which the interactions between the fields 
are removed. This will allow us to define the free Hamiltonian Pf and to display 
the necessary commutation relations. The energy-momentum tensor of the non¬ 
interacting fields is defined as Tf u . Use of Eq. leads to the result 

TT = a'W - - m¥] + W* - 9 — [8„ x d°x - A 2 ], (2.4) 

with 


To = V ± 0 • V ± 0 + M 2 0 2 + V ± y • V ± x + h Z X 


2 ,2 


(2.5) 


The scalar nucleon fields can be expressed in terms of creation and destruction 
operators: 

d 2 k ± dk + 9(k+) 


= 


(2n) 3 P^/2H 


ai(k)e~ ik ' x + <4(fc)e ifc ' a 


( 2 . 6 ) 
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where i = 1, 2 is a particle index, k ■ x — \{k~x + + k + x~) -k^-x^ with k~ = — 
and k = (fc + ,fcj_). Note that k~ is such that the particles are on the mass shell, 
which is a consequence of using a Hamiltonian theory. The 6 function restricts k + to 
positive values. Likewise, the scalar meson field is given by 


= f d^k ± dk+ 0(Jc+) j ^ k y-ik-x + a Uk)e ik ' x ] , (2.7) 

V ' J (2tt) 3 / 2 v / 2F : L X ' J V ' 

2 +fe 2 

where k~ = k+ - , so that the mesons are also on the mass shell. The non-vanishing 
commutation relations are 


Mfc), 4(*0] = S(k ± - k’ ± )8(k + - k ,+ ), (2.8) 

where a — 1, 2, % is a particle index. The commutation relations are defined at equal 
light-front time, x + = 0. It is useful to define 

5 (2 ’ +) (fc - k') = 5(k± - k' L )8(k + - k ,+ ), (2.9) 

which will be used throughout this work. 

We write a ket in the two-distinguishable-particle sector of the Fock space as 

I Ah, k 2 ) = a\(k 1 )al(k 2 )\0). (2.10) 


This implies that the identity operator in this Fock space sector can be written as 


I 2 — / d 2 ki t ±dki / d 2 k 2 ,±dk2\ki, k 2 )(ki, k 2 \. 


( 2 . 11 ) 


The derivatives appearing in the quantity T 0 + are evaluated and then one sets 
x + to 0 to obtain the result 


p o = 


+ _)_ a \(k)a 2 {k )) + r cd x (k)a x (k) 


fi 2 + k 2 ± 


( 2 . 12 ) 


k+ \ V k+ 

with j k = f d 2 k±dk + 6{k + ). Eq. (|2 .1 2|) has the interpretation of an operator that 

r fc + fcj - for the nucleons and M ^ + fc 


counts the light-front energy k (which is M A fcj - for the nucleons and M A |) fc - L for the 


mesons) of all of the particles. 














10 


We now consider the interacting part of the Lagrangian, Cj. An analysis similar 
to that for the non-interacting parts yields the interacting part of the light-front 
Hamiltonian Pf] 


Pf 



x | 2 a\(k + k')a x (k')ai(k) + a x (k + k')ai(k')ai(k) 


Hermitian conjugate 


(2.13) 


The interaction Hamiltonian is self-adjoint since the Hilbert space is the Fock space. 
The total light-front Hamiltonian is given by P~ = P 0 - + gPf. 


2.1.2 Hamiltonian Bound-state Equations 

We will be studying the bound states of two distinguishable nucleons. The technology 
of time-ordered (old-fashioned) perturbation theory is used to construct the light-front 
time-ordered perturbation theory (LFTOPT) for our Hamiltonian. We start with the 
light-front Schrodinger equation in the full Fock space, 

(Po+gPr)\rl> F ) = \1> F )p-, (2.14) 

where P 0 - + gPf is the Hamiltonian in the full Fock-space basis, | ipp) is the wave 
function in the full Fock space, and P - is the light-front energy of that state. Recall 
that Pq“, the non-interacting part of the Hamiltonian, is diagonal in the momentum 
basis, while Pf, which contains the interaction, has only off-diagonal elements. 

A serious drawback of this equation is that the wave function \if; F ) has support 
from infinitely many sectors of the Fock space, since Pf changes the total number of 
particles. However, the components of the wave function with many particles will be 
small compared to the two-particle component if the coupling constant is not too large 
and the exchange particle y is massive. We will construct the two-particle light-front 
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Schrodinger equation which the two-particle component of the wave function satisfies. 
From this construction, we will obtain the rules for the LFTOPT. 

We start by introducing the projection operators V and Q. The operator V 
projects out the sector of Fock space with two distinguishable nucleons and no mesons, 
while Q = I — V projects out all the other sectors. We define 

Wf> = (2.15) 

QlfM = H’q), (2.16) 

so that | i/jp) — If/') + IV'q)- Since the free Hamiltonian does not change the number of 
particles, [P, P 0 ~] = [Qi Pp] = 0- The interaction Hamiltonian changes the particle 
number, so it cannot connect the two-particle sector to itself, thus VPfV = 0. 

Using these projection operators, Eq. (|2.14|) can be broken up into two parts, 

W> + g'PPFQtyo) = I WP-, (2.17) 

(P 0 ” + gQPfQ) \^q) + gQPfV\if>) = \^ Q )P~. (2.18) 

Eliminating the \Fq) and using the expression of the identity given in Eq. ( j2.11| ) we 
obtain the two-particle effective light-front Schrodinger equation 

d 2 pi,±dpf J d 2 p 2 ,±dpt(k 1 ,k 2 \ [P 0 _ + V(d,P~)] |Pi,P2>(PnP2l^> 

= (ki,k 2 \il>)P~, (2.19) 



where P 0 and the potential V act in the two-nucleon basis. The two-particle potential 
is given by 


Q 

V(g, P~) = g 2 VP7 --- P7V. 

; y 1 P~ - P 0 _ - gQPfQ 1 


( 2 . 20 ) 


Note that Eq. ([2.19|) is similar to Eq. ( |2.14| ), except for two main differences. Here 
we have a two-nucleon wave function, which makes it simpler. However, the potential 
depends on the light-front energy P”, thus making it more complicated. 
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The denominator in the definition of the potential is non-diagonal in the full Fock 
space, so the matrix inversion that it represents is highly non-trivial. This problem 
is avoided by expanding the inversion in powers of the coupling constant g to get 


V(g, P~) — VPj 


g 2 Q 


P--P 


E 

0 71=0 


gQ 


P~ - Pn 


P7V. 


( 2 . 21 ) 


This can be simplified further by noting that in the two-nucleon sector of our the¬ 
ory, every meson emitted must be absorbed, so there must be an even number of 
interactions. Thus, the full potential can be written as the sum of n meson exchange 
potentials, 


V(P~,g) — J29 2n V{2n)(P~), (2.22) 

n =1 


where V( 2 n) is the potential due to the exchange of n mesons, given by 


WO = P { p TpS 


Pn 


2n—l 


PfV. 


(2.23) 


To see how to write a sum of diagrams for the potential, we express what what 
Eq. (12.23|) represents in words. We start off with two particles, then the interaction 
occurs. There are two possibilities of what can happen; nucleon 1 or 2 can emit 
a meson. Each possibility has a separate diagram. After the interaction, there is 
propagation with the light-front Green’s function, 


Glf(P-) = p _ ]_ . (2.24) 

until another interaction occurs, and so on. We simply sum up all of the possible 
orderings of the interaction to get the full potential. The n th order potential is simply 
the sum of all possible diagrams with n-meson exchanges. 

Each intermediate state in Eq. ( |2.23| ) has more than two particles, so the diagrams 
are two-particle irreducible with respect to the two-particle Green’s function G 2 lf = 
'PG'lf'P- We can represent G 2 lf by its diagonal matrix elements, 


G 2LF (k 1 ,k 2 ;P ) - p _ k _ _ k 


( 2 . 25 ) 
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In the diagrams we draw, the nucleons will be represented by solid lines and the 
mesons by dashed lines. Although the states we will be considering consist of two 
distinguishable nucleons, we will not label the nucleon lines. Energy denominator 
terms are represented by vertical, thin, dotted lines. We will be using the quenched 
approximation (so there are no nucleon loops) and neglect the mass and vertex renor¬ 
malization diagrams (so the physical masses and coupling constant are used, and each 
meson emitted from one nucleon must be absorbed by the other nucleon). It is not 
expected that these restrictions will lead to qualitatively different results than the 
true full solution when the states are not too deeply bound. The quenched approxi¬ 
mation is reasonable when the masses of the nucleon fields are large compared to the 
binding energy. Use of the physical masses and coupling constant are reasonable as 
well when the momenta are not too large. 

A truncation must be made for the potential in Eq. (|2.22| ), since in a Hamiltonian 
theory the infinite sum of graphs cannot be calculated. However, all the graphs are 
required to obtain manifest Lorentz invariance. We stress again that we will compare 
various truncations of the light-front Hamiltonian to other calculations which do not 
include renormalization diagrams. This is because we want to determine the effect of 
truncation on the light-front Hamiltonian. The differences between our calculation 
and those which include the self-energy graphs [^,^], which may be large for deeply- 
bound states, are not considered here. 

The rules for drawing the n-meson-exchange graphs that correspond to this ap¬ 
proximation are: 

1. Draw all topologically distinct time-ordered diagrams with n mesons. Use solid 
lines for the nucleons and dashed lines for the mesons. 

2. Delete all graphs which couple particles to the vacuum. In the massive theory 
we consider here, these diagrams always vanish since the vacuum has zero plus 
momentum and massive particles always have positive plus momentum. 



14 


3. Our quenched approximation and use of the physical masses and coupling con¬ 
stant requires us to delete all graphs that have nucleon loops or have mesons 
that are emitted and absorbed from the same nucleon. 


Once the diagrams are drawn, we use the following rules to convert the sum of 
diagrams into the potential (kif, k2f\V(2 n )(P~)\ku, k 2 i): 

1. Overall factor of -— hfci/+fc 2 / fc n frO . This delta function says that the to- 

2(27T)3 y/k+k+k+k* 

tal light-front three-momentum is conserved. We define the light-front three- 
momentum P = k\f + fe 2 /- 


2. To each internal line, assign a light-front three-momentum q t where 1 <i < N 
and N is the number of internal lines. The light-front energy for particle i with 
mass rn t is qi = 1 + I,X . It is useful to define z t = q[ / P + . 


3. 

4. 

5. 


O(z') 

A factor of T ; for each internal line. 

Z T 

An extra factor of p A +p- for each internal meson line. 


A factor of , _ p —rr between consecutive vertices, where the sum is over only 

(T Ai Qi ) 

the particles that exist in the intermediate time between those vertices. 


6. Use light-front three-momentum conservation to eliminate all the independent 
momenta. 

7. Integrate with f 2 p+ q p-(2n) 3 over remaining free internal momenta. 

8. Symmetry factor of | when two nucleons are created or destroyed at the same 
time. 


With these rules, one can calculate the effective potential for any order. 
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2.1.3 Further Development of the Light-front Schro dinger Equation 


Once the potential is calculated, we can plug it into Eq. (|2.19|) , which we write as 


J d 2 k u^dkfi J d 2 k 2 i,±dk^ i (k 1 f,k 2 f\[P 0 + V(g(P ),P )] \k u , k 2i )(k u , k 2i \if) 

= (kif,k 2f \^)P~, (2.26) 


where P~ is an arbitrary light-front energy and g{P~) is the coupling constant which 
yields the bound-state wave function with P~ as the bound-state energy. We call 
this g(P~) the spectrum of the light-front Schrodinger equation for the corresponding 
wave function. 

The total momentum P = kif + fc 2 / is conserved by the potential given in 
Eq. ( [Egg ), so the wave function in Eq. ( |2.26| ) can be parameterized by the total 
momentum. To make the calculations easier later, we choose to be in the center-of- 
momentum frame, where the components of the total momentum can be written as 
Pi = 0 and P + = P~ = E. The energy, E, is the same as the mass of the bound 
state. In terms of the binding energy B, E = 2 M — B. In the center-of-momentum 
frame, the wave function is parameterized by E, so we can define 


(ki f , k 2f \if M ) = S^ + \k lf + k 2f - PMM, (2.27) 

(k\f, k 2 f\V(2n){P )\ku, k 2 f) = 8^ 2,+ \kif + k 2 f — ku — k 2 i)V( 2n ')(E ; kif ; k\f). (2.28) 

With these, Eq. ( p.26| ) effectively becomes a one-particle equation, where particle 2’s 
momentum is determined by k 2 f = P — kif. The minus component (the light-front 
energy) of particle 2 is defined by the requirement that the particle 2 is on mass shell, 
so kff = (M 2 + k\f L )/ky. We also define 

x = kff/P + = XBj, (2.29) 

where XBj is the Bjorken x variable, so that ky/P + — 1 — x. Likewise, we write 
the Bjorken variables that correspond to the momenta ku and q v as y = K/P + and 
z = qf/P + . 
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Using Eqs. ( |2.28| ), and the fact that the plus momentum of both nucleons 

is positive, we can write the light-front Schrodinger equation Eq. ( |2.26|) as 


J d 2 k lit± J dkf i V(g(E),E]ki f ]k li )'ip(k li ) = ^(k lf )(E » k lf - k 2f ). (2.30) 

It is useful to convert from light-front coordinates ki = (kf, k±) to “equal-time” 
coordinates fc ET = (fcj_, k 3 ), using an implicit definition of k 3 |39| 

E 


K = 


2/c°(fc E T 
k°(k ET ) = sJlVP + k 2 ET , 


[k°(k ET ) + k 3 ] , 


(2.31) 

(2.32) 


or an explicit transformation 


k° 


M 2 + k 2 ± 
4a; (1 — x) 1 


k 3 = (kf — k 2 ")— = (2x — l)k°. 
E 


(2.33) 

(2.34) 


Often the explicit dependence of k° on fc E T will not be shown. It is worth em¬ 
phasizing that this is just a convenient change of variables; ij){k E t) is not the usual 
equal-time wave function. With this transformation, we can express fcf and k 7 as 


k 3 \ 2 

kl ~ 

k 3 \ 2 (k o y 


kt> ~ ' 1 + jfco 


E 


(2.35) 

(2.36) 


We may also write 


*-< ‘4)t 

k 3 \ E 


H = 11 


k° 


(2.37) 

(2.38) 


Using Eqs. ( |2.35[ - |2.38l ), we can rewrite Eq. ( p.30|) as 
2 ktkt 


d 3 ki 


ET- 


k? 


V(g(E), E, fc /;ET ; fe iiET )^(fei, et) = V#/,et) [-E 2 - (2 k° f ) 2 ] . (2.39) 
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Now consider the exchange of the particle labels 1 and 2. This causes 


ki,± —>► k 2) ± = — fepj., 

(2.40) 

kf —> H = E — kf, 

(2.41) 


which means that k 3 as defined in Eq. Q2.31|) transforms as k 3 —> — k 3 , so k E t —> — ^et- 
Consequently, exchange of particle labels 1 and 2 is the same as parity for fc ET . 

For scattering states, the OBE potential computed for on-shell nucleons and used 
in the Weinberg integral equation |4l[ (which is essentially the scattering analogue of 
Eq. ( 2.26 )) leads to manifestly rotationally invariant results when written in terms of 
the relative momenta ]18|. The similarity between that rotationally invariant result 
and the usual equal-time result implies that this equal-time momentum k can be in¬ 
terpreted as the relative momentum of the two particles. For bound states, this exact 
simplification does not occur, since for bound states the potential is, of necessity, 
evaluated off the energy shell (but on the mass shell). However, we expect that the 
OBE potential, written in terms of the relative momentum, is approximately spheri¬ 
cally symmetric for lightly-bound states. Thus, the wave functions are approximate 
eigenfunctions of the “relative angular momentum”, which we define by the operator 
L = x x k. Our “relative angular momentum operator” is not the same as the true 
orbital angular momentum operator which is obtained from the Lagrangian via the 
energy-momentum tensor in a way similar to the Hamiltonian. 


Since the two nucleons are identical except for the particle label, the effective 
potential commutes with parity to all orders in g 2 . Furthermore, the light-front 
Hamiltonian is explicitly invariant under rotations about the three-axis. These con¬ 
siderations allow us to classify the wave functions as having eigenvalues p of parity 
(V) and m of the three-component of the angular momentum operator (J 3 ). We label 







18 


the wave function as (/ci|t/Vp); where 

<*i|J 3 |C, P ) = (2-42) 

(ki\V\ C, P ) = (felO* ( 2 -43) 

= (^i|V4,p)P- (2-44) 

With this, we may write the OBE truncation of full uncrossed Hamiltonian as 

^obe(#> E) = P 0 + g 2 V OBE (E), (2.45) 

which gives OBE light-front Schrodinger equation 

foes (sSKf®, e) |«° BE } = |*™ e )-B, (2,46) 

where E is an arbitrary energy, |^p BE ) is the n th wave function with parity p and 
J 3 quantum number m, and , 9 obe P (E) is the coupling constant which yields that 
bound-state wave function with E as the bound-state energy. 

For the OBE+TBE truncation, we have 

P~ BE (g, E ) = P 0 ~ + g 2 V OBE (E) + g 4 V TBE (E), (2.47) 


which gives TBE light-front Schrodinger equation 


Ptbe(&be P (E),E 


I / n,TBE\ _ | / n,TBE\ j? 

\Vm,v / “ IWm.p )E, 


(2.48) 


where the quantities here are defined in an analogous way to Eq. (|2.46| ) . By comparing 
the spectra go BE P (E) and g^^f(E), we can see what effect adding the TBE potential 
to the OBE potential has on the coupling constant for a given bound-state energy. 

The quantum numbers m and p of the wave function can be used to rewrite 
Eq. ( p.39p as 


dE 


'i,ET 


r ‘ 7r / 2 2kt i kt i k‘f Er r sindj 

dOi - -jjj - V Pt m{kf } ET, Of] &i,ETj Qf)lpp,m(ki, ET, Oi 


'0 


— 4>p,m(kf, et, 6f) [ E 2 — 4(fc°) 2 ] , 


(2.49) 
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where 


Vp t m(kf t ET, Of] ki t ET, Oi) 

= r Vm(kf, ET, Of-, ki^T, Oi) + pV m (kf : ET, Of] fcj,ET, ^ ~ Oi) 


Vm{kf, ET, Of] k if ET, Oi) 


2n d<fif d(f>i 


2vr 


e im(^-^ )y ( fc/jET;fc . ET )_ 


(2.50) 


(2.51) 


We call V (kf t ET, Of] k^ et, #i) the azimuthal-angle-integrated potential. The light-front 
vectors can also be expressed in terms of the azimuthal angle (j), the angle between k± 
and the :r-axis. This allows the azimuthal-angle-integrated potential to be written in 
light-front coordinates 


bm(A if, k\f } j_, kn, kn t ±) 


2t r 


d(j)f d(j>i 

2 t r 


e im(0 / -^)y( fei/;fei .)_ 


(2.52) 


All of the simplifications of Eq. ( 2.49 ) based on physical considerations have been 
addressed. However, further rearrangements need to be done before Eq. fl2.49|) is fit 
to be solved on the computer. Since these involve only numerical techniques, they 
are relegated to Appendix |B|. 


2.1.4 Bethe-Salpeter Equation 

The Bethe-Salpeter equation ||-|7j] provides a way of describing bound states based on 
Feynman propagators and kernels constructed from covariant quantum held theory, 
and thus is manifestly covariant. The equation for the bound state of nucleons 1 and 
2 can be written as 


GKi) = 1 ), (2.53) 

G is the free two-particle propagator, which is the product of two one-particle prop¬ 
agators, G = S 1 S 2 , 1) is the four-dimensional Bethe-Salpeter amplitude, and K is 


the sum of all two-particle irreducible two-to-two Feynman graphs, shown in Fig. 2T 
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t- -a- r - 1 t 1 -r 





Figure 2.1: The full kernel for the Bethe-Salpeter equation of the Wick-Cutkosky 
model. 


—T 

^ladder ^ 

l 


Figure 2.2: The ladder (one-boson-exchange) kernel for the Bethe-Salpeter equation. 


There are no nucleon exchange graphs since we treat only the case where the two 
nucleons are distinguishable. 

We consider the ladder approximation to the Bethe-Salpeter equation, 

/GA'ladderV’ = V’, (2.54) 

which is obtained by replacing K in Eq. fl2.53[) with g 2 Aj a dder- We define A'i a dder = 
~2 Ai a dder > where A'iadder is the graph due to one-boson-exchange, shown in Fig [2~^ . 
This definition of A'i a dder makes it independent of the coupling constant. Making 
this approximation leaves the Bethe-Salpeter equation covariant. This implies that 
the Bethe-Salpeter amplitudes ip have definite angular momentum l, and hence the 
energies of the states are degenerate for different m projections of the same angular 
momentum. 

This equation can be simplified in the center-of-momentum frame. In that frame, 
once the total energy is defined as P°, the four-momentum of the second particle is 
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Figure 2.3: The crossed (two-boson-exchange) kernel for the Bethe-Salpeter equation. 


given by — fcf, and thus the Bethe-Salpeter equation is effectively a one- 

particle equation. The remainder of the discussion in this section will be done using 
the center-of-momentum frame. We then can rewrite Eq. ( p.54| ), taking into account 
the explicit symmetries of the equation, as 


/ lbse (-^) G(E)Ki adAer (E)il) n j trn — 


1 2 


(2.55) 


where E is an arbitrary energy, Gn,z,m is the n th Bethe-Salpeter amplitude with an¬ 
gular momentum l and three-projection m, and <7 lbse(-^) * s the coupling constant 
which yields that bound-state Bethe-Salpeter amplitude with E as the bound-state 
energy. We call this g(E) the spectrum of the ladder Bethe-Salpeter equation for the 
corresponding Bethe-Salpeter amplitude. The Greens function G(E) and the ladder 
kernel /Gadder (A) are functions of the energy in the c.m. frame and are implicitly 
effective one-particle operators. 

We may also take the kernel to be K = /Gadder + /Grossed, where /Grossed = 
■jzK crossed, and /Grossed is the crossed diagram is shown in Fig. p~3| . Since this kernel 
is also covariant, the discussion given above for the ladder kernel also applies for the 
ladder plus crossed kernel. In particular, the bound states calculated by using kernel 
with the Bethe-Salpeter equation will be angular momentum eigenstates. 

The comparison we draw is only between the ladder (or ladder plus crossed) Bethe- 
Salpeter equation and the light-front Hamiltonian that corresponds to that approxi¬ 


mation. The exact nature of the correspondence is discussed in section 2.1.2. Since 
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we are not looking at the solutions to the full theory, for our purposes it does not mat¬ 
ter that there are sizeable differences between the solutions to the full Bethe-Salpeter 
equation and the ladder approximation when the coupling constant is large p2| , p^ -|35| . 


2.1.5 Comparison Between the Light-front Schrodinger Equation and the Bethe- 
Salpeter Equation 

The solution to the untruncated light-front Schrodinger equation in the uncrossed 
approximation is equivalent to the solution of the ladder Bethe-Salpeter equation. 
When the full uncrossed Hamiltonian is truncated, differences will be introduced. 
Thus, here we think of the ladder Bethe-Salpeter equation as the exact theory which 
the truncated light-front Schrodinger equations approximate. As more graphs are 
included in the truncated light-front Hamiltonian potential, the agreement with the 
Bethe-Salpeter equation will obviously be better. The question we wish to answer is 
how well the spectra g(E) for the two different truncations [Eqs. (|2.46| ) and (|2.48| )1 
approximate the spectra for the “exact” theory, the BSE results. 

The lack of manifest rotational invariance of the truncated Hamiltonian the¬ 
ory causes a breaking of the degeneracy of the spectra of the truncated light-front 
Schrodinger equations for different m states — unlike the case for the Bethe-Salpeter 
equation. The wave functions from the Hamiltonian approach are classified by their 
dominant angular momentum contribution l, so that we can compare the spectra for 
different m projections of the same total angular momentum l. By doing this, we can 
compare how the degeneracy of the spectra is broken for the OBE and the OBE+TBE 
truncations, and also compare to the spectra obtained from the ladder Bethe-Salpeter 
equation. 

We want to approximate our potential V so that Eq. (|2.26|) is physically equivalent 
to the ladder Bethe-Salpeter equation. This approximation of V will be called the un¬ 
crossed approximation. By physically equivalent, we mean that the spectra of the po¬ 
tential V should reproduce the spectrum for the states of the Bethe-Salpeter equation, 
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excluding the so-called “abnormal” states PBI,ED|. It is well known how to reduce the 


Bethe-Salpeter equation to a physically equivalent Hamiltonian (Schrodinger-type) 
equation. For an extensive discussion of this issue of defining the potential equivalent 
to a sum of Feynman graphs in the equal-time case see, for instance, Klein [Q, Phillips 
and Wallace |43|], Lahiff and Afnan [[h|], and, for examples on the light front, Ligterink 
and Bakker [45j and Schoonderwoerd and Bakker [46]. The general procedure to get 
the effective potential due to n boson exchange takes two steps. First, write the sum 
of all Feynman graphs obtained from iteration of the Bethe-Salpeter kernel with n 
boson exchanges. Then, write that sum in terms of LFTOPT graphs, and discard all 
graphs which are not two-particle-irreducible with respect to the light-front two-body 
propagator, 


Glf (P-) = p _ * p _ . (2.56) 

The graphs which remain after this procedure constitute the Wbe- 

As an example, we construct the TBE potential. When the ladder Bethe-Salpeter 
equation is used, only one Feynman graph contributes, the box diagram arising from 
the iteration of the Feynman OBE kernel. This gives six non-vanishing LFTOPT 
diagrams. (The other diagrams vanish because the vacuum is simple on the light 
front.) 

The first four diagrams are iterations of the OBE potential and are reducible with 
respect to Glf- The last two are two-particle-irreducible and thus constitute the TBE 
potential, Vtbe- 


2.2 XVick-Cutkosky Perturbative Potentials 

A truncation must be made of the expansion of the potential given in Eq. ( |2.22| ), since 
it is not feasible to calculate the infinite sum of graphs for the potential in this Hamil¬ 
tonian theory. For this investigation we consider three truncations of the potential 
derived from the field theory. First, the OBE potential and the TBE potentials are 
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Figure 2.4: The relation between the Feynman box diagram and non-vanishing light- 
front time-ordered perturbation theory diagrams. 


calculated. Then, we note that a subset of the TBE diagrams, the stretched-box dia¬ 
grams, correspond to the truncated potential derived from the ladder Bethe-Salpeter 
equation. Thus, three truncated potentials are obtained that have a physical inter¬ 
pretation. 

The matrix elements of these potentials are written in the two-particle momentum 
basis, denoting the momentum of the incoming particles by ku and k 2i , and the 
outgoing particles k\f and k 2 f. For simplicity, we choose to work in the center-of- 
momentum frame. By inspecting the rules for converting a light-front time-ordered 
diagram into a potential given in section EE3. and looking at Eq. (|2.28|) , we find that 
each piece of the effective-one-particle potential V(E; k±f, ku) is proportional to 


E 


2(2tt) 3 


U+ U+ £.+ £.+ 


(2.57) 


(The factor of E in the numerator is included to simplify later equations.) This term 
will by suppressed in all of the potentials written in this chapter. 


2.2.1 0BE Potential 

We start by drawing all the allowed and non-vanishing time-ordered diagrams with 


one meson exchange. These diagrams are shown in Fig. 2J3. The light-front time- 
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Figure 2.5: The diagrams for the OBE potential. 


ordered perturbation theory rules given in section 2.1.2 are used to calculate the 


potential due to OBE potential, 


Robe {E\ kif ; bn) 


{MV 6(x-y)/\x-y\ 

\EJ E - k^-k^ f - uj~(k lf - k u ) 

_i_ Q(y - x)/\y - x\ 

E - k~ f - k^i - u-{k u - ki f ) 


(2.58) 


We have introduced the notation that meson with light-front three-momentum q has 
a light-front energy given by 


u; (q) 


V + q± 



(2.59) 


The azimuthal-angle integration of Vobe is discussed in Appendix |C.1| . 

The potential given in Eq. ( |2.58| ) can also be used for scattering states. In that 
case, E = /p + /p = V + V , which allows the potential to be written as 

FoBE (-Escat; k\'. P]') = JJ-S (2.60) 


The scattering potential is the same as the usual equal-time OBE potential. This 
must be the case, since the scattering potential is also given by covariant Feynman 
diagrams, which have the same form independent of the form of dynamics. 

Returning to the bound-state regime, we note that the OBE potential can be 
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written in terms of the equal-time coordinates. A reorganization of Eq. 


yields 


Vobe (E; kf t et; ki et) 
2 


M 

~E 


E 


0(x - y) 


(K f ~ - k~ f ) -y 2 - 0 k i± - k f± y 

_ Q(y- x) _ 

(k\i ~ kt f )(E - k~ f - k. 2i ) - [E - (fcu. - k fl _)‘- 

Using the relations in Eqs. (|2.35| - |2.38|) , we hnd 

' E 2 - AM 2 


(2.61) 


k 3 k 3 

rv j rvj 


k 2 - k 2 

\ ET ft '/,ET 


+|| (*? - *?) 2 - (d - *, 3 ) 2 . 

Ivf Atj 


(2.62) 


Under the exchange kf ki , the only thing that changes in Eq. ( j2.62| ) is that the 
first term picks up a minus sign. This observation allows Eq. (|2.61|) to be rewritten 
as 

E 

E J ~ *~ 


'M\ 2 

U OB E(-E;fc/,ET;fci,ET) = ( — ) -j^3— 


pJ TtT 

'v K i 


k3 f k ? / .0 \2 


A + 


Qe t - h 2 


(2.63) 


where 


A _ A 2 - AM 2 l2 l2 

^ — 2 k iET kf ET: 


Qet ~ kf t et k? 


ET- 


(2.64) 

(2.65) 


Note that q ET is not the light-front energy of the meson, since in a Hamiltonian 
theory only the light-front three-momenta are conserved; the four-momenta are not 
conserved. Equation (|2.63|) will be useful in the context of approximations based on 


the physical arguments that we will discuss in section 2.3 


2.2.2 TBE Potentials 


As in the previous section, we start by drawing all the allowed, non-vanishing time- 


ordered diagrams with two meson exchanges shown in Fig. 2.6. The diagrams are 






















27 


classified according to the behavior of the intermediate particles. The total TBE 
potential is given by the sum of all the diagrams, so 


E, 


TBE 


= V , 


TBE:SB 


En 


TBE:SX 


+ E 


TBE:TX 


E 


TBE:WX 


E 


TBE:ZX- 


( 2 . 66 ) 


In the diagrams for the TBE potential in Fig. |2.6| , the intermediate loop momenta 
can be parameterized by ki m or k 2m . The dependent variable is defined by the relation 
P = ku + fc 2 j. The Bjorken x variable that corresponds to fe lm (fc 2m ) is labeled with 
z (1 — z). We use the Feynman rules to calculate all of these potentials, starting with 


1 / TBE:Sb(£'; kif] fell) 

/M \ 4 r d 2 k m ± T / -i 9(z - y)6(x - z) 

\EJ J 2(2tt) 3 [J 0 Z z(l - z)(z - y)(x - z) 


1 

E - k^ m - k~ f - u~(k lf - fe lm ) 

1 

E k^i U) (fel f felm) W (felm fell) 

1 

E-k^- k7 lm - a;-(feim - fei*)_ 


+ -(1 > 2}. 


(2.67) 


The symbol {1 *-*■ 2} means that all labels 1 are replaced with 2 and vice versa, as 
well as replacing the Bjorken variables x, y, and z with 1 — x, 1 — y, and 1 — z. This is 
a way of explicitly stating the symmetry of the potential under exchange of particles 
1 and 2. A detailed discussion of the evaluation of the loop integral in Eq. (|2.67| ) is 
given in Appendix |C.2. 












(e) V TBE:ZX — 



+ 



Figure 2.6: The components of the TBE potential, (a) the stretched-box, (b) 
stretched-crossed, (c) T-crossed, (d) wide-crossed, and (e) Z-crossed diagrams. Here, 
kirn — k i f T k\i ki m and /C 2 m — ^ 2 / T k-h k 2 m - 
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It is straightforward to calculate the other parts of the TBE potential, 


VtEE:Sx(E; fcl/i fell) 
4 /* 


M 

~E 


ra_L 


2(2vr) 




0(x — z)0(2 — y) 


X- 


z(x — z)(l + z — x — y)(z — y) 

1 


E — k 


1 m 


k 2 f u> {k\f 


k 


Ira ) 


l 

_ 

E k 2 f U) (fel / felm) ^ {k\ m fell) 

1 

X-—- 

E - k^ - k 2m - u~(k lf - k lm ) 


+ {!<-> 2}, 


( 2 . 68 ) 


where we have denoted the light-front energy of particle 2 by k 2m , given by 

(2.69) 

(2.70) 

The rest of the TBE potential is given by 


^2 m — e 


(P + k 


Ira 


k\i — k\f ), 


■w - ^ 


I'tbe:Tx(E; kif] ku) 

/M \ 4 r d 2 k ml _ r r 1 6{x - z)9(l + z - x - y)9(y - z) 
= \E ) J 2(2tt) 3 L/o * z(x -z)(l + z-x- y)(y - z) 

1 

_ 

E k^ m k 2 j u {k\f ^im) 

1 


x--- 

E k^ m k 2m U1 (^1/ ^lm) ^0 (kli felm) 

1 

X--- 

E - Ki - Km - u~(k lf - k lm ) 


+ {!<-> 2}, 


(2.71) 



30 


VtBE:Wx(-E; kif] fell) 

M\ 4 f d 2 k m ± f l 6{x - z)6{ 1 + z-x- y)0(y - z) 
E ) J 2(2tt) 3 [J 0 Z z(x - z)(l +z - x - y){y - z) 

1 


x 


E k^ m k 2 j ^ (fei f feim) 

1 


x ■ 


E - k lm - k 2m - u~(k lf - k lm ) - uj-{k u - k lm ) 
1 


x- 


E - k lm - k 2i - u (k u - ki m ) 


+ {!<-> 2}, 


(2.72) 


VtBE:Zx(-E; kif] kli) 

(mV fd 2 k m± 

\E J J 2(2vr)3 


-l 


dz- 


9(x + y — 1 — z) 


x 


x 


x 


Uo z(x-z)(x + y- l-z)(y-z) 

1 

E - k^ m - k 2f - iv~{k lf - k lm ) 

1 

E - V m - k 2f - k 2i - e~ (feu + feu - P - k 


1 m ) 


E - k lm — k 2i — uj {ku ~ ki m ). 


+ {!<—> 2}. 


(2.73) 


The loop integrals in the expressions for the TBE potentials and the azimuthal-angle 
integrations are discussed in Appendix |C.2| . 


2.2.3 TBE:SB Potential: Connection to the Ladder Bethe-Salpeter Equation 


As discussed in Section [2.1.5| , a truncated kernel of the Bethe-Salpeter equation is 
physically equivalent to a Hamiltonian potential which does not include all the graphs 
that the full theory allows. In particular, consider the Bethe-Salpeter equation when 
the ladder kernel is used. The physically equivalent light-front potential will not 
include any graphs where the meson lines cross, so to order g 4 , the potential is given 
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by 5 ,2 Vobe + 5 ,4 Vtbe:SB- Therefore, by considering the TBE:SB truncation, we can test 
how well the light-front Hamiltonian approach approximates the full ladder Bethe- 
Salpeter equation. This concept is discussed more throughly in Refs. |36],[I^,|45| . 


2.3 Wick-Cutkosky Non-perturbative Potentials 

The potentials discussed in this section are derived from the OBE field theory poten¬ 
tial, but additional approximations are made to simplify the expressions. 


2.3.1 Symmetrized-mass Approximation 


Krautgartner, Pauli and Wolz ^7|, and Trittmann and Pauli |19|] studied positronium 
with a large coupling constant in light-front dynamics. The one-photon-exchange 
potential they obtain has a colinear singularity due to the sum of the instantaneous 
photon exchange graph and a gauge-dependent factor from the spin sum. They argue 
that the singularity is not physical, and therefore must be canceled by higher-order 
terms in the potential. The effect of those terms can be simulated by choosing the 
bound-state energy so that the coefficient of the singular term vanishes. They find 
that this condition is met when the light-front energy P~ in the one-photon-exchange 
potential is replaced with the operator u, expressed here in the two-particle basis, 


P- 


v {kif, k 2 f, k u , k 2i ) = - (k u + k 2i + k lf + k 2f ) 


(2.74) 


2 V-1* ' -a ' • v lf ' 2 f) 

This approximation is called the symmetrized mass [IT?] , the average of the total P~ 
in the initial and final states. It is important to note that this approximation affects 
not only the singular term, but also the energy denominator in the rest of the OBE 
potential. The modified denominators simulate the effects of the non-perturbative 
higher-order terms that are not included explicitly in the OBE potential. Potentials 
obtained with this approximation are similar to those given by the unitary trans¬ 
formation method ]4^,|4D|, where the potentials depend explicitly on the initial- and 
final-state energies. 
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In our model, there are no singularities associated with the OBE graphs because 
we deal only with scalar fields. However, we may use their approximation to obtain a 
new light-front OBE potential that should incorporate some non-perturbative effects. 
Recalling that the only place where P~ occurred in Eq. ( |2.58|) was in the denominator, 
the E in the denominator of the OBE potential is replaced with u; to get 


k\ i) 

'mV 

~e) 


9{x - y)/\x - y | 


I {~ k u + k 2i + k if~ k 2f) ~ w “( fc i/ - k ii) 
0(y ~ x)/\y - x\ 


\ i +k 


H k 2i 


k !f + k 2 /) ^ ( k li k l /) 


M 

~E 


E 


(2.75) 

(2.76) 


H k tf ~ k ti) {~ k u + k 2i + k if~ k 2f) - V - (kf± ~ k i± ) 2 ' 

Writing the light-front variables in the denominator in terms of the equal-time vari¬ 
ables, as prescribed in Eqs. ( |2.35| - [2.38D , we find 


1 k 3 k 3 

(-*r,++ k, - Ki) = -(*?- kn 2 + 


{k° f ~ fc°) 2 - (2-77) 


Thus, Eq. (|2.76| ) can be rewritten as 


V UJ (kf t ET', k 


ET, K'i.ETj — 


M 

~E 


E 


k :f H , 

TJT TJJ \ 

Kf K i 


(2.78) 


$(k°f - k i ) 2 - (fc/,ET - fci.ET) 2 - /i 2 

This result can also be obtained more directly by considering the Erst term in 
Eq. (|2.62D . Recall that the E 2 that appears in the denominator is written as P + P~ 
in an arbitrary frame, so in the symmetrized-mass approximation, the E 2 term is 


replaced with Euj. This causes the A term in the denominator of Eq. ( 2.63 ) to 
vanish, so that the equation reduces to Eq. ( |2.78| ). Also, note that by writing this 
new potential, we attempt to incorporate physics from higher-order graphs than just 
the OBE graphs. 


The singularity structure of the symmetrized-mass potential is easily analyzed. 
When scattering states are used, in the center-of-momentum frame the total energy 















of the state is E — 2 k® = 2 k®, so the relations in Eqs. (|2.35| - |2.38|) become 

kf = k°± k 3 , 
kf = k° -F k 3 . 
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(2.79) 

(2.80) 

Using these relations, the symmetrized mass is a ) — E. Tims, for scattering states, 
this potential is same as the OBE scattering potential and the singularity structure 
is the same. 


2.3.2 Instantaneous and Retarded Approximations 


For our bound states, k ET -C M 2 , so that the — k®) 2 will be much smaller than 

the other terms in the denominator of Eq. (j2.78p . Therefore, we may approximate 
the symmetrized-mass potential Vb by the instantaneous potential, 

f M\ 2 -E 

hlnst(fc/,ET; &?:,Et) = ( — J 77 ;-y- , 9 - (2-81) 


\ E J (kf, et — k it et ) 2 + h 2 
Alternatively, we may also argue that since the energy difference term is small, we 
can also approximate V u by the retarded potential, 


hpiet(fc/,ET; k 


ET, K'i,ETJ ~ 


M 

~E 


MV 


E 


(*? - - <*/, et 

E 


k / 


i,ET 


) 2 


k 


(2.82) 

(2.83) 


E J (kf t et — ^7,et) 2 — h 2 ’ 
where kf t et and k it et represent four-vectors, defined by the equal-time three-vectors 
and the condition that A; 2 ET = kf ET = M 2 . These potentials resemble the three- 
dimensional Blankenbecler-Sugar [II]] or Gross [11| quasi-potentials. 

Both of these approximations are reasonable if the energy difference between the 
initial and final states is small, which is valid for lightly-bound states. The instanta¬ 
neous potential is a better approximation of the symmetrized-mass potential, since, 
if we expand the symmetrized-mass potential to second-order in perturbation theory 
about k® = A;°, we get V u = Unst- Also, note that these potentials are explicitly 













34 


rotationally invariant in terms of our equal-time parameterization, which provides 
significant computational advantages. 


2.3.3 Three-dimensional Reduction of the Bethe-Salpeter Equation 


We now consider a non-perturbative approximation used by Wallace and Mandelzweig 
|50],[5l]]. The basic idea is to first make an approximation of the Bethe-Salpeter equa¬ 
tion, then reduce that modified Bethe-Salpeter equation to the physically equivalent 
Hamiltonian equation. This approach was used by Phillips and Wallace [[IIJ for the 
model we use, however, they obtained an equal-time Hamiltonian, while we seek a 
light-front Hamiltonian. Before we do this, we first review the basic mechanics of the 


three-dimensional reduction, as presented in Sales et al. |38j and specialized to our 
particular case. We postpone the discussion of the approximation until section |2.3.4 . 

The Bethe-Salpeter equation can be written in matrix form as T = KGoT, or 
explicitly in function form in the momentum basis as 

d 4 ku 


n*if,p) = 


(2tt) 


K(k lf ,k li -P)G 0 (k li -P)r(k li -,P). 


(2.84) 


Here, T is the four-dimensional vertex function, K is the four-dimensional kernel, 
and Go is the two-particle four-dimensional Green’s function. The momenta are P, 
the total four-momentum, and k\, the four-momentum of particle 1. Particle 2’s 
momentum is implicitly P — k\. The four-dimensional Green’s function is given by 


G 0 (ki, P) = idik^d^P-k i), 


(2.85) 


where d is the one-particle Green’s function. On the light front, d can be written as 

1 \ 1 


dip) = 


p+J p — Sign(p+)e ( p )’ 
where the light-front energy e~ is given by 

^ M* + pl . 

e (p)= nn— 


( 2 . 86 ) 


(2.87) 
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The real part of e is a positive definite quantity, and rj is positive infinitesimal. 

as 


The Bethe-Salpeter equation can be rewritten [52 


T = WG 0 T, 


where G$ is an auxiliary Green’s function, and W is defined by 


( 2 . 88 ) 


W = I< + K(G 0 - G 0 ) W. 


(2.89) 


The advantage of this rearrangement is that we are free to choose the form of the 

is 


auxiliary Green’s function, G 0 . The choice of G o advocated in Ref. |38 


- h (2 ’ + )ffci ^ ) 

G 0 (k lf , k u - P) = G 0 (k lf ; P) . m g o(fcih P), (2.90) 


9o(k lf ;P) 


where 


9o{ki,P) = 


dk x 

2 ( 2 ? 


■Go(^i; P ). 


(2.91) 


There is an extra factor of 2 in the denominator of Eq. ( |2.91 ) when compared to the 
equal-time formalism. This is due to the Jacobian of the light-front coordinates. 

Using the definition of Go given in Eq. ( |2.90| ), we can integrate the modified 
Bethe-Salpeter equation, in Eq. (|2.88|) , over the light-front energy to get 

d 2 ku ^ dk J 


7(*i/; P) = 


(2tt) 


l -w(ki f , ku ; P)g 0 (ku; P)^(ku] P), 


(2.92) 


where 


w(k lf , ku; P ) = 
l(k i; P) = 


9o(k lf ;P) 


{G 0 WGo)(k lf ,k u ;P) 


9o(ku;Py 


dk { 


■G 0 (fc 1 ;P)r(fc 1 ;P). 


(2.93) 

(2.94) 


go(ki;P)J 2(2tt) 

The functional (/) is defined by its action on an arbitrary function f(k±f, ku), where 
k\f and ku are four-vectors, as 




(2.95) 
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We proceed by calculating the specific form of g 0 , 

Q(ktMK) i 

2 kf 1^2 P~ — ki — 



(2.96) 


where kj = e~{ki). As <70 is a three-dimensional quantity, it is clear that k~ is not 
the independent minus component of a momentum four-vector. With this expression 
for g 0l we can specialize Eq. (|2.92|) to the center-of-momentum frame and obtain 


(E~k lf ~ k 2f ) ip(k lf ; E) 


= / d 2 k 


li,J 


dk 


w(k lf , k u ; E ) 


1 i 


2(2n) 3 s fk^kp4M l 


-^{kii\ E), 


(2.97) 


where 


\E) = ^=£ 7 (^ 1 ; E). 


(2.98) 


By comparing this equation to Eq. (|2.30|) , we find that the full light-front two-nucleon 
effective potential that corresponds to the kernel K , after suppressing the coefficient 
given in Eq. ( |2.57|) . is 


V (kif, k lf ; E) = —w(k lf ,k lf ;E). 


(2.99) 


Thus, we can calculate light-front potentials directly from the Bethe-Salpeter equation 
using this method. 

The potential V can be expanded in powers of the coupling constant, as done 
in LFTOPT. We find that when the auxiliary Green’s function given in Eq. (|2.90| ) 
and the OBE kernel are used, the lowest order parts of the potential (as calculated 
in Ref. [Q) are the same as our OBE and TBE:SB potentials. Thus, we conclude 
that this method produces the physically equivalent Hamiltonian theory to the Bethe- 
Salpeter equation being used. We will use this in the next section to derive a Hamil¬ 
tonian potential for a situation where LFTOPT cannot be applied. 
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2.3.4 The Modified-Green’s-function Approach 


Now that the technology for the three-dimensional reduction has been reviewed, we 
derive an approximate kernel for the Bethe-Salpeter equation. We will follow the 
approach of Phillips and Wallace [£h| and works cited therein. The idea is to start 
with the Bethe-Salpeter equation where the kernel is truncated to only include ladder 


(one-boson-exchange, see Fig. |2]2|) and crossed (two-boson-exchange, see Fig. 
parts, 


F — (Aladder T A" cross )GoF. 


( 2 . 100 ) 


An uncrossed approximation is used where the crossed part of the kernel is approx¬ 
imated by Across ~ A'ladderGc'A'iadder- Our job is to find a valid modified Green’s 
function, Gc- Using this uncrossed approximation, 


r ~ (Aladder T Ai a dderGc , Ai a dder)GoF. 


( 2 . 101 ) 


One can attempt to rewrite Eq. Q2.10C| ) as an equation linear in A' ]ad derj to obtain 
the modified-Green’s-function Bethe-Salpeter equation, 


MGF 


= iu 


ladder 


+ Gc^Fmgf- 


( 2 . 102 ) 


By iterating this integral equation for we obtain 


MGF 


A ladder T ^ ^ A ladder (Gc A ] 


ladder ) 


n= 1 


G n r 


O 1 MGF- 


(2.103) 


The part of Eq. ( 2.102|) that plays the role of the kernel includes the uncrossed approx¬ 
imation of the original kernel A) a dder + AUdderGcAUdder as well as many more terms. 
We note that the higher-order terms approximate some of the higher-order terms that 
should be included in the full kernel, such as three-boson-exchange diagrams where 
several meson lines cross. However, this approach undercounts the higher-order terms 
which it approximates, and also leaves out some terms completely. Therefore, this 
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new Bethe-Salpeter equation will give results that are closer to the full solution than 
Eq. (|2.1001) , but will not give the exact solution. The articles by Wallace and Man- 
delzweig []50|, |5T| demonstrate that this approach, by effectively summing an infinite 
set of interactions, gives the correct one-body limit, which is something that the usual 
Bethe-Salpeter equation with a truncated kernel cannot do. 

The modified Bethe-Salpeter equation in Eq. (|2.102 ) is reduced to a Hamiltonian 
equation via the technique discussed in the previous section. The equal-time Hamil¬ 
tonian has been derived by Phillips and Wallace j43|. They found that this modified- 
Green’s-function approach gave a spectra that lies closer to the full ground-state 
spectra than the other approximations they considered. We will use the light-front 
reduction to obtain the light-front potential for the Hamiltonian equation physically 
equivalent to Eq. ( [2. 102|) . 

To clearly see what role Gc plays, we compare the crossed and uncrossed Feynman 
Using the Feynman rules, 

d k\ m 1 


graphs in Fig. 

-^crossed '-*■ 
-^-uncrossed '-*■ 


(2vr) 4 (, k lf - k lm ) 2 
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2 d(ki rn ^d(k2m) /, 
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hm ) 2 

1 
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(k\i ki m ) /i 


(2.104) 

(2.105) 


where d is the one-particle propagator given in Eq. ( |2.86|) , and k 2 m = ■*fc 2 m + & 2 j + & 2 /- 
The only difference between these two graphs is that the crossed one has d(/c 2 m) while 
the uncrossed one has d(k- 2 m ). 

We want an approximate one-particle Green’s function dc that only depends on 
qi and P, so that 


d c (ki, P) « d(k 2 ). (2.106) 

Substitution of dc(k\]P ) for d(k 2 ) in the crossed graph causes the graph to become 
uncrossed. The penalty for this simplification is that a modified Green’s function 
propagates in the intermediate state, namely id(ki)dc(ki, P). It is important that 
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Figure 2.7: The two-boson crossed and uncrossed Feynman graphs. Here, = 
— Mm + k 2 i + k 2 f ■ 


this approximation be invariant under relabeling particle labels, so we explicitly sym¬ 
metrize by defining 

G c = \ [dikjdciku P) + d{k 2 )d c (k 2 ; P )]. (2.107) 

How should we approximate dc? Since we are interested in obtaining a bound 
state, a low-energy approximation is chosen. Specializing to the center-of-momentum 
frame in this limit, the external momenta are half the total momentum, so Mi = 
k 2 f = P/2 and Mm = P — Mm = Mm- This approximation is similar to the one used 
by Phillips and Wallace. Thus, we define dc{k\\ P) = d(M) so 

G c (ki] P) = Gi(ki, P) + G 2 (M; P), (2.108) 

where we define G\ and G 2 by 

G 1 (k l -P)= l -d(k 1 )\ (2.109) 

G 2 (k l] P)= l -d(P-k 1 ) 2 . (2.110) 

This approximation for Gc is valid for this model for the energy range we study, as 
discussed in Appendix ||. 
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We can write the modified-Green’s-function Bethe-Salpeter equation as 


where 


TmGF — -^ladder Go r M GF) 


Go — Go + G i + G 2- 


( 2 . 111 ) 


( 2 . 112 ) 


This is used as the starting point of three-dimensional reduction discussed in Sec¬ 


tion |2.3.3| , where Go is considered as the Green’s function. Before doing the reduction, 
note that the two poles of G\ and G 2 lie in the same half plane for each function, so 


/ fjU~ ~ 

= g 0 (ki,P). 


(2.113) 

(2.114) 


Proceeding with the three-dimensional reduction of Eq. ([2.111 ) in the center-of- 
momentum frame, we obtain 


(E - k lf - k 2f ) ip(k lf \E) 


= / d 2 k 


1*,-L / u " v l i 

Jo 


dkt. g(fcl f* u;g) jfr u ;E), 

2(27t) 3 y k^jkyk^k^i 


(2.115) 


where 


w{k lf ,k u -P) = . 1 . {G Q WG Q )(k lf ,k li ]P) 1 


go(k lf ;P) 


go(ku]Py 


k 1 -E) = 


dk x 


^kfk+ J 2(2tt) 

and the modihed kernel W is given by 


G 0 (fci; P)r M GF(fci; P), 


(2.116) 

(2.117) 


W = Ki 


ladder 


AT, 


ladder 


Go — Go 


Go(kif,ku]P) — G 0 (kif] P) 


8^ + \k 


1/ 
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W, 


Go(ku; P ). 


(2.118) 

(2.119) 


g 0 (k lf ;P) 

It is a feature of the light front that go = go, so that the uncrossed approximation 


only affects the potential, and Eq. ( 2.115 ) has the same form as Eq. ( 2.92 ). In the 
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equal-time calculation [£T|] g Q ^ g 0 , so the approximation changes both the Green’s 
function as well as the potential. 

We now expand W in powers of the coupling constant, and keep only the lowest 
order term, inlander- According to Eq. ( |2.99|) , the light-front potential that corresponds 
to this truncation of the kernel is the modified-Green’s-function (MGF) potential 
VmGF, 


VmGF^I/, hi] P) = 


-(G 0 Aladde r G , o)(fc 1 /, hi] P)- 


E go(ki f ] P) u “ u/v 'g 0 (ku;P) 

The one-boson-exchange kernel /Gadder is given by the Feynman diagram, so 


( 2 . 120 ) 


A ladder (^1/, 1G E) 

_ (iMf 

(■ hf - hi) 2 - X + irj 

1 

u+ _ u+ 

ft l / ft li 

where the light-front energy of the meson is given by 

u~( q ) = - - h- ( 2 . 122 ) 

By examining the locations of all the poles in the k~ integrals for Vmgf, we find 
the integrals are non-vanishing only when both x and y are between 0 and 1. The 
sign functions in the denominator of /Gadder naturally divide Vmgf into two parts, one 
for x < y and the other for x > y. The integrals in Vmgf are straightforward, but 
quite lengthy and tedious. Therefore, we show only the final answer, 



M 2 


(hf - hi) - Sign (ktf ~ kfju (hf - p u ’ 


( 2 . 121 ) 
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where 


k+ 

Nf, 12 = ~ ^ 1 / — k 2 f), (2.124) 

K V 

k + 

JV/,2i = #( E -k- f -k~ f ), (2.125) 

K l/ 

Di = E - k^- ky - u>~(kif - k 2 i), (2.126) 

D 2 = E - k~ f - k^ - 0J~{ku - k 2f ). (2.127) 


The expressions for W .12 and iVj, 2 i are obtained by replacing / with i in Nf t i 2 and 
N f, 21- 

What is the physical interpretation of this modified-Green’s-function potential? 
The first term multiplying each 6 function gives the OBE potential we derived before 
from the perturbation theory. There the D in the denominators corresponds to one 
meson exchange. The second and third terms multiplying the 6 functions, with D 2 and 
D 3 in the denominators appear to be effective two- and three-meson-exchange terms. 
Since time-ordered perturbation theory does not apply to the modified Bethe-Salpeter 
equation that we use, the exact nature of these terms is not easy to understand. 
However, it is clear that these terms increase the strength of the potential, and should 
mimic the higher-order diagrams that are not being included explicitly. 

The only dependence on the direction of the perpendicular components of kf and 
ki comes from the D 1 s. This allows the azimuthal-angle integration of Vmgf to be 
done easily, as shown in Appendix |C.l . 


2.4 Results for the Wick-Cutkosky Model 

For our numerical work, we pick the meson mass to be 0.15 times that of the nucleon, 
so ji = 0.15M. This is chosen so that our ground state can be considered a toy model 
of deuterium, and also to facilitate comparison with the results of Nieuwenhuis and 


Tjon [221, Phillips and Afnan [201, and Schoonderwoerd, Bakker, and Karmanov [14 . 
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Nieuwenhuis and Tjon used the Feynman-Schwinger representation (FSR) of the two- 
particle Green’s function 0] in the quenched approximation without the mass and 


vertex renormalization terms [221. Their result is to be considered the full solution 


that the Bethe-Salpeter and Hamiltonian equations approximate. For the Bethe- 
Salpeter equation, computation of the bound-state energies for models similar to ours 


have been done for the ladder |53| and ladder plus crossed |35| kernels over 30 years 
ago. More recent results are found in Refs. [^(||54|, where the solutions are compared 
to those given by the FSR approach. 


2-4-1 Ground State Results 

We want to examine the ground state energies obtained for each of the light-front 
potentials derived in this chapter, and see how they compare the energies calculated 
in other approaches. Our main interest is in finding out how well the light-front 
potentials approximate the Bethe-Salpeter results, and also in how well both the 
perturbative and non-perturbative potential agree with the Feynman-Schwinger rep¬ 
resentation approach. 

Consider first how the light-front Hamiltonian approach fits in with the other ap¬ 
proaches. As discussed in section [2.2.3| , different light-front potentials can be derived 
from Bethe-Salpeter equations with different kernels. We have mentioned that the 
OBE+TBE:SB potential should approximate the ladder Bethe-Salpeter equation, and 
similarly the OBE+TBE potential should approximate the Bethe-Salpeter equation 
when the ladder plus crossed kernel is used. The best that these truncated Hamilto¬ 
nians can do is approximate their respective Bethe-Salpeter equations. 

With this in mind, we evaluate the coupling constant versus bound-state energy 
curves (which we will call the spectrum) for the Hamiltonian equation with the OBE 
potential, the OBE+TBE:SB potential, and the OBE+TBE potential. For the range 
of values we use here, we find numerical errors in the value of g 2 are less than 2%. 
Our results (without error bars) are plotted along with the results obtained with 
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the ladder BSE [[2I|, and ladder plus crossed BSE |^TJ in Fig. We note that 
the OBE+TBE:SB potential agrees well with the ladder BSE, and the OBE+TBE 
potential agrees with the ladder plus crossed BSE. This is the best that a Hamiltonian 
can do, so this result is interpreted as evidence that, in general, the higher-order 
diagrams are very small for the ground state on the light front. 


If all one wanted was a way to approximate the spectra for Bethe-Salpeter equa¬ 
tions with different kernels, one could just use the truncated potentials that the 
LFTOPT provide. However, the true goal is to approximate the spectra for the full 
ground state, which in this model is given by the FSR approach |22|. We expect that 
the non-perturbative potentials should give a better approximation of the full solu¬ 
tion than the perturbative potentials, since the non-perturbative potentials attempt 
to incorporate physics from higher-order diagrams, although this is not immediately 
clear by looking at the forms of the potentials used. We plot the results for all of 
the light-front potentials described in this chapter, the three truncated potentials 
(OBE, OBE+TBE:SB, and OBE+TBE) and the four non-perturbative potentials 
(symmetrized mass, retarded, instantaneous, and modified Green’s function) along 
with the results for the full theory in Fig. [2.9| . For deeply-bound states, there is con¬ 
siderable disagreement between the perturbative results and the full results, while the 
non-perturbative results are in better agreement, with the modihed-Green’s-function 
(MGF) potential achieving the closest agreement. For lightly-bound states, the results 
for all of the potentials appear converge to each other, close to the full result. 


For the modihed-Green’s-function potential, only the Erst term of the expansion 
in g 2 was kept. In principle, higher-order terms could be calculated. However, since 
there was fairly good agreement between the OBE potential and the ladder Bethe- 
Salpeter equation, it is expected that the MGF potential will give results that are 
close to the ladder Bethe-Salpeter equation using the modihed-Green’s-function. 
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E/M 

Figure 2.8: The spectra calculated using the OBE, OBE+TBE:SB, and OBE+TBE 
potentials are displayed in this figure, along with the spectra for the Bethe-Salpeter 
equation for the ladder and ladder pins crossed kernels. The curves for the lad¬ 
der (ladder plus crossed) Bethe-Salpeter equation and the light-front OBE+TBE: SB 
(OBE+TBE) potentials are very close to each other, almost indistinguishable in this 
figure. E is the energy of the ground state of two nucleons, and M is the mass of 
the nucleons. The meson mass is n — 0.15 M. The binding energy B is related to the 
bound state energy via E = 2 M — B. 
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Figure 2.9: The spectra for the seven light-front potentials (OBE, OBE+TBE:SB, 
OBE+TBE, retarded, instantaneous, symmetrized mass, modified-Green’s-function) 
are plotted here along with the spectra for the full solution (FSR). The instantaneous 
curve lies very close to the symmetrized mass curve, hence the same line style is used 
for both. 
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2-4-2 Excited State Results 


Now, we are ready to address the problem of rotational invariance of the light-front 
potentials, and see how the bound states are affected by the breaking of rotational 
invariance. We expect a splitting of the bound-state energies for different m states, 
and we also expect that the states will not be pure angular momentum states. In 
order to examine these effects, we must have excited states, preferably those with 
non-zero angular momentum. 

Furthermore, since the results of the previous section show that the ground-state 
energies calculated with the OBE+TBE potentials are in excellent agreement with 
the Bethe-Salpeter results, we expect that the excited states calculated with the light- 
front potentials will at least be similar to the Bethe-Salpeter results. This would mean 
that the excited states should almost be angular momentum eigenstates and could be 
classified easily. 


First, we consider the solution for the excited states of the Bethe-Salpeter equa¬ 
tion. The technology for doing these bound-state Bethe-Salpeter equation calculations 
was developed over 30 years ago [|33|, Q. Since the eigenstates of the ladder Bethe- 
Salpeter equation are also eigenstates of the total angular momentum, there is exact 
degeneracy in the energies of the different m states for the same angular momentum. 
For the range of parameters used in this study, we find that the numerical errors in 
g 2 are less than 0.5%. The numerical errors are largest for the most deeply-bound 
states (E ps 1.85M) with the largest coupling constants (g 2 /An pa 50). 

The solutions to the ladder Bethe-Salpeter equation form bands where excited 
states with different orbital angular momentum are approximately degenerate with 
each other, as shown in Fig. p.lOj . This approximate degeneracy is due to the fact 
that when /i —> 0, this model can be shown to have the same degeneracies as the 
non-relativistic hydrogen atom [[ 


Thus, in that limit, all states with the same 


principal quantum number have the same energy. When /j ^ 0, that degeneracy is 
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broken, but only slightly, as we can see from Fig. m Because of this, we will label 
our states using atomic spectroscopy notation. 

Next, we consider the two light-front Schrodinger equations given by Eqs. 


and (|2.48D . These equations are solved numerically (for each parity and several m 
values) for the spectrum g(E) for a range of bound-state energies E. The symmetries 
of the light-front Hamiltonian allow us to classify the states according to m and the 
action under parity, so as an example, we calculate the spectra with even parity and 
rn — 0. For the range of values we use, we find the numerical errors in g 2 are less 
than 2%. The errors are largest for the most deeply bound states with the largest 
coupling constants, as was the case for the Bethe-Salpeter equation. 

We plot both the OBE spectra [in Fig. [2.111 (a)] and the OBE+TBE:SB spectra 
[in Fig. [2.111 (b)] along with the Bethe-Salpeter spectra for the even parity and rn — 0 
states. Based on the energies, we see that the n — 1 state is expected to be the Is 
state, the n — 2 state the 2s state, and the n = 3 and n = 4 states the 3s or 3d 
states respectively. We see approximate agreement between both of the truncated 
Hamiltonian results and the Bethe-Salpeter results for lightly-bound systems (where 
E « 2 M). 

Our states are not manifest eigenstates of total angular momentum J 2 . However, 
if our approximation of the Hamiltonian potential was good enough, we would be able 
to identify the n — 3 and n — 4 states of the light-front Hamiltonian calculation with 
3s and 3d states of the Bethe-Salpeter equation unambiguously. This would determine 
the angular momentum of the light-front Schrodinger equation wave functions. From 
Fig. |2TT], we see both that the addition of the TBE:SB potential brings the n = 3 
and n — 4 states closer to the Bethe-Salpeter results, and that in the lightly-bound 
region the n = 3 and n — 4 states can be identified as 3s and 3d states. For more 
deeply-bound states, this identification cannot be made. 

An alternative approach to assigning angular momenta labels would be to perform 
a partial-wave decomposition of the wave functions in the real angular momentum ba- 
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Figure 2.11: Spectra for the four lowest wave functions of even parity and m — 0 (a) for 
the OBE equation (short-dashed line) and the Bethe-Salpeter equation (solid line), (b) 
for the OBE+TBE:SB equation (long-dashed line) and the Bethe-Salpeter equation 
(solid line). The curves are labeled by n, which indicates that the curve belongs to the 
n th eigenvector. The curve for the Is state for the Bethe-Salpeter equation and the 
curve for the first OBE+TBE:SB wave function are very close together, but distinct. 
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sis. The wave functions conlcl then be labeled by the angular momentum component 
which is dominant. In light-front dynamics, this is difficult to do since the perpendic¬ 
ular components of the real angular momentum operator (.J x and J y ) are dynamical, 
which makes the angular momentum operator as complicated as the Hamiltonian. We 
are encouraged to look for an alternative operator which is easier to use, yet approxi¬ 
mates the behavior of the real angular momentum. We shall use the “relative angular 
momentum”, defined by L = x x k, where the vectors are equal-time three-vectors. 
The relative momentum k is defined by Eq. (|2.31|) and x is canonically conjugate to 
k. The “relative angular momentum” L is used to help analyze our solutions, and is 
not the same as the real angular momentum that can be derived from the Lagrangian 
using the energy-momentum tensor. However, we will see that the “relative angu¬ 
lar momentum” gives results that are expected from the real angular momentum, so 
our use of the “relative angular momentum” in place of the real one appears to be 
justified. 

A partial-wave decomposition is performed on the wave functions (represented in 
the relative momentum basis) to obtain the radial wave functions Ri >m for all “relative 
angular momentum” states Yj m . Since the potential is not manifestly rotationally 
invariant when written in terms of the relative momentum, our wave functions will 
have support from many different partial waves. Not all partial waves are allowed; 
only those partial-wave states with the same J 3 and parity quantum numbers as the 
wave function give non-vanishing radial wave functions. We have 


OO 


(*#”,„) = ]Ti T(«,<P)RZ(k)- 


(2.128) 


We define the fraction of the wave function with “relative angular momentum” l as 



(2.129) 


The fractions /” are a measure of the amount of “relative angular momentum” state 
l in the n th eigenfunction. 
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Figure 2.12: The fractions of each angular momentum state in several wave functions 
is plotted as a function of the binding energy E/M. The wave functions chosen are 
the four lowest coupling constant wave functions with m = 0 and even parity. The 
results for the OBE truncation are shown on the left, (a-d), and the OBE+TBE:SB 
truncation are shown on the right, (e-h). Note that the OBE wave functions, in 
general, have more support from a larger number of partial waves. 
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To illustrate these “relative angular momentum” fractions, we perform this anal¬ 
ysis on the four lowest coupling constant wave functions with even parity and m — 0 


for both the OBE and OBE+TBE:SB truncations, the same as in Fig. |2.11| . We 
show /” as a function of E in Fig. |2.12| . (The plots for the m = ±1, ±2 states show 
similar behavior, with different l values.) Note that more higher “relative angular mo¬ 
mentum” states contribute for deeper bound states. This is the same behavior that 
would be expected if the real angular momentum was used to perform the partial 
wave decomposition. Recall that the real angular momentum will commute with the 
full potential if no truncation is made. However, the potentials we use are truncated, 
and neglecting the higher-order terms breaks the rotational invariance of the poten¬ 
tial. Also, both the binding and the importance of the higher-order diagrams increase 
with the coupling constant. Thus, both truncations we consider will break rotational 
invariance more when the states are more deeply bound. In the absence of having 
the real angular momentum fractions, we will use the “relative angular momentum” 
fractions. 

Examination of the /" curves in Fig. |2.12| shows us that, as postulated earlier, the 
n = 1 and n — 2 states shown in Figs. 2.12 (a,e) and |2.12| (b,f) are predominately s- 
wave states, which we label the Is and 2s states respectively. For the OBE truncation, 
the n = 3 state shown in Fig. |2.12| (c) is predominately d-wave (labeled the 3d state) 
and the n — 4 state shown in Fig. j2.12| (d) is predominately s-wave (labeled the 3s 
state), with little mixing between the two states. For the OBE+TBE:SB truncation, 
the n = 3 and n = 4 states are mixtures of both the 3s and 3d states, as seen in 
Figs. |2.12| (g,h). In fact, we see a level crossing between the n — 3,4 states in that 
the n — 3 state is predominantly 3d for lightly-bound systems, but as the binding 
increases, the n = 3 state becomes predominantly 3s. Again, this type of mixing 
would also be found if the real angular momentum was used. 

If this 3s-3d mixing is ignored, then Fig. |2.12| shows clearly that when the TBE:SB 
is included the amount of “relative angular momentum” states mixed in actually 
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decreases as compared to the OBE results. This implies that the wave functions 
of the OBE+TBE:SB equation are better “relative angular momentum” states than 
the wave functions of the OBE equation. Thus the TBE:SB potential restores some 
rotational invariance to the OBE calculation. 

We attempt to classify the wave functions as states with definite angular momen¬ 
tum. If a wave function has the most support from the “relative angular momentum” 
state l, we classify that state as having angular momentum /. This procedure desig¬ 
nates the first two wave functions for both the OBE and the OBE+TBE:SB cases as 
s-wave states, which is clearly the right thing to do. For the third and fourth wave 
functions of the OBE+TBE:SB equation, there are points where the fraction of the 
s-wave state equals the fraction of the d -wave state. Here it no longer clear that such 
a state should be assigned a dehnite “relative angular momentum” value; we do so 
regardless and analyze the consequences later. 

We now examine the breaking of rotational invariance of the two truncations 
based on the comparison of states with different m projections of the same angular 
momentum. Since the s-wave state has only m = 0, there is nothing to analyse in that 
case. Further discussion of the ground-state s-wave appears in Ref. [[57J. In Figs. |2. 1 


2.14|, and |2.15| , we plot the Bethe-Salpeter bound-state spectra along with the spectra 
for the states constructed with the OBE and the OBE+TBE:SB potentials. The 
different m states for the Bethe-Salpeter equation are exactly degenerate as a result 
of the rotational invariance of the equation. The curves from the Hamiltonian theory 
do not exhibit the exact degeneracy in m; the degeneracy is broken whether the 
OBE or OBE+TBE:SB potentials are used. However, we see that the degeneracy is 
partially restored when the TBE:SB potential is included, in that for a given binding 
energy the spread of the coupling constants is always smaller. 

In Figs. |2.13| and |2.14| , we plot the spectra for the 2 p and 3 p states, respectively. 
In both figures, we see that the m — 0 and m — ±1 curves move closer together after 
addition of the TBE:SB potential. However, the average of the two curves for each case 
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does not move much. In fact, for the 3 p state in Fig. [2.14 , the m — ±1 curve moves 
farther away from the ladder Bethe-Salpeter curve after addition of the TBE:SB. We 
also note that in Fig. p,14| the spread of the spectra is larger than in Fig. |2.13| for both 
the OBE and the OBE+TBE:SB light-front Schrodinger equations. We attribute this 
to the neglect of the higher-order graphs. Since the states with different m values 
would be degenerate if all the higher-order graphs were included, not including them 
causes a breaking of the degeneracy. Because the coupling constant is larger for the 
3 p states than for the 2 p states with the same binding energy, the omission of the 
higher-order graphs causes a larger breaking of the degeneracy for the 3 p states. 

We also consider the first d-wave states. We see in Fig. |2.15| that, as in Figs. |2.15 


and |2.14| , the states with different m values move together after addition of TBE:SB. 
The effect of the level crossing of the 3s and 3d states shown in Figs. |2.12| (g) and 
j2.12| (h) is seen here in that there are two disjoint curves for m = 0. When the level 
crossing occurs, at E/M rs 1.96, the wave function which originally was the 3s state 
becomes the 3d state and vice versa. 


There are some problems that this level crossing causes. In the intermediate region 
where the 3s and 3d states are about equally mixed, it is probably not physically 
sensible to call the state a 3s or 3d state. However, plotting the curves gives us 
an indication of what the wave functions are doing, and for this case we see that 
the m = 0 curve is always bounded by the m = ±1 and m = ±2 curves. A more 
restrictive classification scheme would leave a gap (in E/M) between the two m = 0 
curves and it would not be so clear that the bounding of the m = 0 curve which we 
observe here occurs. 

Finally, we note that the spread of the curves in the 3d case is approximately 
the same as the spread in the 3 p case for each truncation. This tells us that the 
higher-order graphs have the same qualitative effect in both states. 
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E/M 

Figure 2.15: The first d -wave state, the 3d state. The Bethe-Salpeter (BS) result is 
plotted with a solid line. In (a) the OBE results are the short-dashed lines, and in (b) 
the OBE+TBE:SB results are the long-dashed lines. Note that there are two separate 
m — 0 curves when the OBE+TBE:SB are used, one extending a quarter of the way 
over, and the other hiding under the m = ±1 curve with the OBE+TBE:SB. 
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2.5 Conclusions 


We have examined several aspects of bound states in a massive Wick-Cutkosky model 
using several light-front potentials. First, we used three different truncations of the 
effective potential derived from the perturbative held theory and four approximations 
that attempt to incorporate non-perturbative physics. For each of these potentials, 
we calculate the coupling constant that gives the ground-state for a given bound- 
state energy (the spectrum), and compare to the spectra from different approaches 
found in the literature. We find fairly good agreement between all the methods for 
lightly-bound systems. 


For the full range of binding energies studied, the results for calculations including 
one- and two-boson-exchange potentials agree with the Bethe-Salpeter equation using 
the physically equivalent truncation of the kernel within the numerical errors. (This is 
a consequence of examining the ground state. For the excited states more higher-order 
light-front time-ordered graphs are required to get the same level of agreement ||36|| .) 
The agreement for the case with the stretched box diagrams (OBE+TBE:SB) has 
been shown previously by Sales et al. [[8|; the result with the crossed-box contribution 
(OBE+TBE) is new. This excellent agreement with the Bethe-Salpeter results has 
an undesirable consequence: The BSE results are known to be a poor approximation 
of the full solution for deeply bound systems, so the truncated Hamiltonian approach 
cannot provide a good approximation to the full solution in that regime. 


The non-perturbative potentials based on physical considerations give a better ap¬ 
proximation of the full solution than the potentials obtained from LFTOPT. For all 
binding energies, the modified-Green’s-function potential achieves the closest agree¬ 
ment with the full solution of all the potentials considered here. However, there is 
still considerable disagreement between approximate potentials and the full result for 
deeply bound states. We interpret this as an indication that the approximations, 
while incorporating some non-perturbative physics, do not go far enough. In the 
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weakly-bound regime, which is of relevance for deuteron calculations, the spectra for 
all of the potentials are close together, indicating that light-front dynamics provides 
a good description of lightly-bound systems. 

After examining the ground state thoroughly, we considered the low-lying excited 
states of the Wick-Cutkosky model. We considered two truncations of the light-front 
Hamiltonian, the OBE and OBE+TBE:SB truncations. Using these truncations, a 
“relative angular momentum” operator was used to study the partial-wave decom¬ 
positions of the bound-state wave functions. We found that the “relative angular 
momentum” operator acting on those states yield behavior similar to that expected 
from the real angular momentum. This result encourages the use of this operator to 
classify the states according to their angular momentum values /, and to study the de¬ 
generacy of the spectra for states with the same l value but with different projections 
m. We found less breaking of the degeneracy when the OBE+TBE:SB potential was 
used than when the OBE potential was used. Both of these findings indicate that the 
OBE+TBE:SB truncation of the Hamiltonian breaks rotational invariance less than 
the OBE truncation alone. 

However, there is still some discrepancy between our truncated Hamiltonian spec¬ 
tra and the ladder Bethe-Salpeter spectra. In general, our OBE+TBE:SB results for 
the p- and d -wave states show deeper binding than the Bethe-Salpeter results. Not 
surprisingly, this disagreement is largest for the most deeply bound-states where the 
coupling is largest. This difference would be removed if all the higher-order pieces 
of the potential were included. These higher-order pieces are also needed for the full 
restoration of rotational invariance. 
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Chapter 3 

REALISTIC NUCLEAR MODELS 


In this chapter, we start with an effective nuclear Lagrangian [l8|,|)5|] for nucleons 
and six mesons which incorporates a chiral symmetry. We derive the light-front 
Hamiltonian from the Lagrangian and address some of the difficulties encountered. 
Then the Hamiltonian is used to construct the Feynman rules for nucleon-nucleon 
potentials. The issues of the breaking of rotational invariance and chiral symmetry 
by certain subsets of the potentials are considered and addressed. We use the Feynman 
rules and the symmetry considerations to construct new one-meson-exchange (OME) 
and two-meson-exchange (TME) light-front nucleon-nucleon (LFNN) potentials. 


The potentials are used in numeric calculations for two families of nuclear models. 
For our first set of calculations, we assume that the pion is the only meson that 
interacts with the nucleons, giving us a special case of the full nuclear model. This 
pion-only model is inspired by the non-relativistic one-pion-exchange model used by 
Friar, Gibson, and Payne [Q. We generalize their model to form the basis of our 
pion-only light-front model, which includes relativity automatically. 


Our second set of calculations uses the full model with all six mesons ( 7 r, cr, p, 
u, p, and h) for the interaction. This model is an extension of the light-front model 


used by Miller and Machlcidt [55]. A new feature of this model is that light-front 


energy dependent denominators are used in the potentials; the denominators used in 
Ref. |5§] are energy independent. 
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3.1 Model and Formalism 


Our starting point is a nuclear Lagrangian |T8[ which incorporates a non-linear chiral 
model for the pions. The Lagrangian is based on the linear representations of chiral 
symmetry used by Gursey |57]]. It is invariant (in the limit where m n —*■ 0) under 
chiral transformations. 

The model prescribes the use of nucleons -0' (or 0, the free nucleon field, as 
discussed in Section |3.3.2|) and six mesons: the ir, 8 (also known as the do(980)), o 
(also known as the /o(400 — 1200)), p, p, and u mesons. The coupling of each meson 
to the nucleon is governed by the combination of the meson’s spin and isospin. The 7r 
and r) are pseudoscalars, the p and u are vectors, and the S and a are scalars. Under 
isospin transformations, the 7r, p, and 5 are isovector particles while the p, u, and a 
are isoscalar particles. 

The use of scalar mesons is meant as a simple representation of part of the two- 
pion-exchange potential which causes much of the medium range attraction between 
nucleons 0,^]. It can also be interpreted as the effect of fundamental scalar mesons 

1M1- 


The Lagrangian C is based on the one used in Refs. 18,55, 53|. It is given by 


fp'"' ' 


m 


V “ • P M - + 


m; 


> P 1 


+ i/ 2 Tr(^f/<W f ) + ^m 2 / 2 Tr (U + I/ f - 2) 

+ ^(d p ad^a - mla 2 ) + \{d p 8 • d^S - m]8 2 ) + ^(d^d^p - m 2 p 2 ) 


+ 0 


2 ° ' 2 
- g p p ■ r - g w u p ) — U(M + g a a + g s 8 ■ r + ig^ 5 p) -0', (3.1) 


where the bare masses of the nucleon and the mesons are given by M and m a where 
a = 7r,p,cr, <5, p, u. We have defined V^ u = d ll V v — d u V tl for V = p, u. The unitary 
matrix U can be chosen to have one of the three forms Up 

^ 1+U5T-W2/ , 

1 - 77 5 T • 7T/2J 


U 3 = y/l - 7T 2 // 2 + iy 5 T • 717/, (3.2) 
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which correspond to different definitions of the fields. Note that each of these defini¬ 
tions can be expanded to give 


U 


1 + ^75 


r • 7r 

~T~ 



(3.3) 


In this work, we consider at most two meson exchange potentials, so we consider U 
to be defined by Eq. (|3.3p. 

In the limit where m n —* 0, this Lagrangian, is invariant under the chiral trans¬ 
formation 


ip' -> e* 75 ’ r ' a, 0 / , 


U 


,—vysT-a 


Ue 


—t'YsT-a 


(3.4) 


In this model the other mesons are not affected by the transformation because they 
are not chiral partners of the pion. This is in contrast to the Lagrangian given in 
Refs. ||T8l , |63|| , where the mass and scalar interaction terms for the nucleon were written 
as MU + g s (f) instead of U(M + g s (f>). 


3.2 Non-interacting Nucleon-Nucleon Theory 


The light-front Hamiltonian is derived from this Lagrangian using the same approach 
used in section O- the approach used by Miller |[TF| and others pd| |H7j| . The basic 
idea is to write the light-front Hamiltonian ( P ~) as the sum of a free, non-interacting 
part and a term containing the interactions. We consider the free part first. 


3.2.1 Free Field Expansions 

The solutions for the free fields are similar to those obtained by using equal-time 
dynamics. In fact, the solutions are formally related by a change of variable, and so the 
most obvious difference between the two is due to the Jacobian. The held equations 
have the general form (when Lorentz, spinor, and isospin indices are suppressed) of 


(3.5) 
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where a = tt, rj, a, 5, p, u, ip. Note that in the exponentials, 


kJ'Xn = - (k + x + k x + ) — k_i ■ x±. 

In particular, the solutions for all the mesons and the nucleon held are 
d 2 k ± dk+6(k H 


TT(X) = 

rj(x) = 
6(x) = 
a(x) = 


/^(x) = 


IV^ix) = 


(2vr) 3 / 2 v / 2 k+ 
d 2 k ± dk+6(k + ) 
(2vr) 3 / 2 v / 2 k+ 
d 2 k ± dk+6(k + ) 

(27t) 3 / 2 v / 2 k+ 

d 2 k±dk + 9(k + ) 


[a n (k)e- ik ^ + a\{k)e +ik "*'*] , 
[a v (k)e~ ik ^ + a\(k)e +ik ^] , 
a s {k)e - ik ^ + a\{k)e +ik ^ 
[a a (k)e~ ik ^ + al(k)e +ik ^] , 


(27r) 3 / 2 v / 2fc q 

d 2 k.dk u^ y ^ s)e +^H , 

(2vr) 3 / 2 v / 2fc+ fl t^ 3 L 

rf 2 fc ± dfc+^+) y + at(fc, S )e + ^l , 

(2vr) 3 / 2 v / 2F s t^ 3 L 


Hx) = V2M f 

J ( 2~i :i - v'2/, ,+ 


(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 


x Xi [“(*!. A)6(fe)e-“" z ' + «(fc,A)rf , (*)e+“" I ‘']xi,. (3.13) 

A=+,— £ 3 =+,— 

The polarization vectors are the usual ones. The most general of the commutation 
relations is 


ti {k,s),alAk',s') = <5 0li at5,jt5 s y<5 t2 ' +) (A: - fc')> 


(3.14) 


where a, i, and s denote the meson type, isospin, and spin, and 5^ 2,+ ^ is defined by 
Eq. ( |2.9|) . The anti-commutation relations are 


{6(Jfe, A),6 t (fc',A')} = {d(k, A),d t (fc , ,A / )} =6 X}X ,6 {2 ' + \k-k'). 


(3.15) 


All other (anti-)commutation relations vanish. The spinors are normalized so that 
u(p, X')u(p, A) = 5 X i X . For more information on the definition of the spinors, see 
Appendix [E. 
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3.2.2 Non-interacting Hamiltonians 

The general form of the non-interacting Hamiltonian for each meson is 

k 2 , + m 2 


P 0 (a) = / d 2 k^dk + 9(k + )a) a (k)a a {k)- 


k + 


(3.16) 


For the vector mesons (p and u>), there is an implicit sum over the meson spins. 
Explicitly, this means that for vector mesons a' v (k)av(k) —*■ ^ s=13 a' v (k, s)av{k, s). 
Likewise, the sum over the isospin of the isovector mesons ( 7 r, S, and p) is implicit. The 
sum over isospin can be made explicit by writing a\(k)aj(k) —> JT =13 a\ i (fe)aj ! i(fe). 
The non-interacting Hamiltonian for the nucleons has a similar form as well, 

k =5 + M 2 


p-(^) = / d 2 k ± dk + 9(k + ) 


J2 b\k,X)b(k,X) + d\k,X)d(k,X) 

LA=+,- 


k+ 


(3.17) 


These equations are what one expects, since a free particle with momenta k± and k + 

-\~TTL' ^ 

has light-front energy k~ = —. 


3.3 Interacting Nucleon-Nucleon Theory 

The interaction Hamiltonians are derived from the Lagrangian in Eq. (|3.1| ) using the 
techniques presented in chapter [^. However, there are some additional complications 
due to the structure of the interactions. 

One complication is that the chiral coupling of the pion field to the nucleons 
through the U matrix generates vertices with any number of pions. This is addressed 
simply by expanding the U matrix in powers of j, and considering the interaction 
Hamiltonians order by order. 

Another complication is due to the fact that both the vector mesons and the 


fermions have components which depend on other components of the field l|lq. pf- - 
Tip. Vector meson fields have four components, but only three degrees of freedom. 
Likewise, fermion fields have four spinor components, but only two degrees of freedom. 
When the dependent components are expressed explicitly in terms of the independent 
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components, we obtain new (effective) interaction Hamiltonians for instantaneous 
vector mesons and fermions. A complete derivation is given by Miller in Ref. JR|, we 
illustrate only the main points of the derivation here. 


3.3.1 Expanding the Pion Interaction 


We start by Taylor-expanding U in powers of j, after which the derivation of the one- 
meson-interaction Hamiltonian P'J 1 is straightforward. (The prime indicates that it 
is in terms of i/j', not if). We derive the expressions for Pf x in the section |3.3.2| .) The 
result is 


P'i,i= I d 2 x±dx i/j(x) g P ^Pn,i(x)Ti +g^uj^x) + g s Si(x)Ti +g a a(x) 

+ 9n(i'Ys)Ti7ri(x) +g x (i^ 5 )x(x) ip'(x). (3.18) 


We have defined a dimensionless coupling constant g n = y. To save space and to 
generalize, we define 


Rv = 


T = 

-I rv 


iy 5 if a is a pseudoscalar meson (n, rj) 
1 if a is a scalar meson (5, a) 

y M if a is a vector meson (p, u>) 

Ti if a is an isovector meson ( 7 r, 5, p ) 

1 if a is an isoscalar meson ( 77 , cr, a;) 


(3.19) 


(3.20) 


and denote the meson fields by <h a . This allows us to write 


Pja = y, j d 2 x ± dx ^\x)g a T a T a ^ a {x)^\x), 

ot=7r,ri,<j,S,p,uj 


(3.21) 


where the appropriate sums over the meson indices are implicit. 


The next step is to derive the two-meson-interaction Hamiltonian which arises 
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from chiral symmetry, P'j 2c : 


P’i 2c — / d 2 x±dx if i x ) 


9n 


+ 


TiT j 'Ki(x)'K j {X) 

(i^ b )TiTTi(x)<i){x) 


2 M 

9tt9(I> (2^,5 

~M 


i(x)x{.x) 'fp'(x) 


-s a I d 2 x±dx ip , (x)[T^T^ 7T (x)][T a T a ^ a (x)]if , (x), 


(y=7T,T],(J,S 


M 


(3.22) 

(3.23) 


where s a is a symmetry factor, equal to \ when a = 7r, and 1 otherwise. When the 
contact interaction is used to calculate diagrams, an additional factor is picked up for 
the TiTi contact term (due to indistinguishability) which cancels the symmetry factor 
s., r . Note that this contact interaction involves only scalar and pseudoscalar mesons. 


3.3.2 Elimination of Dependent Fermion Components 

“It does not do to leave a dragon out of your calculations, if you live 
near him.” 

— J.R.R. Tolkien 


We are now ready to express the dependent components of in terms of the 
independent components, and address the problem of instantaneous nucleons. The 
generation of instantaneous interactions is a general feature of theories with inter¬ 
acting fermions in light-front dynamics. This allows us to use a simplified model 
to demonstrate how these instantaneous nucleons arise. In particular, we want to 
postpone the discussion of the complication introduced by the vector mesons until 
the next section. To this end, we choose to remove all mesons except the a from the 
Lagrangian given in Eq. (|3.1| ). (The cr is chosen since it has the simplest coupling to 
the nucleon.) The equation of motion for the nucleons is then 


ipf>' = (M + gao)^'. 


( 3 . 24 ) 
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Applying the projection operators A± = l7°7 ± (defined in Appendix |A[) to Eq. (|3.24]) 
splits it into two equations, 


id ip' + = [«_l • p ± + (3(M + g a a)}^'_, (3.25) 

id + ^'_ = [a.± ■ p ± + /3(M + g a a)}^' + , (3.26) 

where 0 ± = A ± 0 '. This split is useful because in Eq. (|3.24|) , all four components of 
the nucleon field are interrelated, while in Eqs. (|3.25|) and (|3.26|) , the two components 
of 01 are related to the two components of 01, and vice versa. 


First, notice that Eq. (|3.25|) involves d ~, a dynamic operator in light-front dy¬ 
namics. Dynamic operators should be avoided since they involve the interaction, and 
are therefore complicated. We use Eq. (|3.25|) to avoid that complication and relate 
the components of 0. Secondly, to keep the relation as simple as possible, we do not 
attempt to invert the spinor matrix on the right side of Eq. (|3.25 ). Requiring that 
the equation for the dependent components be both a kinematic equation and sim¬ 
ple equation forces us to choose 0(_ as the independent components. The dependent 
components, 01, are defined by 


0 - = [ax • P± + (3(M + g a cr)} 0 +. 


(3.27) 


Notice that the dependent components consist of a non-interacting part and an inter¬ 
acting part. We separate these parts by defining 0 (without a prime) to be the free 
nucleon field, and £ to be the part of 0' that is due to interactions. So 


where 


pi 


0' = 0 + f, 

(3.28) 

7+ 

Pl9aO-W+ = 2 + (^<t)0+. 

(3.29) 

7+ 

: ^ + 2 p+ (w^)0+ 

(3.30) 

: ^ + 2 p + ^ 

(3.31) 


This allows us to write 
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The last equation follows is obtained from noting that ( 7+) 2 = 0, which implies that 
= 7 +, 0 +. 


Plugging Eq. (|3.31|) into Eq. (|3.18|) and removing all mesons except the cr meson, 
we obtain 


P i,i ~ Pi, i + Pi,2i 
Pf A = [ d 2 x ± dx-^(x)g a a(xMx), 


Pi 2 = / d 2 x±dx ip(x) 


g a (7{x) 


7 


2 p+ 


9acr(x) 


lf(x) 


(3.32) 

(3.33) 

(3.34) 


We interpret the factor as a type of nucleon propagator that joins any two 
meson interactions (although having two of these propagators adjacent to each other 
causes the interaction to vanish since (7 + ) 2 ). Because this propagator does not allow 
for an energy denominator (as it is already between two potentials), is called an 
instantaneous propagator. 

Thus, when constructing the diagrams for the light-front potentials, we must also 
include instantaneous propagators for the nucleons in addition to the usual propaga¬ 
tors. (This is one of the dragons mentioned earlier.) 


3.3.3 Elimination of Dependent Vector Meson Components 


Like the nucleons, the vector mesons have a dependent component that contains 
interactions and must be eliminated. This process is complicated somewhat by the 
fact that the dependent nucleon components must be eliminated at the same time. 
The salient points of the combined elimination of the dependent nucleon and vector 
meson components are discussed in detail by Miller Jl 8 | . 


The result is that the vector meson field must be redefined and an instantaneous 
vector meson propagator is generated in addition to an instantaneous nucleon prop¬ 
agator. However, when the nucleon-nucleon potential is calculated, the redefinition 
of the vector meson field exactly cancels the contribution of the instantaneous vector 
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meson. The result is that the potentials can formally be calculated using the original 
vector meson field. 

In this work, we use that result to simplify our derivations of nucleon-nucleon 
potentials by formally using the native form of the vector meson field. Thus, we find 
an interaction Hamiltonian similar to the one shown in Eq. ( |3.34| ), 

P T1 = P r] + Pj 2 , (3.35) 

Pf y 1 = J d 2 x±dx~'ip(x)g a r a T a $ a (x)'ip(x), (3.36) 

ot=7r,r),cr,5,p,uj 

PT;. 2 = J d 2 x_ L dx~'i/j(x) [g ai T ai T ai ^ ai (x)] ^ 

ol\ i&2='K,r),cT,5,p,uj 

[g a2 T a2 T a2 ^ a2 (x)]p(x). (3.37) 


We may continue to interpret as an instantaneous nucleon propagator, since in 
the derivation in Ref. Jl8| it is clear that the potential vanishes when there are two 
adjacent instantaneous propagators. In Refs. |l8],^5|, sign of the coupling of the vector 
mesons in the equations equivalent to Eq. ( |3.37p has the wrong sign; the coupling of 
the mesons in Eq. ( |3.37| ) must be the same as in Eq. ( |3.36| ) . 

Also note that in principle the same prescription has to be applied to the contact 
interaction, although to the order of two mesons, this simply has the effect of removing 
the primes. We find that 


Pi,2c ~ 


E . 

cx=7r,rj,a,S 


Sa 

M 


d 2 x ± dx p(x)[g n T 7T T n $ 7T (x)][g a T a T a $ a (x)\p(x). (3.38) 


Note that Pj 2c and P T 2 have forms that are very similar. In fact, we can obtain Pj 2c 
from Pf 2 by making the following changes: 


1. Replace ^ with 

2. Replace aq with 7 r. 

3. Restrict the sum on a 2 to 7 r, 77 , a, 5. 
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3.3.4 Interaction Hamiltonians in Momentum Space 


We look at the matrix element of the interaction Hamiltonians between initial and 
final states given by |&j, A*, t,) = Aj)x Ti |0) and (k f ,X f ,T f \ = {0\b(k f , X f )xt f 
where both the bra and the ket have units of [M -3 / 2 ]. We find that 


PTM’*)= E 


g a 2M 


d 2 q^dq + 6(q ] 


=^S, P ^2(27r) 3 y/kfl4 J y/2(2ir)3y/F 
u(k f ,X f )r a u(ki , A i)xl f TaXn 

[ d 2 x ± dx-e +ik f x »e~ ik ? x » [( i a (q)e~ iq ^ + al(q)e +iq ^] F a (q), 


(3.39) 


where q = ( q + , q ± ) and there are implicit ^-functions on q + for each particle to ensure 
that the light-front energy is positive. These are occasionally suppressed to simplify 
the equations. We include a meson-nucleon form factor F a to phenomenologically 
account for the fact that the mesons and nucleons are composite objects. 

Now, we evaluate P~ at x + = 0 and use 


I d 2 x ± dx-e< k f- ki - q ^ = 2(2vr ) 3 8 {2 ' + \k f -k t - q), (3.40) 


to write 

PFAf,i)= E 

a=n,r],a,s,p,uj y / 2(27r) 3 W kk+ 

[a a (q)6(k+ - k; ) - <J a (q)9(-(kf - kf))\ F a (q). (3.41) 

We define q = sign (kj — kf)[kf — ki ), which ensures that the q + of the meson is 
positive. 


Next we consider the two-meson interactions. First take the interaction with the 
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instantaneous propagator given by Eq. (|3.37|). Plugging in the field definitions 


gives 


-P/Tat/'O = E 


g ai g ao 2M f d 2 q 1± dqfd(qf) f d 2 q 2± dq^9(q^) 


2(2tt Y^kpt ' W 2 VH J (2vr)3/ 2 v^+ 

A / )r oa 7 + r« 1 «(fc i> \)x\T a 2 T aiXn 


d 2 x ± dx~e +ik f x ^e- ik ^ 


'a a 2 (q 2 )e~ iq ^ + at 2 (q 2 )e +i ^] F a 2 (q 2 ) 

(3.42) 


a Ql (qi)e-<^ +at tl (q 1 )e + ^J F ai { Ql ) 

y g ai g a2 2M ^ kf ' X f) V ^l + ^AK y_ 

2(2vr) 3 ^kfkf 

' d 2 k m j_dk+6(k+) 


ai -,OL2—n,r],a,S,p,uj 


X 


t T T 

Tf- L a2- L OL 1 


Xn 


zk+yyt 

\ a a 2 (kf — k m ) 9 (kj: — k m ) + o-Q 2 (&m — kf) 6 (k m — A^)] F a 2 ( y q 2 ) 
k(^ - - fc * + ) + - A;+)] ^(qq) ( 3 . 43 ) 


Note that the momenta qq and q 2 are implicitly functions of the momenta kf , k m , and 

k . 

We could also write out the contact interaction given by Eq. ( p.38|) in momentum 
space, but it is related to Eq. ( |3.43|) by replacing e< ' p 9p ' , + with and restricting the 
allowed values of the a’s. 


3.4 Feynman Rules for Nucleon-Nucleon Potentials 

Now that we have the one-meson- and two-meson-exchange expressed in momentum 
space, we are now ready to write out the Feynman rules for diagrams in our model. 
For simplicity, the only diagrams considered are those where a meson emitted by 
one nucleon is absorbed by the other nucleon. Since we are only interested in the 
two-nucleon to two-nucleon potentials, we follow the same approach as outlined in 
section |2.1.2| to derive the rules. We denote a “normal” nucleon propagator by a solid 
line with an arrow, an instantaneous nucleon propagator by a solid line with an stroke 
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across it, mesons of all type by a dashed line, and energy denominator terms by a 
vertical, light, dotted line. 


1 . 

2 . 

3. 

4. 

5. 


Overall factor of 


4 M 2 S ± ’+(P f -P i ) 
2(2 *) 3 S /k+k+k+k+ ' 


Usual light-front rules for p± and p + momentum conservation. 


Factor of for each internal line, including any instantaneous nucleon lines. 

A factor of ——^—— for each energy denominator. 

p -E» <ii 

Each meson connects the two nucleons, and each end of the meson line has a 
factor of g a T a T a F a (q). The indices of the isospin factors on each end of the 
meson are contracted together. The Lorentz indices of the gamma matrices are 
contracted with — g^ v for the vector mesons. 


6 . For each contact vertex, multiply by a factor of jj. If the vertex is a 7r — tt 
vertex, symmetrize the T n T n = TiTj by replacing it with Sij. 

7. Factor of = ^ A w(/c, A)u(fc, A) for each propagating nucleon and A- for an 
instantaneous nucleon. 


8 . Integrate with AM 2 f d 2 (-^p + over an F hiternal momentum loops. 

9. Put the spinor factors for nucleon 1 and 2 between mm’s and the isospin factors 
between the initial and final state isospin. 


From this list, it is useful to summarize what needs to be done to convert a graph 
with an instantaneous nucleon to one with a contact interaction: 

1. Replace e UU7 + J_. 


2. If both mesons are pions, replace T{Tj with S tJ . 

These rules make it easy to write down what various potentials are. 
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K = 


T 


-a— r ~ —7— -a- r 

\ / \ \ i 

+ x + \ + \ i 

/ x \ \i 

— i -. —.— 



+ ... 


Figure 3.1: The first several terms of the full kernel for the Bethe-Salpeter equation 
of the nuclear model with chiral symmetry. 


3.5 Nucleon-Nucleon Perturbative Potentials 


The meson exchange potentials have the same basic form as in Section [2.1.2| , however 
we must include the contact interaction and instantaneous nucleon propagators for 
the nuclear model used here. First, we discuss how to include the contact diagrams 
from the standpoint of the Bethe-Salpeter equation to maintain chiral symmetry, then 
we begin to calculate the light-front potentials. 


3.5.1 The Bethe-Salpeter Equation 

The kernel of the Bethe-Salpeter equation m for this nuclear model is richer than 
the one presented in section |2.1.4| for the Wick-Cutkosky model. This is due mainly to 
the presence of the contact interactions which are generated by the chirally invariant 
coupling of the pion to the nucleon. Several of the lowest-order pieces of the full kernel 
K are shown in Fig. ECU (Note that for Feynman diagrams, it is useful to combine 
the “normal” nucleon propagator with the instantaneous nucleon propagator |ff^ [74 
and denote the combination with a solid line.) 


As discussed in section |2.1.4| , each of these Feynman diagrams is covariant. This 
means that we may choose any of the diagrams from K to construct a new kernel K', 
and the infinite series of potential diagrams physically equivalent (Sec. |2.1.5|) to K' 
will also be covariant. 


In practice, this means that when deciding which two-meson-exchange potentials 
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to include for calculating the deuteron wave function, we may neglect the crossed di¬ 
agrams. Although including only the box and contact two-meson-exchange diagrams 
may affect the exact binding energy calculated, we should find a partial restoration 
of rotational invariance. We reiterate that the focus of this work is on understanding 
the effects of the breaking of rotational invariance and how to restore it; our goal is 
not precise agreement with experimental results. 


Chiral Symmetry 


We also want to keep the potentials chirally symmetric as well. Whereas Lorentz 
symmetry is maintained by using a kernel with any Feynman diagrams (with poten¬ 
tially arbitrary coefficients), chiral symmetry relates the strength of the mi contact 
interaction to the strength of the pion-nucleon coupling. 


Chiral symmetry tells us that for pion-nucleon scattering at threshold, the time- 
ordered graphs approximately cancel [[1^]. Furthermore, upon closer examination, we 
find that all the light-front time-ordered graphs for the scattering amplitude vanish 
except for the two graphs with instantaneous nucleons and the contact graph. These 
graphs are shown in Fig. [3.2| . Using the Feynman rules, and denoting the nucleon 
momentum by k and the pion momentum by q, we fold that 


Mu = C 2 ^ q+) u(k'hMk), (3.44) 

M x = C ^ q+) u(k')^u(k\ (3-45) 

M c = C-^-u(k')u(k), (3.46) 

where the factors common to all amplitudes are denoted by C. 


For threshold scattering, we take k = k' = M and q = q' = m n . In that limit, we 
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- \ | / - - r" | x . - -V- 

(a) ( b ) (c) 

Figure 3.2: The non-vanishing diagrams for pion-nucleon scattering at threshold: (a) 
Mu, (b) Mx, and (c) Me- The mesons here are pions. 


find 


Mu = C 


M x = C 


M c = C 


/ ^hj T ^i,j,k^~k 

2(M + m n ) 

2 (M — m n ) 

/ ~°i-3 


M ’ 


(3.47) 

(3.48) 

(3.49) 


where C = Cu{k')u{k). In the limit that m n —> 0, the sum of these three terms 
vanishes. The term in these equations proportional to Tk is the famous Weiuberg- 
Tomazowa term |75|.[76|. 

The fact that the amplitudes cancel only when the contact interaction is included 
demonstrates that chiral symmetry can have a significant effect on calculations. In 
terms of two-pion-exchange potentials, this result means that the contact potentials 
cancel strongly with both the iterated box potentials and the crossed potentials. This 
serves to reduce the strength of the total two-pion-exchange potential, which should 
lead to more stable results. 

However, since we do not use the crossed graphs for the nucleon-nucleon potential, 
we must come up with a prescription which divides the contact interactions into two 
parts which cancel the box and crossed diagrams separately. We do this by formally 
defining two new contact interactions, so that 


M 

M 


r = - M 

Cu 2 (M + m u ) 

c x — -Me — Mc v - 


c, 


(3.50) 

(3.51) 
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Figure 3.3: The two diagrams which contribute to the OME potential for each meson. 


With these definitions, we find that at threshold and in the chiral limit, 

Mu + M Cu = 0, (3.52) 

Mx + M Cx = 0. (3.53) 

This indicates that we can incorporate approximate chiral symmetry without includ¬ 
ing crossed graphs simply by weighting each graph with a contact interaction by a 
factor of 


M 


2(M J r m 7 r) * 

3.5.2 OME Potential 

The one meson exchange (OME) potential connects an initial state with two nucleons 
to a final state with two nucleons, has one meson in the intermediate state, and has 
the meson emitted and absorbed by different nucleons. With these restrictions, along 
with the fact that in light-front dynamics the interaction does not allow for particles 
to be created from the vacuum, we End that there each meson has only two diagrams 
for the OME potential. These diagrams are shown in Fig. [To. 


Note that chiral symmetry does not impose restrictions when considering the 
OME potential alone. It does affect the TME potentials considered, as discussed in 
the following section. 

The Feynman rules derived in the previous section are used to derive the potential 
for these diagrams. We factor out an overall factor of iM <5( -- ' +A ' ; = L, ~^~- that is 

2 \J k lf k li k 2f k 2i 
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common to all the two-nucleon potentials and suppress it from now on. Then we get 
n 2r r T 

jol~*~ a,l -*■ o 


Id 


a, 2 


OME,a 


(2tt) 


u(ki f , A 1/ )r a u(k u , Ai i)u(k 2f , A 2 f)r a u(k 2 i, A 2 i)F^(q) 

0(K f - k+i) 


(K } - K)( p - k u - k 2f) ~ m l~ - k if ,±) 2 

_ (KK - K f ) _' 

(K - Kf)(P- - Kf - k 2i) -ml- (ku 7± - k lf)± y 

To simplify the potential, we define 


a = [6{k+ - kfi) x (Kf ~ Ki)(P~ 

o(K-Kf) x (K-Kf)(P- 

b = 2 k it ±kf 7 ±, 


k li k 2 /) P 

Kf - K )] - k h 


- k 2 
K f,- L> 


(3.54) 


(3.55) 

(3.56) 


so 


^ _ 9 a Pa,l ' Pa,2 ( „\2^( k ^fi ^1 /)T 0 ,'u(fcij, A\i)ll{k 2 f , A2/)T Q , , u(fc2i, A 2 ?) 


/ OME,a 


^a(Q) 


. (3.57) 


(27t) 3 [a - ml + bcos((j)f - 4>i)] 

Now we consider the precise form to use for the meson-nucleon form factor. We 


assume that it has a n-polc form |5£|, so that the denominator of the meson-nucleon 
form factor has the same form as the denominator of the potential in the scattering 
regime. In particular, A a playing the role of m a for the form factor. For simplicity, we 
choose that the denominator of the form factor has the same form of the denominator 
of the potential. Thus, 

.2 


Fa(q) = 


A« - 


(3.58) 


7 a - Al + bcos((j)f - 4> i ) / 

Inserting the explicit expression of for the meson-nucleon form factor into the 
potential results in 

n 2 r r t 

T r _y a ± a,l J «,2 /A 2 ^,2\2n a 

'OME,a - -- A " _ ^ 

^( k ifi di/)T Q , , it(fcij, Xu)u(k 2 f, A 2 f)T a u{k 2i) A 2 j) 


[a - ml + bcos((j)f - </>*)] [a - A. 2 a + 6cos(0/ - </>*)] 


‘2tl 0 . 


(3.59) 
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Now, since the light-front potentials have exact rotational invariance about the 
z-axis, they conserve the J z quantum number m. Thus, these potentials connect only 
states with the same value of m. This means that, in general, light-front potentials 


may be written as 

OO 

V{<j>fAi)= J 2 (3.60) 

m=—o o 

where V m is the potential in the magnetic quantum number basis. This relation can 


be inverted to obtain V m in terms of V((pf, (pi), 
glT a ,i ■ T at 2 


V m 

^OME.q 


(2it) 

d(pf 


-(A 


l - <) 2n ° 




Aij)rAii) 14(^2/ 5 A2/)r a '?i(fc22 ? A 2 i)]^Q 


, „ „ . (3.61) 

J 2tt [a - m\ + bcos((p f )} [a - A 2 a + &cos(0/)] 2n “ 

To perform the integration, we need an expression for (pj dependence the uu matrix 
elements. These are calculated in Appendix |F]. Summarizing, we can write 


u(kif, X lf )r a u(k u , \ u )u(k 2 f, A 2 /)r Q u(fc 2 i, Am) = y^ Cj(T a , r a )e l3<l>f , (3.62) 

j 

where the C depends implicitly on all the variables on the left-hand side except (pf 
and (pi. 


Thus, 


AaME,ct — 9a{T a , 1 ‘ Ta,2 

f d( pf 


(Aq - rn 


2 \2tIq: 


(2tt ) 


£c j (r„,r«) 


a i(m+j)(j> f 


i2 riot 


(3.63) 


The cosine integral is denoted by /, and is 
for the integral, we obtain 

KjME.a = 9a(T a ,l ' T a ,2j ^ \ 3 


27T [a - m 2 a + bcos((pf)} [a - A 2 a + bcos((pf)]‘ 

calculated in Appendix 0. Substituting / 


( A a - m a 


2 \ 2n a 


(2 


,n 


Cj(r a , T a )I(m + j, a - ml, l, a- A 2 a , 2 n a , b ). 

n 


(3.64) 


All of the terms in Eq. ( |3.64| ) are known, and the potential can now be calculated 
numerically. 
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3.5.3 TME Potentials 


The two-meson-exchange (TME) potentials consider here are the box diagrams (see 
Fig. [).4| ) and the contact diagrams (see Fig. |3.5| ). We do not consider the crossed 
diagrams because they are not needed to restore rotational invariance, as shown in 
Section |3.5.1| . We derive a few TME potentials to illustrate the general form that 
they take. 


Stretched Box Potential 


We use the Feynman rules to write the stretched-box potential shown in Fig. |3.4| (b), 


T r a f>OLi _ JOLf^OLi 

K TME:SB — 


5 , a f 5 , a i 4M 2 [T Q/i 1 • T af)2 ][T ai ,i ■ T Qii2 ] f d 2 k 2m ±dk. 


2m 


i/(f r a/ - 


(2vr) 3 
hm + M 


2 M 


-r„. 


2(2tt fk+ m kt m 

hm + M^ 


\ X l( r ^2m + M v V 

/lm X 2/\ 1 a/ 2]\/[ a i I / 2m 


d(ktmW(kt m )0(k+ - k+je(kl m - k+) 


(ktm - Kf)( p ~ k i m ~k 2f )~ mi 2 


a f 


(k 


1 fl. ^lml 


P -k u - k 2f -q f -q i 

_ FM 2 _ 

{kti - ktJ(P~ ~ k~ - k 2m ) - m 2 . - (fe lm± - k li± ) 2 
+ {1^2}. (3.65) 


To compress notation, we defined i/(= u(kif,\if), )u = u(ku, Aij), and so on. The 
symbol {1 <->■ 2} means that all labels 1 are replaced with 2 and vice versa. This is a 
way of explicitly stating the potential is invariant under exchange of nucleons 1 and 
2 . 

We also used the following relation to simplify Eq. ( |3.65| ): 


A 




X T2f Ta f ,jXr 2m Xt2 m Tai,iXr 2i X.t\ f ^ otf.jX.Ti m X r\ rn ^ on,iXru 


A 


T n 


A 


— { r f\[T af ,i ■ T Q/ , 2 ][T Qiil • T aij2 ]|rj) 
= [T af ,i • T af>2 \[T au i ■ T a J. 


(3.66) 
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Figure 3.4: The TME potentials for (a) Vtme : m (the Mesa potential), (b) Vtme:Sb (the 
stretched box potential), (c) Vtme:Sbi (the stretched instantaneous potential), and (d) 
Vtme:Sbii (the stretched double instantaneous potential). Note that the graphs on the 
right side are obtained from the graphs on the left side by relabeling 1 <->• 2. 
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Figure 3.5: The TME potentials that include the contact interaction for (a) Vtme : c 
( the contact potential), (b) Vtme:Sbc (the stretched contact potential), (c) Vtme:Sbic 
(the stretched instantaneous contact potential), and (d) Vtme:Sbcc (the stretched 
double contact potential). Note that the graphs on the right side are obtained from 
the graphs on the left side by relabeling 1^2. 
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Mesa Potential 


We consider the Mesa potential next, Vtme:M, shown in Fig. |T4](a). This is one of 
the potentials which is due to the instantaneous nucleon, a feature not found in the 
Wick-Cutkosky model. Using the Feynman rules results in 


K TME:M 


9a f 9aA M2 [Ta f , i ■ T af )2 \[T aiA ■ T Qi)2 ] f d 2 k 2m ±dkt. 


2m. 


(2tt) 


2 i. 2 ^Y k t m k tm 


i f\ 


7 


(( r Q/ —r a . ))i m x 2 /(( r Q/ 


^2m + M 


2 M 


r «, > 


/2m 


- 4J^2m - k+ ) 

F af {q f ) 2 


( k tf ~ Hm)( P ~ k lf- k 2m) ~ m l f - ( fc l/± - fe l™i) 2 2M 


u l/ "'2mJ ""af 

FaMiT 


( k ti ~ k 2m){ P ~ k li - k 2m) ~ m l t ~ ( fc lm± ~ fcli±) 2 

+ {!<-> 2}. (3.67) 


Contact Potential 

The last potential which we calculate explicitly is the contact potential Vtme:C, shown 
in Fig. P75|(a). The rules relate this potential to Vtme:M in a simple fashion. If both 
the initial and final mesons are not pions, then the potential is 

i ' T af: 2][T au i • T Qi , 2 ] f d 2 k 2ml _dk% m 


\ 7 ’ a f’ a i _ 

V TME:C — 


1/ 


( ^r a/ g )im x 2f( 


(2vr) 3 J 2(2vr) 3 k+ m M 

h m + M, 


2 M 


-r Ql ) 


/2m 


e( k tje( k + - k L)o( k L - k t f ) 

F *f{<lfY 


( k 2f ~ k 2m)(P ~ k lf~ k 2m) ~ »Y f ~ (k lf± ~ k lm± ) 2 2 M 
_ F ai {q-f _ 

( k 2i ~ k tm)( P ~ ~ k li ~ k 2m ) “ "Y H ~ ( fe l m± ~ k liA _) 2 

+ {!<—> 2}. 


(3.68) 


To get the contact potential for the pions, the following change has to be made: 


T aj ,i ■ T a 2 }[T au i • T a . i2 \ —> t 2 — 3. 


(3.69) 
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This is due the additional symmetry that the two pion vertex has. 

Before this potential is used in calculations, it must be multiplied by a factor of 
2 (M+m ) ( as discussed in Section |3.5.1| ) if the crossed diagrams are not included. 


Common Features and Simplification 


Each potential in Figs. |3.4| and |3.5| can be written schematically in the following, 
general form: 


Vtme — 


gHM 2 T\A 2 - m 2 ) 4n “ f d 2 k 2m± dk+ m 


(2tt) 3 J 2(2tt) 3 

lf{(T af MlT ai )) lm X 2f{(^'a r M. 2 T ai )) 2m 




a f~ m l f + b f cos (0/ - 0 m ) 


~ n f 


2 Tin. 


a f ~ K f + b f COs(0/ - 0m) 

[a m T b mf COS(0J 0m) T b m i COS(0m 0i)] 

[a* - + bi cos(0 m - fa)] [a, - A 2 . + bi cos(0 m - 0*)] ^ , (3.70) 


where the expression for F a in Eq. ( |3.58| ) was used. 

The TME potentials, like the OME potential, have exact rotational invariance 
about the z-axis, and thus conserve m. We project the potential onto states of 
definite m by setting fa to zero and integrating the potential with f dfife 1 ™^*/2n. 

Now we have to obtain the uu matrix elements explicitly to perform the azimuthal 
integrations. Appendix [F] shows that we can write 

Q-i ) ) 1 m X 2 f ((r af M 2 r ai )) 2m ] ^ =0 

= ^ Cj >k (T af ,Mi, T ai ; r af ,M 2 , r (3.71) 

j,k 

Note that the C is implicitly a function of the momentum and helicity of the initial 
and final states, but is independent of the azimuthal angles 0/ and fa. 

Then, after a change of variables, 0 m —> 0', 0/ — 0 m —> 0j, and 0/ —> 0)- + 0', the 
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potential can be written as 



f( k lnv k 2m,q-fi<lt) 


f 



^ f 2 * dftj c i(m+jW f c i(m+k)ti 


l 


—2 n t 


f 


(3.72) 


The azimuthal-angle integrals taken together are denoted as I, and the method for 
calculating the (p' m and <f>'f integrals is discussed in Appendix Summarizing, we 
can write 



f( k lnv k 2rm<l}i<lt) 


(3.73) 


This potential is evaluated using the numerical integration techniques discussed in 
Appendix [C]. 

3.5.4 Further Development of the Light-front Schrodinger Equation 

The potentials derived here possess a high degree of symmetry. To solve the light- 
front Schrodinger equation efficiently, these symmetries should be explicitly exploited, 
as was done in section |2.1.3| . In addition to the invariance of the potentials under 
parity, there are additional symmetries due to the conservation of nucleon helicity 
and invariance under time reversal [Q . 
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We follow Machleidt’s approach for taking advantage of the symmetry of rotation- 
ally invariant potentials with helicity |59|]. However, since the light-front potentials 
derived here do not have full rotational invariance, Machleidt’s approach must be 
modified. The symmetry properties of helicity matrix elements are rederived in Ap¬ 
pendix |G] without assuming full rotational invariance of the potentials. These results 
allow a modified version of Machleidt’s approach to be combined with the exploita¬ 
tion of parity (using the transformation from light-front coordinates to equal-time 
coordinates) discussed in section |2.1.3| . In particular, the potentials are initially cal¬ 
culated in the |pet, 6, M, X±, X 2 ) basis, although the relations in Appendix |G.5| are 
used to transform to the |p ET , J, M, L, S) basis, which is useful to solve for the wave 
functions. Note that the potentials connect states with different J values in general. 
Once the symmetries have been explicitly expressed, we can discretize the Schrodinger 
equation as done in Appendix [Bj. 

Note that we apply the transformations to the potential in order to simplify the 
calculation of the wave functions. Once the wave functions are obtained, we may 
apply the inverse of the transformations to the wave functions to express them in the 
helicity basis (|pj_,p + , Ai, A 2 ) ) or, by using Eq. ( |G.111|) , the Bjorken and Drcll spin 
basis ( \p ± ,p + ,m 1 ,m 2 )). 


3.6 Results for the Nucleon-Nucleon Potentials 


The next step towards numerically calculating the bound states for these potentials is 
to choose the parameters (meson masses, coupling constants, etc.) for the potentials. 
We consider two models: a pion-only model and a light-front nucleon-nucleon model 
which incorporates all six mesons used in the Refs. [H|, 55 . 
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3.6.1 Pion-only Results 


As a starting point, we look for inspiration from the non-relativistic one-pion-exchange 
(OPE) model used by Friar, Gibson, and Payne (FGP) [Q. They use the conventional 
OPE potential along with a n-pole pion-nucleon form factor to regulate the high- 
momentum (short-range) part of the potential. By tuning the form factor parameters 
to reproduce the physical binding energy of the deuteron, they find that the other 
the deuteron observables are also reproduced fairly accurately. This leads to the 
conclusion that this OPE potential model is a good starting point for investigating 
the deuteron observables. 

The basic parameters of the FGP model are the mass of the nucleon (M = 
938.958 MeV), the mass of the pion (m, = 138.03 MeV), and the pion decay constant 
(/o = 0.079), which corresponds to a coupling constant of |^ = L ^Y° = 14.6228 in 
our formalism. The family of ?r-pole pion-nucleon form factors they consider have the 
form 



A 2 — m 2 
A 2 + q 2 


n 


(3.74) 


The parameter A is fit to reproduce the deuteron binding energy for a given value of 

n. 


We can obtain the FGP model from our light-front nucleon-nucleon potential by 
considering only pion exchanges and performing a non-relativistic reduction. To un¬ 
derstand our results, it is important to first consider how well the light-front potential 
and wave functions approximate the non-relativistic potential and wave functions. 

It is useful to start by reviewing some of the properties of the deuteron. First, the 
deuteron is very lightly bound. Second, although a majority of the deuteron wave 
function resides in the non-relativistic regime, it has a high-momentum tail that falls 
off rather slowly [Q. Recalling our experience with the Wick-Cutkosky model, where 
we found that mass splitting between states with different m values is small for lightly 
bound states, it may appear plausible that masses of the rn — 0 and m = 1 states of 





the deuteron should be approximately the same when calculated with the light-front 
OPE potential. However, the high-momentum tail of the deuteron enhances the effects 
of the potential’s relativistic components. Since the breaking of rotational invariance 
is a relativistic effect, as shown in Eq. (|2.63| ), this implies that mass splitting may be 
large for the deuteron calculated with the light-front potential. 

To clearly understand how a large mass splitting might arise in the light-front pion- 
only model, we take a step back and consider a scalar version of the pion-only model, 
in which we assume that the pion has a scalar coupling to the nucleon. This allows for 
a more direct comparison with the Wick-Cutkosky model, since the main difference 
between the scalar-pion-only potential and the Wick-Cutkosky potential is the pion- 
nucleon form factor. Since the denominator of the pion-nucleon form factor has the 
same form as the denominator of the potential, the form factors do not significantly 
change the rotational properties of the scalar-pion-only potential. Another difference 
is that factors of uu appear in the numerator of the scalar-pion-only potential, but 
the effect of these terms is small. 


In Table |3.1| , we show the binding energies for deuterons calculated with the scalar- 
pion-only model. Two pion-nucleon form factors are considered, the first one has 
A = 1.0 GeV, for which the coupling constant was fit to give the correct binding energy 
for the non-relativistic potential, and the second form factor uses A = 1.915 GeV, 
which was fit to give the correct binding energy for the light-front potential. Those 
form factors are used to calculate the binding energies for the non-relativistic and 
light-front potentials, as well as the instantaneous and retarded potentials, which are 
relativistic and defined by analogy with Eqs. ( |2.81|) and ( |2.83| ). We find that the 
binding energies for all of the potentials have the same order of magnitude, and that 
the light-front potentials have consistently lower binding energies, which confirms the 
behavior observed in Fig. [2.9| . In addition, the binding energies of the m — 0 and 
rri — 1 light-front potentials are essentially degenerate for both form factors. 

Now we ready to consider the (pseudoscalar) pion-only model. The only differ- 








Table 3.1: Binding energies for deuterons calculated with several potentials and two 

different pion-nucleon form factors for the scalar-pion-only model. The n parameter 

2 

of the form factor is 1, the coupling constant is f- = 0.424, and the pion mass is used 
as the mass of the exchanged meson. 


Potential 

A = 1.0 GeV 

A = 1.915 GeV 

Non-relativistic 

2.2236 MeV 

4.2539 MeV 

Instantaneous 

2.1581 MeV 

4.0862 MeV 

Retarded 

2.3111 MeV 

4.5224 MeV 

Light-front, rn — 0 

1.2027 MeV 

2.2296 MeV 

Light-front, m = 1 

1.2027 MeV 

2.2294 MeV 


ence between this potential and the scalar-pion-only potential is that the numerator 
contains factors of ui^u instead of uu. Although this is formally a small change, it 
has a large effect on the binding energies. The pseudoscalar coupling generates a ten¬ 
sor force, which is more sensitive to the relativistic components of the wave function 
than the scalar force. In general, this means that the differences in binding energies 
obtained with different potentials, such as those shown in Table m will be larger for 
the pion-only potential. 


In Table |3.2| , we show the binding energies for deuterons calculated with the 
pion-only model. Two pion-nucleon form factors are considered, the first one has 
A = 1.01 GeV, which was fit for the non-relativistic pion-only model [5(|, and the 
second form factor uses A = 1.9 GeV, which was fit to give the most reasonable 
binding energy for the light-front potentials. We find that the binding energies vary 
greatly depending on which potential is used. In fact, the light-front potentials do 
not bind the deuteron with the first form factor, and with the second form factor, 
the mass splitting between the m = 0 and m = 1 states is very large. These facts 
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Table 3.2: Binding energies for deuterons calculated with several potentials and two 

different pion-nucleon form factors for the (pseudoscalar) pion-only model. The n 

2 

parameter of the form factor is 1, the coupling constant is |- = 14.6, and the pion 
mass is used as the mass of the exchanged meson. The negative binding energy for 
the light-front potentials in the first column indicates that the states are not bound. 


Potential 

A = 1.01 GeV 

A = 1.9 GeV 

Non-relativistic 

2.2244 MeV 

227.019 MeV 

Instantaneous 

0.0146 MeV 

28.192 MeV 

Retarded 

1.3590 MeV 

130.472 MeV 

Light-front, m = 0 

-0.0269 MeV 

0.788 MeV 

Light-front, rn — 1 

-0.0246 MeV 

8.856 MeV 


indicate that the deuteron wave functions are very sensitive to subtle changes in the 
relativistic structure of the pion-only potentials. 

Now that we understand that relativistic effects are important for the OPE po¬ 
tential, we can start to analyze the TPE potential. We reiterate here that the main 
objective of this work is to investigate this breaking of rotational invariance and how 
it is restored. Obtaining results which are in agreement with experimental data is 
of lesser importance; we must have states that transform correctly under rotations 
before we can address the experimental data. Naively, one might think that we can 
choose any set of two-pion-exchange graphs which is a truncation of a rotationally- 
invariant infinite series of graphs. In particular, one might think that using the box 
graphs shown in Fig. [Tl| (the two-pion-exchange graphs without the graphs mandated 
by chiral symmetry, denoted as the non-chiral-two-pion-exchange (ncTPE) potential) 
would be adequate for our analysis of rotational invariance. 

Another choice of a set of two-pion-exchange graphs is sum of the box graphs 
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shown in Fig. |3.4| and the contact graphs shown in Fig. |3.5| . To incorporate chiral 
symmetry as accurately as possible without including the crossed two-pion-exchange 
potentials, we weight the contact vertex with a factor of 2 (M+ m ) , as ex pl a i n ed in 
Section |3.5.1| . Note that in the sum, which we call the two-pion-exchange (TPE) 
potential, chiral symmetry provides a cancelation between the box diagrams and the 
contact diagrams. This indicates that the results obtained with the TPE potential 
will be more stable than those obtained with the ncTPE potential. 

To check this stability and the restoration of rotational invariance of the pion- 
only model, we calculate the energy of the deuteron using the OPE, OPE+TPE, and 
OPE+ncTPE potentials. We verify that the results are independent of the choice of 
the pion-nucleon form factor by considering three different choices for the form factor. 

For the Erst form factor, we choose n = 1 and hnd that A = 1.9 GeV gives a 
reasonable range of binding energies. We use the form factor to calculate the lowest 
energy bound state with arbitrary total angular momentum. In addition, we take the 
J = 1 component of the potential and use it to perform the same calculation. We 
expect to obtain similar results since the deuteron is a spin one object. 

The results for the first pion-nucleon form factor are shown in Table |3.3| . Note 
that the results do not depend much on whether a restriction to spin one is applied. 
This shows that even though the eigenstates of the light-front Hamiltonian are not in 
principle eigenstates of the angular momentum, numerically they are almost angular 
momentum eigenstates. However, this does not mean that rotational invariance is 
unbroken. Rotations relate the difierent m states, and from the mass splittings, we 
see that the states do not transform correctly. 


Table |3.3| also shows that both the mass splitting and the difierence in the per¬ 
centage of the D-state wave function decrease when TME diagrams are included. In 
addition, the wave functions are almost completely in the J = 1 state. As for the 
D-state probability, we first note that it increases with the binding energy. It is con¬ 
sistent (given the range of values it and the binding energy take) with the value of 
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Table 3.3: The values of the binding energy, percentage of the wave function in the 
D state, and the percentage of the wave function in the J — 1 state for the m = 0 
and m — 1 states for different potentials. The values obtained for the full potential 
(which connects states of arbitrary total angular momentum) are shown along with 
the values obtained for potentials restricted to the J = 1 sector which are shown in 
parentheses. The pion-nucleon form factor uses n = 1 and A = 1.9 GeV. 


Potential 

Binding Energy (MeV) 

D state (%) 

J = 1 (%) 


m=0 

m=l 

Diff 

m=0 

m=l 

m=0 

m=l 

OPE 

-0.7884 

(-0.7129) 

-8.8561 

(-7.6970) 

8.0677 

(6.9841) 

4.09 

(3.97) 

12.27 

(11.99) 

99.99 

(100) 

99.78 

(100) 

OPE 

-0.6845 

-2.5606 

1.8761 

3.98 

8.08 

99.98 

99.88 

+TPE 

(-0.5886) 

(-2.0644) 

(1.4758) 

(3.79) 

(7.61) 

(100) 

(100) 

OPE 

0.0107 

0.0087 

0.0020 

0.15 

0.66 

100.00 

100.00 

+ncTPE 

(0.0140) 

(0.0155) 

(-0.0015) 

(0.11) 

(0.58) 

(100) 

(100) 


6% reported in Ref. [Q for the FGP model. 

Another thing to note in Table [P3] are the effects of using the ncTPE potential 
is used. As mentioned earlier, this potential is greater in magnitude than the TPE 
potential and should have a larger effect on the binding energy. In fact, the effect 
is so large that it serves to unbind the deuteron. (Strictly speaking, this indicates 
only that the binding energy is very small or zero; the error in the binding energy 
calculation increases as the binding energy approaches zero.) Because the ncTPE 
potential has such a large effect by itself, it is impossible to determine what effect it 
has on the rotational properties of the state. Only the TPE potential can be used to 
analyze the restoration of the state’s rotational invariance. 


To make sure that the results found are independent of the pion-nucleon form 
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Table 3.4: A pion-nucleon form factor with n 


1 and A = 2.1 GeV is used to calculate 


the same quantities as shown in Table [O 


Potential 

Binding Energy (MeV) 

D state (%) 

J = 1 (%) 


m=0 

m=l 

Diff 

m=0 

m=l 

m=0 

m=l 

OPE 

-1.7357 

(-1.6103) 

-14.9116 

(-13.2111) 

13.1759 

(11.6008) 

5.60 

(5.49) 

14.30 

(14.14) 

99.98 

(100) 

99.71 

(100) 

OPE 

-1.6296 

-4.5609 

2.9313 

5.62 

9.89 

99.96 

99.85 

+TPE 

(-1.4551) 

(-3.8203) 

(2.3652) 

(5.43) 

(9.53) 

(100) 

(100) 

OPE 

-0.0219 

-0.0180 

-0.0039 

0.69 

1.33 

99.99 

99.99 

+ncTPE 

(-0.0045) 

(0.0056) 

(-0.0101) 

(0.38) 

(0.86) 

(100) 

(100) 


factor, we performed the same calculation with n = 1 and A = 2.1 GeV (shown in 
Table |3.4|) and with n = 2 and A = 2.9 GeV (shown in Table |3T5|) . The results in these 
two tables are qualitatively the same as the results in Table |3.3| , which demonstrates 
that these results are robust. 


The deuteron masses for the different m states, as shown in Tables |T3|, fh4|, and 


are shown in Fig. |3.6| as a function of which potential was used for the calculation, 
ft is easy to see that the relatively large mass splitting obtained with the OPE po¬ 
tential is significantly reduced by including the TPE potential. The apparent further 
reduction of the splitting by neglecting the chiral part of the TPE potential (leaving 
the ncTPE potential) is an artifact of the deuteron becoming unbound. Since the 
unbound states form a continuum, the lowest energy unbound states will always be 
essentially degenerate, regardless of how badly rotational invariance is broken. 
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Pion-Only Model Mass Splittings 


n=1,A=1.9GeV n=1,A=2.1 GeV n=2, A=2.9 GeV 



Figure 3.6: The values of the binding energy for the m — 0 and m — 1 states for 
different pion-only light-front potentials. Three different pion-nucleon form factors 
are used: (a) n — 1 and A = 1.9 GeV, (b) n — 1 and A = 2.1 GeV, and (c) n — 2 and 
A = 2.9 GeV. 
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Table 3.5: A pion-nucleon form factor with n = 2 and A 


2.9 GeV is used to calculate 


the same quantities as shown in Table [O 


Potential 

Binding Energy (MeV) 

D state (%) 

J = 1 (%) 


m=0 

m=l 

Diff 

m=0 

m=l 

m=0 

m=l 

OPE 

-1.2060 

(-1.1068) 

-10.7746 

(-9.4180) 

9.5686 

(8.3112) 

4.86 

(4.74) 

13.06 

(12.82) 

99.98 

(100) 

99.75 

(100) 

OPE 

-1.0781 

-3.3983 

2.3202 

4.78 

8.95 

99.97 

99.87 

+TPE 

(-0.9494) 

(-2.7863) 

(1.8369) 

(4.59) 

(8.52) 

(100) 

(100) 

OPE 

-0.0002 

-0.0028 

0.0026 

0.30 

0.90 

100.00 

99.99 

+ncTPE 

(0.0069) 

(0.0106) 

(-0.0037) 

(0.20) 

(0.68) 

(100) 

(100) 


3.6.2 Light-Front Nucleon-Nucleon Model Results 


Having addressed the simple pion-only model in the previous section, we are now 
ready to consider the full nuclear model where the nucleon-nucleon interaction is 
mediated by all the six mesons shown in Table fi.6| . For numerical work, use the 
parameters for the light-front nucleon-nucleon (LFNN) potential from the work of 
Miller and Machlcidt [^5J. Those parameters were fit for a potential that used a 
retarded propagator for the energy in the potentials. Since the potentials used in 
this paper have energy dependent denominators, the parameters must be modified 
somewhat. We choose to vary the coupling constant for the a meson. The parameters 
are given in table [3.6| . 

As with all the other deuteron models presented in this paper, the light-front 
one-meson-exchange (OME) potential breaks rotational invariance and causes a mass 
splitting of the deuteron states with different magnetic quantum numbers. We expect 
that the splitting will be removed somewhat by including higher order potentials. 


The first step is to determine which two-meson-exchange potentials to use. One 
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Table 3.6: The parameters for the ir, rj, p, ui, S, and a mesons from Ref. |p|. For 
the meson type, “iv” and “is” stand for isovector and isoscalar, while “ps”, “v”, 
and “s” stand for pseudoscalar, vector, and scalar. The cr meson is also known as 
/o(400 — 1200), and the 6 meson is the ao(980). 


Meson 

type 

mass [MeV] 

A [GeV] 

n 

dl 

47r 

/ 

9 

7T 

iv,ps 

138.04 

1.2 

1 

14.0 


V 

is,ps 

547.5 

1.5 

1 

3.0 


p 

iv,v 

769. 

1.85 

2 

0.9 

6.1 

UJ 

is,v 

782. 

1.85 

2 

24.5 

0.0 

5 

iv,s 

983. 

2.0 

1 

2.0723 


a 

is,s 

550. 

2.0 

1 

f a x 8.9602 



choice is to use only the two-pion-exchange potentials, TPE and ncTPE, as defined 
in the previous section. However, we expect to get better results using the two- 
meson-exchange diagrams generated by all the available mesons, including the contact 
diagrams for the pions, which we denote as the two-meson-exchange (TME) potential. 
In addition, we can also investigate the effect of leaving out the contact potentials for 
the pions, resulting in the non-chiral two-meson-exchange (ncTME) potential. 

We do not include diagrams with a contact interaction between the nucleon, a 
pion, and another meson. This is because, as mentioned in section |3.5.1| , the infinite 
series of the box diagrams is rotationally invariant and the contact diagrams are not 
needed to achieve rotational invariance. Furthermore, they are not required to control 
the convergence of the series, since there is no strong cancellation between the contact 
diagram and the instantaneous diagrams. 

The first step in analyzing the bound states is to determine what range of / CT 
gives reasonable results. We iteratively solve for the binding energy of the deuteron, 
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Table 3.7: The values of f a required to give the physical value of the deuteron mass 
for a given potential and state with J z = m. The percent of the wave function in the 
D-state and in the J = 1 state are also shown. 


Potential 

u 

% D state 

%J= 1 


m=0 

m=l 

Diff 

m=0 

m=l 

m=0 

m=l 

OME 

1.2407 

1.2125 

0.0282 

2.87 

3.55 

99.99 

99.97 

OME 

+TPE 

1.2829 

1.2819 

0.0010 

2.96 

3.23 

99.99 

99.97 

OME 

+TME 

1.2968 

1.3079 

-0.0111 

2.95 

3.28 

99.99 

99.96 

OME 

-fncTPE 

1.3064 

1.3121 

-0.0057 

2.99 

3.16 

99.98 

99.97 

OME 

+ncTME 

1.3198 

1.3397 

-0.0199 

2.98 

3.21 

99.98 

99.96 


varying f„ until the binding energy matches the physical value of the binding energy, 
for each of the potentials. The results are shown in Table We find that a value of 
f a in the range 1.2 to 1.3 will give reasonable results for the binding energy. Note that 
D-state probability (about 3%) is lower in this model than for the energy-independent 
light-front used in Ref. |5I|, were a value of 4.5% is found. This is expected since the 
f a is greater than 1 in this model, meaning that the scalar interaction is strengthened 
relative to the tensor interaction, leading to a decrease in amount of the D-state 
present. 


We choose two values of f a , one from the low end of the range (1.22) and one from 
the high end (1.2815) for our investigations. Using two values helps ensure that our 
results are robust. 
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First, we examine the bound states for f a = 1.22. The results for several different 


choices of the TME potentials are shown in Table |3.8| and the binding energies are 
plotted in Fig. [J/7|. In addition to the two-meson-exchange potentials mentioned 
above, we also consider the ir-a plus 7r-u; Mesa potential. The reason for considering 
this potential is that Carboncll, Desplanques, Karmanov, and Mathiot (T^jj have shown 
that it helps restore rotational invariance of the deuteron. 

In particular, they have used manifestly covariant light-front dynamics to analyze 
the deuteron. They start with a deuteron wave function calculated in equal-time dy¬ 
namics, then use a light-front one-pion-exchange potential (expanded to lowest order 
in powers of jj) to calculate the perturbative corrections to the wave function. They 
find that the resulting wave function has an unphysical dependence on the orientation 
of the light-front plane, which would manifest itself as a breaking of rotational invari¬ 
ance in our formalism. They also use the n-o and n-u) Mesa potentials (expanded to 
lowest order in powers of jj) to calculate the correction to the wave function. When 
the wave function corrections are combined, they find that the directional dependence 
of the longest-range part of the deuteron wave function cancels exactly. 

This implies that for our model, using the n-cr plus tt-lu Mesa potential (which we 
denote by ir-(cr-uj)) should partially restore the rotational invariance of the deuteron, 
assuming that the breaking of rotational invariance is due primarily to the one-pion- 
exchange potential. Note that since we solve for the deuteron wave function self- 
consistently and to all orders for our potentials, we do not expect to find exactly the 


same result as Ref. 13 


The first thing to notice about the data in Table |3]8| is that the results are es¬ 
sentially the same regardless of if arbitrary angular momentum or is used or if the 
potential is restricted to the J = 1 sector. This was also seen in the previous section 
for the pion-only model. It means that the wave functions are numerically approxi¬ 
mate to angular momentum eigenstates. 

Next we notice the splittings between masses and D state percentages for the 
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m = 0 and m = 1 states. The dependence on m implies that the states do not 
transform correctly under rotations. All of the two-meson-exchange potentials used 
reduce the splittings by similar amounts, by about 60% for the binding energy and 
by about 70% for the percent D state. Note also that the mass splittings for the 
pion-only model were much larger. 

Examining the effects of the individual two-meson-exchange potentials, we see that 
TT-(a-uj)) potential does reduce the mass splitting, but it does not fully remove it. This 
is expected since the OME potential includes more than just the pion potential, and 
the potential is relativistic. 

Next, we compare the ncTME and ncTPE potentials to the TME and TPE po¬ 
tentials. The non-chiral potentials reduce the binding energy more than the chiral 
potentials, as we expected from our experience from the pion-only model, ffowever, 
unlike for the pion-only model, we find that the chiral and non-chiral potential have 
fairly similar effects. 

Finally, notice that the mass splitting for the TPE potential is much smaller than 
for the other two-meson-exchange potentials. By itself, this does not imply that the 
rotational properties of the deuteron calculated with that potential are significantly 
better than those from other two-meson-exchange potentials. The individual poten¬ 
tials that make up the TME potential are fairly large in magnitude, but vary in sign. 
This means that using any subset of those potentials may result in either a larger 
or smaller mass splitting. In this case, it is smaller. To investigate this further, we 
examine the currents for the TME and TPE deuteron wave function in Chapter [|. 

To verify that our results are independent of the value of / CT , we recalculate the 
deuteron properties for each of the potentials with f a = 1.2815. The results are 
summarized in Table [B~Tj| . and the binding energies are shown in Fig. [1 .8| . The change 
in f a increases the binding of the states, but the rest of the results are qualitatively 
the same. 
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Table 3.8: The values of the binding energy, percentage of the wave function in the 
D state, and the percentage of the wave function in the J — 1 state for the m — 0 
and m — 1 states for different potentials. The values obtained for the full potential 
(which connects states of arbitrary total angular momentum) are shown along with 
the values obtained for potentials restricted to the J = 1 sector which are shown in 
parentheses. The cr coupling constant factor is f a = 1.22. 


Potential 

Binding Energy (MeV) 

D state (%) 

J = 1 (%) 


m=0 

m=l 

Diff 

m=0 

m=l 

m=0 

m=l 

OME only 

-1.7653 

(-1.7141) 

-2.4200 

(-2.2655) 

0.6547 

(0.5514) 

2.73 

(2.72) 

3.61 

(3.54) 

99.99 

(100) 

99.96 

(100) 

OME 

-1.9236 

-1.7021 

-0.2215 

2.80 

3.38 

99.99 

99.96 

+7 r-(a-oj) Mesa 

(-1.8625) 

(-1.5340) 

(-0.3285) 

(2.78) 

(3.25) 

(100) 

(100) 

OME 

-0.4948 

-0.2646 

-0.2302 

1.97 

1.76 

99.99 

99.98 

+ncTME 

(-0.4611) 

(-0.2034) 

(-0.2577) 

(1.93) 

(1.59) 

(100) 

(100) 

OME 

-0.6620 

-0.4825 

-0.1795 

2.16 

2.09 

99.99 

99.98 

+ncTPE 

(-0.6232) 

(-0.4132) 

(-0.2100) 

(2.13) 

(1.99) 

(100) 

(100) 

OME 

-0.7861 

-0.6060 

-0.1801 

2.25 

2.31 

99.99 

99.97 

+TME 

(-0.7460) 

(-0.5191) 

(-0.2269) 

(2.22) 

(2.19) 

(100) 

(100) 

OME 

-0.9981 

-0.9155 

-0.0826 

2.42 

2.57 

99.99 

99.98 

+TPE 

(-0.9531) 

(-0.8253) 

(-0.1278) 

(2.39) 

(2.48) 

(100) 

(100) 
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Figure 3.7: The values of the binding energy for the m — 0 and m — 1 states for 
different nucleon-nucleon light-front potentials. The a coupling constant factor is 

fa = 1 . 22 . 


LFNN Model (f a =1.22) 



OME jt-(a-a>)Mesa ncTME ncTPE TME TPE 

Potential Used 
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Table 3.9: The a coupling constant factor used here, f a = 1.2815, distinguishes this 


table from Table 


Potential 

Binding Energy (MeV) 

D state (%) 

J = 1 (%) 


m=0 

m=l 

Diff 

m=0 

m=l 

m=0 

m=l 

OME only 

-3.3500 

(-3.2818) 

-4.4546 

(-4.2622) 

1.1046 

(0.9804) 

3.09 

(3.06) 

3.97 

(3.92) 

99.99 

(100) 

99.96 

(100) 

OME 

-3.6331 

-3.2408 

-0.3923 

3.10 

3.85 

99.99 

99.95 

+7r-(cr-u;) Mesa 

(-3.5520) 

(-3.0194) 

(-0.5326) 

(3.10) 

(3.77) 

(100) 

(100) 

OME 

-1.3766 

-0.9901 

-0.3865 

2.67 

2.64 

99.99 

99.97 

+ncTME 

(-1.3230) 

(-0.8807) 

(-0.4423) 

(2.65) 

(2.54) 

(100) 

(100) 

OME 

-1.6532 

-1.4693 

-0.1839 

2.81 

2.88 

99.99 

99.97 

+ncTPE 

(-1.5939) 

(-1.3578) 

(-0.2361) 

(2.79) 

(2.81) 

(100) 

(100) 

OME 

-1.8617 

-1.6032 

-0.2585 

2.85 

3.05 

99.99 

99.96 

+TME 

(-1.8021) 

(-1.4696) 

(-0.3325) 

(2.83) 

(2.97) 

(100) 

(100) 

OME 

-2.1915 

-2.2137 

0.0222 

2.95 

3.23 

99.99 

99.97 

+TPE 

(-2.1267) 

(-2.0831) 

(-0.0436) 

(2.94) 

(3.17) 

(100) 

(100) 
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LFNN Model (f o =1.2815) 



Potential Used 

Figure 3.8: The values of the binding energy for the rn — 0 and m — 1 states for 
different nucleon-nucleon light-front potentials. The a coupling constant factor is 
fa = 1-2815. 
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Chapter 4 

FORM FACTORS OF THE DEUTERON 

“It is theory that decides what can be observed.” 

—Albert Einstein 

In Chapter |3], we considered several different truncations of the light-front nucleon- 
nucleon (LFNN) potential and used them to solve for the deuteron wave function. In 
this chapter, we use those wave functions to solve for the deuteron current, which is 
used to calculate the deuteron electromagnetic and axial form factors. We can analyze 
the rotational properties of the deuteron current like we did for the wave function. 

In this chapter, we first outline the covariant theory of the electromagnetic form 
factors for spin-1 objects, like the deuteron. Then we recall the features of light-front 
calculations (including the breaking of rotational invariance) of the form factors. After 
that, we review the covariant and light-front tools for calculating axial form factors. 
The formalism is then applied to calculate the electromagnetic and axial currents and 
form factors for the light-front deuteron wave functions. 

One notable feature of this calculation is that it is done entirely with light-front 
dynamics. The covariant Lagrangian generates light-front potentials, which generate 
light-front wave functions, which are used in a light-front calculation of the deuteron 
current and form factors. This is different from other approaches which use deuteron 
wave functions calculated from equal-time dynamics, then transformed to the light 


front [77-84 . 


4.1 Electromagnetic Form Factors 
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Figure 4.1: The Feynman diagram for one-photon-exchange electron-deuteron scat¬ 
tering. 


4-1.1 Covariant Theory 

In the one-photon-exchange approximation, shown in Fig. |4.1| . the amplitude of the 
scattering process ed —* ed is just the contraction of the electron and deuteron cur¬ 
rents, multiplied by the photon propagator, 

(P r , \)^(k',m'\JX\k,m), (4.1) 

where 


(p'.A'KIp.A) = e«(p',A') 7 „M(p,A). 


(4.2) 


From Lorentz covariance, current conservation, parity invariance, and time rever¬ 
sal invariance [|85| -|89|1, we infer that the deuteron form factor can be written as 


(k\m'\J^\k,m) = - — (e*) p (fc', m!) J^e a (k, m), 


where the spin-1 polarization vectors satisfy 


e*(fc,m)e p (fc,m') = -d nhm >, 
^2e* fl (k,m)e u {k,m) = - g fiv + 

m 

k fJj e tl (k , m) = 0, 


kfjk v 

My 


(4.3) 

(4.4) 

(4.5) 

(4.6) 
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and the operator is given by 

]% = (K + M 


vUti 2 ) - f^tf) 


Cl" G M 2 ). 


(4.7) 


where = g%g% — g p g^ is the generator of infinitesimal Lorentz transformations, 
q = k! — k, and iq, F 2 , and G\ are functions of q 2 , the invariant mass of the photon. 
Some conventions |79|,^0] denote G\ as — F 3 . 

The F \, F 2 , and G\ form factors are related to the deuteron charge, magnetic, and 
quadrapole form factors, denoted by F(j , F^f , and Fq respectively, by p5t|87| .p0| 


F c = F 1 + -t) [Fi + (1 + rj)F 2 - Gi ], 
Fm = G\ , 

Fq = F 1 + (1 + 7])F 2 - G u 


where 


V = 


-T 

4Mj' 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


Note that since both the initial and final electron and deuteron are on-shell, — q 2 > 0. 
We define Q 2 = — q 2 . 

At zero-momentum transfer, the deuteron charge, magnetic, and quadrapole form 
factors are simply 


F c ( 0) = 1, 

(4.12) 

F m { 0) = —g d = 1.71293, 

(4.13) 

m p 


F q { 0) = M 2 d Q d = 25.8525, 

(4.14) 


where m p is the mass of the proton, /^ is the magnetic moment of the deuteron 
(in nuclear magnetons and Qd is the electric quadrapole moment of the 

deuteron (in units of the deuteron mass, although typically it is measured in e-barns). 
The numerical values of the form factors are determined by using experimentally 
measured values |pS8|.[QT]. fT2 . 
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Another set of form factors are Go, Gi, and G 2 , commonly used in light-front 
dynamics [[79],|80],|3|] . They are simply related to the charge, magnetic, and quadrapole 
form factors, 


Go = F c , 
Gi = -Fi 


Mi 


r , Vs 

G2 = ~y 7 i f q- 


(4.15) 

(4.16) 

(4.17) 


The electron-deuteron cross-section is measured to determine the deuteron form 
factors for large momentum transfers. The unpolarized cross-section has the form 

mu 


da 

dVL e 


CMott 


1 ' M d blIi 2 


A(Q' 2 ) + B(Q 2 ) tail 2 ( y 


(4.18) 


where 


O'Mott — 


a cos -f 
2E e sin 2 % 


(4.19) 


and E e is the electron beam energy, 9 e is the angle by which the electron scatters, and 
a is the fine-structure constant. The structure functions A and B can be expressed 
in terms of the charge, magnetic, and quadrapole form factors, 


A = F 2 + -r? 2 F 2 + - ?/ F 2 ,, 
B = ^ 77(1 + v) f m- 


(4.20) 

(4.21) 


The three form factors iq, F 2 , and G\ cannot be determined uniquely from the 
unpolarized scattering data since that measurement of provides only two structure 
functions. More information is needed, which can be obtained from experiments 

, [S9j . The most 


where the polarization of the electron and/or deuteron is measured 


commonly measured quantity is the tensor polarization observable, T 20 , 
; vFcFq + I VFq + Wl, (1 + 2(1 + rf) tan 2 e f) 


T20 — — 


V~2 (A + B tan 2 y) 


(4.22) 























108 


Note that To o depends on the angle 9 e . For ease of comparison, data for To 0 is usually 
presented for 9 e = 70°. Alternatively, one can eliminate the angular dependence by 
defining T 20 , which is T 20 with F M set equal to zero, 


IvFcFq + I v 2 F* 
V2 {Fq + ^V 2 Fq) ' 


(4.23) 


The extraction of the structure functions from the data is difficult. In practice, 
each cross-section measurement is performed at a different momentum transfer and 
angle, making it impossible to exactly disambiguate, for example, A and B. One 
needs to interpolate between values of B to calculate values for A, and vice versa. 
The Jefferson Lab t 2 o Collaboration used a self-consistent method for obtaining the 


structure functions from the scattering data |94| . They consider a variety of theoretical 
models for the deuteron form factors, then fit the parameters of the models to the 
measured cross-sections. After each model is fit, it is used as the interpolation function 
to disentangle the structure functions. This two-step process helps minimize model- 
dependent effects in the values of A and B. 


4-1.2 Light Front Calculation 

Light-front dynamics is particularly well suited to calculating form factors. One reason 
is that the generators of boosts in the one, two, and plus directions are kinematic, so 
that wave functions calculated with a truncated potential transform correctly under 
boosts. This feature is especially important for form factors at high momentum 
transfer since the wave functions must undergo a large boost. 

Another, more subtle, reason for using the light front is that many of the graphs 
which contribute to the current vanish identically. For example, the three lowest-order 
graphs for the current are shown in Fig. |4.2| . The double line denotes the deuteron, 
and the vertex of the deuteron lines and the nucleon lines represents the deuteron wave 


function. The graph labeled (a) does not vanish and is calculated in section 4.1.3 . 
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Figure 4.2: The lowest-order graphs which contribute to the deuteron current matrix 
element. 


Fig. fh2](b) vanishes in light-front dynamics. To see why, we first note three facts: 
the pins component of each particle in light-front dynamics is non-negative (for mas¬ 
sive particles it must be positive), the plus component of the momentum is conserved, 
and the pins momentum of the vacuum is zero. Combining these facts, we find that 
any vertex which has particles on one side and vacuum on the other must vanish. In 
other words, the vacuum is trivially empty, and no graphs couple to it. 

For Fig. [Ol(c), the coupling of the photon to the nucleon goes like where 

the light-front spinors are defined in Appendix |F]. For /i = +, this matrix element is 


suppressed maximally, and thus J + is the “good” component of the current [p0j [8^, |95 
We calculate only J + , since it is the most stable. 

We do not consider the contribution of higher-order graphs to the deuteron current, 
such as graphs where the photon couples to a meson or a nucleon while a meson is 
present. The omission of the meson-exchange currents does not affect the rotational 
properties of the current, although it does affect the overall values of the deuteron 
form factors. This is acceptable since we are only concerned with the rotational 
properties in this work, not in the detailed results. 

The neglect of the graphs where the photon couples to a nucleon while a meson 
is present may affect the rotational properties of the current. This is because the 
deuteron current shown in Fig. |4.2| (a) is not formally conserved [93], but current 
conservation is necessary (although not sufficient) for the current to have the correct 
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properties under rotation. To construct the conserved current operator associated 
with a given wave function, the current must include diagrams that are related to 
the potential used to calculate the wave function. We expect that these diagrams 
are small since they contain meson propagators, and that neglecting them should not 
significantly affect the conservation of the deuteron current or the rotational properties 
of the current. 


Symmetries of the Electromagnetic Current 


Now, we use symmetries to relate the components of (k',m'\J + (q)\k,m). Although in 
light-front dynamics some generators of Lorentz transformations are dynamic, such as 
the generators of rotations from the x-y plane into the x direction are dynamic, we are 
free to boost in the plane perpendicular to the z-axis, boost in the plus direction, and 
rotate about the z-axis, since the generators of those transformations are kinematic. 
In addition, we will find how the states and the current operator transforms under 
parity and time-reversal. 

The kinematic generators allow us to choose which frame to evaluate the current 
in. We choose the Breit frame ]8(J, where q + = q~ = q v L = 0 and qf = Q. We 
also choose the plus momentum of the deuteron to be M^, since is value of the 
plus momentum in deuteron’s rest frame. To simplify notation, we define the matrix 
elements of J + as 


,(<?) = 



m 


J + (Q) 



(4.24) 


This quantity is represented with / m , m instead of J m / m because J m 'm is used to rep¬ 
resent the matrix elements of J using the instant-form spin basis. This distinction is 
discussed later in this section. 

Next, we consider rotation about the z-axis by n (R z (n)), parity (II), and time- 
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reversal (0). First note that the current operator transforms as 


iUMJ^vr) = J ± , 

(4.25) 

n^n = j t , 

(4.26) 

0J±0 = JT. 

(4.27) 


Obviously, we must use both parity and time reversal to retain the plus component 
of the current. Now look at how the deuteron states transform: 


Rz( tt)| - k_L,m) = | + fcj_,m)(-l)”\ (4.28) 

en| -k ±1 m) = I - l) m . (4.29) 


Although the symmetry operators can introduce extra, constant phase factors, they 
are omitted here since they cancel when calculating matrix elements. These results 
make it easy to find what each of these symmetry operators does to the matrix 
element. Under rotations a rotation of i r about the z-axis, the current matrix element 
transforms to 


- (-1)’”-”' (-yH J+ W) | + T’ m ) ' 


Under parity followed by time reversal, the matrix element becomes 

I&.JQ) - -m\ J + {Q) |+^, -m') . 


(4.30) 


(4.31) 

(4.32) 


Since is Hermitian, we can also take the complex conjugate of the matrix element 
to get 


C.JQ) 


Q± 

2 ’ 


m 


J + (Q) |+y,m'). 


(4.33) 


The appropriate combinations of Eqs. (|4.30| - |4.33|) give JFS 


C.J.Q) = (-i)”A -JQ), 
= (-1 


(4.34) 

(4.35) 
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The same relations also apply for matrix elements. Eqs. ( 4.34 ) and ( |4.35 ) imply 
that of the nine possible matrix elements of J + , there are only four independent 
components. We choose those components to be /^, J^, If_ 1 , and Iq 0 . It is helpful 
to express the matrix elements in a matrix to see the symmetry properties explicitly. 

\ 

(4.36) 



t 

it 

it 

It -1 

T+ = 

m! ,m 


-A + o 

-^00 

Ito 


V 

iU 

-It 

It 


Rotational Invariance and the Angular Condition 

This is as far as we can go with light-front dynamics, but there should be an additional 
redundancy in our matrix elements. We have derived four independent components, 
whereas in a fully covariant framework there are only three form factors. The resolu¬ 
tion of this conflict is that full rotational invariance imposes an angular condition on 
the light-front matrix elements. 

The angular condition can be found by using the deuteron polarization vectors 
for the Breit frame to calculate the current given in Eq. ( f4.3| ) in terms of F \, F 2 , and 
G\. By comparing that result with Eq. (|4.36| ), we obtain linear relations between the 
current matrix elements (//], J^, I ^_ 1 , and Iq 0 ) and the form factors (fq, F 2 , and Gi) 
and also find that in general there is a deviation from the angular condition, which 
we denote with A, given by |83| 


A — + (1 + 2//)/^ + I^~_ 1 2\j2gliQ, 


(4.37) 


where 7] is defined by Eq. ( 4.11 ). Since A vanishes when the deuteron current trans¬ 
forms correctly under rotations, we interpret A as a measure of the extent to which 
the current matrix elements transform incorrectly. 

The form factors are overdetermined by the current matrix elements, which means 
there are many different ways to express the form factors in terms of the current 
matrix elements. In order to make a definite prescription, one the current matrix 
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elements (or a linear combination of matrix elements) is identified as “bad”, and the 
angular condition is used to eliminate it from the expressions for the form factors. 
This classification of matrix elements as “good” or “bad” is similar to the one made 
for choosing which component of the current to use. 

When A is zero, all the prescriptions must be equivalent, while a non-zero A means 
that the form factors depend on which prescription is chosen. Since A is non-zero in 
general, it is important to choose the best prescription to obtain the form factor. 

We consider four different prescriptions |83|] in this work. This allows us to study 
how sensitively the form factors depend on the prescription used. First, Grach and 
Kondratyuk (GK) derive a relation in which /^ is considered bad. Using the angular 


condition to eliminate the for the equations, they obtain [96 


G, 


0 ,GK ~ 


Gi,gk = 21 / 


(3 277 )/^ + 2\[2 7 //^ + Ii_i 

1 „ 


11 


\/2rj 10 

' / 10 


2 ^ ( ,+ 27 / + 

G 2 ,gk - —( * 77/11 + —~ /£—1 


(4.38) 

(4.39) 

(4.40) 


Brodsky and Hiller (BH) use a prescription where is bad, which results in [|97 


Gq } BH 

Gi^bh 

G 2 ,bh 


3(1 + 27/) 
2 


I + 277 


(3 - 277 )/+ + 8 V^Z+g + 2(277 - l)A + -i 


‘00 


/+1 + (277 - 1 ) 


‘10 




2V2 


3(1 + 277 ) 


V^nlto - vloo - (! + 


(4.41) 

(4.42) 

(4.43) 


Frankfurt, Frederico, and Strikman (FFS) start by analyzing the deuteron current 
using an equal-time spin basis, where the current matrix elements are J+, 


Their prescription uses the Cartesian basis to find that J+ is the bad current. Upon 
transformation to the spherical basis and use of the Melosh transformation [82], 98| to 
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relate J+, m to they find 


Go,ffs — 
Gi } ffs = 


3(1 + 7?) 

1 


(2rj + 3)/^ + 2^2rjl± 0 — r) Iq 0 + (2 rj + l)/^] 


G 


1 + rj 

V2 


2,FFS 


3(1 + 77) L 


/++/+_/+ _ 2 ( x 7 ?) /+ 

+J °° J >- 1 ^ ho 

—Tjlti + 2 \Z 2 t// 1 I o — 77 /^ — (77 + 2 )/ 1 t l 1 


(4.44) 

(4.45) 

(4.46) 


Chung, Polyzou, Coester, and Keister (CCKP) choose the canonical expressions for 
the form factors Gq, G 1 , and G 2 in terms of the equal-time current as the starting 


point |£9|,Q. They use rotations and the Melosh transformation to express the equal¬ 
time current matrix elements in terms of the light-front current matrix elements, and 
find that 


0 ,CCKP 


6(1 + 77 ) L 
GifiCKP = Gi^ffs 


(3 - 277 )(/+ + /+) + 10y/2jlt o + (477 - 1)7+ ! 


G 


2, CCKP 


= G 


2 ,FFS- 


(4.47) 

(4.48) 

(4.49) 


Note that all the prescriptions reviewed here are related by the angular condition. 
When A is zero, each prescription gives the same results. 


4-1.3 Impulse Approximation on the Light Front 


We now are ready to relate the deuteron wave function to the current expressed in 


Eq. (|4.24|) . The deuteron wave functions solved in Chapter || are used. We use the 
representation of the wave function in the spin |p ± ,p + , 777.1,777.2) basis, as discussed 
in Sec. |3.5.4|. Note that these spins are expressed in the usual Bjorken and Drcll 


representation [99]. By inserting a two complete sets of states into Eq. (4.24) and 
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making the momentum of particle 1 explicit, we get 




J d 2 p ± dp + \P + ’P-L + 

J (5)(9±) 


m 2 


m 1 ,mi ,7722 


2 (p + ,p ± + y/m', 

P ,Pj. - —, mi, m 2 


+ <7± 

P ,P_L - -y,™i 


Q± 

2 ,m 


(4.50) 


where the spin directions of particles 1 and 2 are labeled mi and m 2 , respectively. Also, 
since the deuteron is an isoscalar combination of nucleons, the isovector component of 
the nucleon current does not contribute and the isoscalar nucleon current is the same 
for both nucleons. This allows us to simply double the isoscalar current of particle 1 
instead of using the isoscalar currents of both particle 1 and 2. 

We want to boost the deuteron wave functions to the rest frame, since that is where 
the wave functions are calculated. Note that the boosts in the x-y plane transforms a 
general vector (k±, k + ) by k± —> k± — leaving k + unchanged. This means that 

the boost that puts the deuteron in its rest frame transforms the individual nucleon 
momentum by 


Q± 


Pi 


d - nf , 


(4.51) 


where we use the Bjorken x -variable x — Note that the spin labels do not change. 
This is because the boost is in the perpendicular direction, so the spin labels (defined 
to point in the z direction) are not affected. 

This allows Eq. (|4.50|) to be written as |k| 




jcPpidp* £ (m'|p+,p ± + (l-*)^,< 


m 2 


^ ,7711,7712 

2 {p + ,P± + ^,m , i 


l) 


P + ,P± - (1 - i)y,mi,m 2 


P + ,P± ~ 


m 


(4.52) 


We have dropped explicit mention of the deuteron’s momentum. 
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Before we continue, we note that we can write the nucleon current matrix elements 
using the u spinors, 


P + ,P± + -7p m i 


^(S)(9±) 


+ 

p ,p± - -y,"u 


= U BD (p + ,p ± + J[s)(*l±) u bd (p + ,p± - I ( 4 - 53 ) 

where we have used the “BD” subscript to denote these as Bjorken and Drell spinors. 
Since the nucleons are on-shell, as they must be in a Hamiltonian theory, we interpret 
«bd (p + ,p ± , m) as wbd (px = P±,x,P y = P±, y ,Pz = where p~ = M . 

However, the matrix elements can be calculated easily if we convert the Bjorken 
and Drell spinors to light-front spinors. This transformation is formally accomplished 
with the Mclosh transformation Rm, given by |79|,[80|], 


Rm — 


p + + M — icr ■ (to x p ± ) 

v/(p+ + M) 2 +p[ ’ 


(4.54) 


where to points in the direction of the light front. Note that the Melosh transformation 
is unitary. 

This transformation converts a Bjorken and Drell spinor to a light-front spinor in 
the following manner, 


ULF<yP + ,P±,rn) = ^2{RM)m,m'U BD (p + ,p ± ,m'). 

m' 

Using this transformation, we rewrite Eq. (|4.52|) as 

) = 


(4.55) 


f d 2 p±dp + ^ 

S( m ' 

P + ,P± + (1 - x)^y, rn”, m 2 ^ (^m)to",to' 1 

/ / 

77lj ,mi ,7712 

m" 



2«lf ( p + , P ± + J ( S )( Q ±) U -lf ( p + , P ±_ ~ ^, m i 

i ( P + ,P± - (1 - x) < Y,rn'",m 2 m^ 


(4.56) 
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4-1-4 The Nucleon Form Factors 

From Lorentz covariance, current conservation, parity invariance, and time reversal 
|86],^7]], the isoscalar part of the nucleon current can be expressed as 


invariance 


•/&,(«) = 7 "ff’(«) + ’(<!). 


Taking the matrix elements of J + with the light-front spinors gives 
u(k', \')J+(q ± )u(k, A) = Fl s> (q)u LF (k', A')7 + u LF (k, A) 


+ F^\c,)^u LF (k', A')7 + 7l u FF (k, A), 


(4.57) 


(4.58) 


where we have parameterized 


k' + = k + , 

^-L = P± T i 

k ’ L 

*±-P± 2 


(4.59) 

(4.60) 

(4.61) 


We use the representation of the spinors given in Appendix |] to simplify the spinor 
matrix elements appearing in Eq. |4.58| . First consider 


u LF (k', X'h + u LF (k, A) = —xi F ,y (M - ^ ■ k r± ) A + + k + A_ 


k+ S 

~ m 8x ' x - 


. _L /-/-La 


7" 


x [A_(M + c*- 1 ■ k L ) + A+fc + ] xlf,a, 


(4.62) 

(4.63) 


Next, 


u LF (k', X')'y + 'y ± u LF (k, A) = j^Xi F ,y ( M ~ a± ■ k ' )A+ + A; + A_ 


. _L ,/-La 


7 7 




x [A_(M + a x ■ fc 1 -) + A + k + ] xlf.a, 


A 


Alf,a'7 7_lAlf,a 5 
£;+ 

M' 


M 

= a 


(4.64) 

(4.65) 


(4.66) 
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Combining these, we find 


k + 




(S)/ 


2 M 


(k',\'\J+ s) (q ± )\k,\) = — X \, 

We can rewrite Eq. (|4.67| ) as 

(k',\'\J+ ) (q ± )\k,\) = + 

d s) to(OV ( t s) (9x)M), 


Xa 


(4.67) 


(4.68) 


which, when inserted into Eq. (f4.56|) , gives 


(4.69) 


where we define 


= / d 2 P±dp + ( m '\P + iP± + (! - x)^,m' 1 ,m 2 ) 


m\, mi, m2 


2(p + ,p ± + 

(p + ,P± - (1 


(4.70) 


for i = 1, 2. Note that both J^ m , m (q ) and , m (g) must satisfy the same equations 
as Jj), m (g) does. In particular, this means that the angular condition should apply 
to m (q) and ,/W m , m (q) independently, so we consider the deviation from the 

angular condition for each. 

There are many parameterizations of the isoscalar nucleon form factors F[ s \q ) 
and F 2 (q). Since the measurement of the electron-nucleon cross section is difficult, 
the data have large errors and are consistent with several different models of the 
nucleon form factors. Some of the models representative of those proposed in the 
literature are: the dipole model, fit 8.2 of Holder [j l00|| , Gari:1985 [|101|| , model 3 of 
Gari:1992 ||102| , |103f] , best fit for the multiplicative parameterization of Merged ||104 


and model DR-GK(l) of Lomon [|105|| . The f\ S) ( q) and F!f s (q) form factors for each 


<s). 


of these models are shown in Fig. 4.3. 
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Nucleon Form Factors for several models 



Q 2 [(GeV/c) 2 ] 


Figure 4.3: The F\ and — F 2 isoscalar nucleon form factors for six different models: 
the dipole model, Holder ||100|| , Gari:1985 ||101|| , Gari:1992 ||102| , |103|| , Merged ||104|| , 


and Lomon [105 
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We can relate the isovector form factors to Ge p , Gm p , Geu, and Gmu, the pro¬ 
ton electric, proton magnetic, neutron electric, and neutron magnetic form factors, 
respectively, with [ |10G| 


jp(S) _ Ge p + Geti + T ( Gm p + GMn) 


F 2 (S) = 


-G 


E P 


2(1 + r) 

Geu + Gm p 


G 


Mn 


2(1+ t) 


(4.71) 

(4.72) 


where r = The value of r is approximately 1 at a momentum transfer of about 
5 GeV 2 , the upper range of momentum transfers that we consider. Since the overall 
magnitudes of the form factors are similar at this momentum transfer, it is important 
to measure each of the form factors with the same accuracy and cover the same range 
of momentum transfers. Currently, the most poorly known form factor is Gehi both 


in terms of the magnitude of the error and in the number of data points ||105 . 

In Section fO, we will find that for momentum transfers greater than about 


2 GeV 2 , the spread in the values of the deuteron form factors due to the break¬ 
ing of rotational invariance on the light front is smaller than the spread in values 
due to using the various nucleon form factors. It is uncertainty of the nucleon form 
factors, not the use of the light front, that limits the accuracy of the deuteron form 
factors at large momentum transfers. Only more accurate measurements of the nu¬ 
cleon form factors, especially Geu , will allow for more accurate deuteron form factor 
calculations. 


4.2 Axial Form Factors 

The formalism used for the axial current and form factor is very similar to that used 
for the electromagnetic current and form factor. Thus, most of the discussion from 
the previous section carries over here with only slight modifications. We highlight 
only the differences. 

The derivation of the symmetries of the axial current matrix elements is almost 
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the same as in Section |4.1.2| , with the exception that under parity, the axial current 
picks up a negative sign 

n = - J\ f. (4.73) 

When this is propagated through the algebra, we find that 

(4.74) 

= (-l)”-”*'/+ ram ,(Q). (4.75) 

Eqs. ( fl.74|) and (|4.75|) imply that of the nine possible matrix elements of J^, there 
are only two independent components. We choose those components to be /+ n , and 
/ ( + )10 . It is helpful to express the matrix elements in a matrix to see the symmetry 
properties explicitly. 

( 


/+ 

(5)772', m 


V 


-+ 

7 + 

0 ^ 


(5)11 

+5)10 


7 + 

+5)10 

0 

I 

(oi + 

0 

(4.76) 

0 

7 + 

+5)10 

J (5)H / 



We have derived two independent components, but an analysis of the covariant 
theory shows that only one deuteron form factor (Fa) contributes for the plus com¬ 
ponent of the axial current | 8 l|]. This implies that the requirement of full rotational 


invariance imposes an angular condition on the light-front axial current matrix ele¬ 
ments. The deviation from the angular condition, denoted by A, given by [ |3T| 1, 

s/Zj, 


A — v 7 + — /+ 

^ — 2 ( 5 ) n (5)10* 


(4.77) 


Since the deuteron axial form factor is overdetermined by the current matrix 
elements, we need to classify the current matrix elements as either “good” or “bad” 
to eliminate ambiguity. We consider two such choices. 

Frankfurt, Frederico, and Strikman (FFS) find that the J + 2 is the bad matrix 
element [^]. After transforming to the spherical basis and using the Melosh trans¬ 
formation, they find that 


Fa = 


2(1 + 77) ' (5)U 


^( 5)11 + V^V 1 * 


(5)10 


(4.78) 
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Frederico, Henley, and Miller (FHM) use the behavior of the matrix elements in the 
non-relativistic limit to determine that the bad element is 77 


(5)10 


Fa 2 ^(5)ii • 


81]]. This means that 
(4.79) 


The current matrix elements are calculated using the nucleon axial current. The 
general form of nucleon axial current is given by ||81|| , 


y 5) „ = 7'V*1 + q^FZ. (4.80) 

Since we choose to /i = + and work in the Breit frame, where q + = 0, Eq. (|4.8Q|) 
reduces to 


y 5 )n = 7 + 7 5 ^. 


(4.81) 


The light-front spinor matrix elements of Eq. ( |4.81| ) can be computed using expres¬ 
sions given in Appendix |F|. The result is 


( k ', X '\ J (5)n( ( l±) IM) = T7 F a((1)X\'^X\- 


M 


(4.82) 


We consider only one model for the nucleon axial form factor since the deuteron 
axial current has such a simple dependence on it. We choose to use the dipole model 

Fa( 0) 


Fa{Q 2 ) = 


1 + 


Q 2 


2 ’ 


(4.83) 


where Ma is the axial mass. For our calculations, we use the value for the axial mass 
determined by Liesenfcld et al. ||107| . 


4.3 Results for the Form Factors 

We use the deuteron wave functions obtained for the light-front nucleon-nucleon po¬ 
tential in Chapter ^ to calculate the deuteron currents and form factors. This gives a 
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solution where light-front dynamics is used consistently throughout. For the poten¬ 
tial, we choose the light-front nucleon-nucleon potential with f a = 1.2815. We have 
verified that the results do not change significantly when f a = 1.22 is used. 

Figure |4.4| show the currents and the associated angular condition for /^, given 
by Eq. ( fl. 70 | ), for several different deuteron wave functions. Results are shown for the 
wave function from the OME, OME+TME, and OME+TPE potentials (calculated 
in Chapter ^j), and the parameterization of the deuteron wave function for the energy 
independent Bonn potential |58|]. The current matrix elements (but not A) are ap¬ 
proximately the same regardless of which wave function is used. This consistency is 
important, since it verifies that the gross features of all the models are the same. 

We find that A for is much smaller than the largest matrix elements when 
using the OME wave function. This means that the matrix elements of the current 
operator transform very well under rotations. This is somewhat surprising, since we 
found earlier that the binding energies for the OME wave functions have a large 
splitting, indicating that OME wave functions transforms poorly under rotations. 

To understand this apparent discrepancy, note that some diagrams are not in¬ 
cluded in the calculation of the current although they are required for manifest current 
conservation. This problem is addressed in Section [4.1.2| . The missing graphs should 
have a small effect on the rotational properties of the current, but they may have a 
large effect on A. In particular, since the OME+TME current is missing more graphs 
than the OME current, it is reasonable that the OME+TME current transforms more 
poorly under rotations than the OME current. 

Comparing the current calculated with the OME wave function to those calculated 
with other potentials, we find that for momentum transfers of more than 1 GeV 2 the 
OME /A current has the best transformation properties under rotation of all the 
currents shown. 

For smaller momentum transfers, the transformation properties of the Bonn and 
OME+TME wave functions are the best. This is expected, since in the limit of 
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no momentum transfer, the current is simply the overlap of deuteron wave 

functions, (m'\m). If the initial and final states have the same mass, the matrix 
element is simply which satisfies the angular condition. However, if the states do 
not have the same mass (which implies that m! ^ m), there will be a non-zero overlap 
between the two states, which violates the angular condition. Since the masses of the 
deuteron states are exactly the same for the Bonn wave function, and approximately 
the same for the OME+TPE wave function, they have a small A at low momentum 
transfer. However, the OME wave functions, having the largest mass splitting, have 
the largest A at low momentum transfers. 

Figures fh5| and |4.6| show the current matrix elements and the angular condition 
for /A and /^, respectively. The general features of these figures are the same as 


in Figure |4H], with one important exception. In both figures, the A for the OME 
wave function has about the same magnitude as the A’s for the other wave functions. 
This means that the rotational properties of /A and /A currents are approximately 
the same regardless of which wave function is used. This is a surprising result since 
it indicates that the rotational properties of the current depend more on how the 
current is constructed than on which wave function is used. 

In Fig. |4.5| , we find that the magnitude of A is almost the same as the magnitude 
of the largest matrix element of /^. This means there is a large deviation from the 
angular condition, and that form factors calculated with this current may depend 
strongly on which matrix element is chosen as “bad”. We show below that this is not 
the case for the electromagnetic form factors. 

We find that A is much smaller than the largest matrix element of the axial 
currents shown in Fig. |4.6| for most values of momentum transfer. This means that 
the deuteron axial form factor will be essentially independent of which matrix element 
is chosen as “bad”, except for within the range of 1.5 to 2 GeV 2 . 

Now we combine the two parts of the electromagnetic current, /A and with 


the nucleon form factors F 1 and F 2 to get the total current. Figure T7 shows the 
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l + (1) Matrix Elements 



012345012345 



012345012345 
Q 2 [(GeV/c) 2 ] Q 2 [(GeV/c) 2 ] 


Figure 4.4: The matrix elements of I(\y m > rn , the component of the electromagnetic 
current which multiplies the nucleon F\ form factor, calculated with the wave function 
from the OME, OME+TME, OME+TPE, and Bonn potentials. 
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l + (2) Matrix Elements 



012345012345 



Q 2 [(GeV/c) 2 ] Q 2 [(GeV/c) 2 ] 

Figure 4.5: The matrix elements of /T m / m , the component of the electromagnetic 
current which multiplies the nucleon F 2 form factor, calculated with the wave function 
from the OME, OME+TME, OME+TPE, and Bonn potentials. 
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I* Matrix Elements 

(5) 




Figure 4.6: The matrix elements of the deuteron axial current including 

the nucleon axial form factor, calculated with the wave function from the OME, 
OME+TME, OME+TPE, and Bonn potentials. 
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currents for F\ and F 2 /^, as well as the sum, J + . The Gari:1985 nucleon form 


1 ( 2 ) 


factors are used [|K31|| - We hnd that F\ 1^ gives the largest contribution to the total 
current, and because A is small for /A, A is also small for the total current, meaning 
that the total current transforms well under rotations. Thus, in spite of the fact 
that A is approximately the same size as the current matrix elements for /^, the 
deuteron electromagnetic form factors should not depend too strongly on the choice 
of the “bad” matrix element. This is especially true for the form factors calculated 
with the OME wave function. 

We calculate the form factors A, B , T 2 o, and Fa using the OME wave function, 
and show the results in Fig. |4.8| . In general, the form factors do not depend strongly 
on which matrix element is chosen as “bad”, in agreement what what we predicted in 
the previous paragraph. The only exception is for the B form factor, and to a lesser 
extent the Fa form factor, near where they crosses zero. This is not too surprising, 
since a small constant shift in any function near a zero-crossing has a large effect in a 
logarithmic plot. Also, we note that the FFS and CCKP choices of the “bad” matrix 
element give the same value for B. This is a consequence of Eqs. ( ff.48| ) and (|4.49| ). 

We also use the OME+TME wave function to calculate the form factors A, B , T 2 o, 
and Fa, which we show in Fig. fOj. We argued earlier that the these electromagnetic 
form factors depend more strongly on which matrix element is chosen as “bad” that 
those calculated with the OME wave function, and that dependence is clear in this 
figure. At low momentum transfers, the dependence on the change is fairly small, but 
as the momentum transfer increases, so does the dependence. The axial form factor 
is not affected as strongly, primarily because each wave function generates an axial 
current which violates the angular condition by approximately the same amount. 

Since there are many different models of the nucleon electromagnetic form fac¬ 
tors, we calculate the deuteron electromagnetic form factors using each of them to see 
what effect the differences have. The results are shown in Fig. [4.10[ At low momen¬ 
tum transfers, all the nucleon form factors give close to the same results. However, 
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Figure 4.7: The matrix elements for F/jT , m , and J+, m calculated with 

the OME wave functions. The Gari:1985 nucleon isoscalar form factors are used for 


F] and F 2 


101 
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Figure 4.8: The form factors A, B , T 2 o, and Fa calculated using the various choices 
of the “bad” matrix element. The definitions of the “bad” matrix elements are given 
in Sections |4.1.2| and [4 . 2[. The OME wave function is used, along with the Lornon 


nucleon form factors ||105|1 for the electromagnetic form factors, and the Liesenfeld 


nucleon form factor [107] for the axial form factor. 
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Figure 4.9: The form factors A, B , T 2 o, and Fa calculated using the various choices of 
the “bad” matrix element. The definitions of the “bad” matrix elements are given in 
Sections |4.1.2| and |4. 2|. The OME+TME wave function is used, along with the Lomon 


nucleon form factors ||105|1 for the electromagnetic form factors, and the Liesenfeld 


nucleon form factor 1107 for the axial form factor. 
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when the momentum transfers is large, we find a large spread in the values due to 
nucleon form factors. In fact, this spread in A and B is larger than the spread of 
values obtained from using different “bad” matrix elements with the OME+TME 
wave functions. In other words, in order to obtain accurate results for A and B at 
momentum transfers over 2 GeV 2 , it is more important to determine which nucleon 
form factor to use than when “bad” matrix to use. The spread in T 2 o is about the 
same size regardless of whether the different “bad” matrix elements or the different 
nucleon form factors are used. 


Finally, in Fig. [4. 1 1|, we compare the A, B, To o, and Fa form factors for the OME 
and OME+TME wave functions to experimental data. The “bad” component was 
chosen according to FFS, and the nucleon form factors of Lomon were used for A, 
B , and T 20 , while the Liesenfeld axial nucleon form factor was used for Fa- The 


data for A is from: Buchanan et al. ||108||, Elias et al. [|109||, Galster et al. ||1 10 


Platchkov et al. [ 111 ], Abbott et al. [112], and Alexa et al. [l]; the data for B is from: 


Buchanan et al. ||108|| , Auffret et al. ||1 13|1 , and Bosted et al. ||114j| ; and the data for 
T 20 is from: Schulze et al. ||115|| , Gilman et al. [|116|| , Boden et al. ||1 1 7|| , Garcon et 
al. |118|] , Ferro-Luzzi et al. [|119|1 , Bouwhuis et al. [|1 2 0| , and Abbott et al. [|121| . 


There is a rather large difference between the form factors calculated with the 
OME and OME+TME wave functions. This difference is due primarily to the fact 
that the OME wave functions are more deeply bound than the OME+TME wave 
functions, and it can be reduced by choosing a different sigma coupling constant f a 
for the OME and OME+TME potentials. However, for our analysis of rotational 
invariance, it is important to keep f a fixed. 

The difference between the calculated form factors and the data is also quite large. 
This is not unexpected, since in our model of the current, meson exchange currents 
are not included. It is known that these have a large effect on the form factors at large 
momentum transfers @00. Including these effects would bring the form factors 
into better agreement with the data. However, we emphasize again that agreement 























































































133 



N 

^ 1 ^ 

B 

1 1 1 1 1 1 

\ 


-Dipole 


\ 

-Hohler 

- 

-Gari:1985 - 


\ 

- - Gari:1992 

- 

\ 

s-® Mergell 

h -□ Lomon 


\ 

^ - - 


V 


- 

1 1 




i 



Q 2 [(GeV/c) 2 ] 


Figure 4.10: The electromagnetic form factors A, B , and T 2 q calculated using the 
various choices of the nucleon isoscalar form factors. The OME wave function is 
used, along with the FFS choice of the “bad” deuteron current matrix. The axial 
form factor is not shown since its dependence on different form factors is trivial. 






134 


with the data is not a priority of this work. Our goal is to gain a better understanding 
of the breaking of rotational invariance by the light front, and how to restore that 
invariance. Only after we have that understanding can we pursue accurate calculation 
of the form factors with light-front dynamics. 
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Figure 4.11: The A, B , T 20 , and Fa form factors for the OME and OME+TME 
potentials, along with data. See the accompanying text for an explanation of the 
data. 
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Chapter 5 

CONCLUSIONS 

“If a thing can not go on forever, it will come to an end.” 

—Herbert Stein 

The issue of rotational invariance in light-front dynamics must be addressed be¬ 
fore one attempt to use light-front dynamics for high-precision calculations. In this 
work, we have sought to find ways to quantify the level to which rotational invariance 
is broken. In addition, we have investigated potentials of different orders, since we 
know that if all orders of the potential were included, rotational invariance would be 
restored. We have used light-front dynamics to obtain one-meson-exchange (OME) 
and two-meson-exchange (TME) potentials for several different nuclear theory La- 
grangians. 

In Chapter |2], the massive Wick-Cutkosky model is used to investigate the break¬ 
ing of rotational invariance in deeply bound states. We studied the pseudo angular 
momentum decomposition of each wave function into different eigenstates of pseudo 
angular momentum. In general, we found that most states are approximately angular 
momentum eigenstates. That is, a partial-wave decomposition of the wave functions 
showes that most states have one dominant (highest percentage) angular momentum 
component. Including the TBE potentials increases the fraction of the dominant 
component in most states. We also showed that the spectra generated by the Bethe- 
Salpeter equation using the ladder and crossed ladder are well approximated with the 
spectra generated by the matching OME+TME potentials. In addition, we found 
that the spectra for the full ground state calculation cannot be reproduced accurately 
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for deeply bound states using either the light-front Hamiltonian or the Bethe-Salpeter 
equation. 

In Chapter |3], a model Lagrangian for nuclear physics which includes chiral symme¬ 
try was used to derive two models, a light-front pion-only model, and a new light-front 
nucleon-nucleon potential. Those potentials are then used to calculate the binding 
energy and wave function for the m — 0 and m = 1 states of the deuteron. We find 
that for both models, the splitting between the m = 0 and m — 1 states was smaller 
for the OME+TME potential as compared to the OME potential gives. We also find 
that the effects of chiral symmetry must be included to obtain sensible results for the 
OPE+TPE pion-only model. 

In Chapter [|, the wave functions obtained in Chapter |3] are used to calculate 
the form factors of the deuteron using only light-front dynamics throughout. In 
light-front dynamics, there are four independent components of the deuteron current. 
However, the requirement of rotational invariance introduces an angular condition 
that the four components must satisfy, reducing the number of physically independent 
components to three. The deviation of the calculated current components from the 
angular condition is denoted by A. We find that A is very small for the deuteron 
wave functions calculated with the OME potential. This is an important result, since 
it means that although in principle the light-front calculation of the deuteron current 
does not transform correctly under rotations, in practice it does quite well. The 
smallness of A means that any reasonable prescription for eliminating the dependent 
component of the current gives essentially the same results; the uncertainty introduced 
by the various nucleon form factors is greater for A and B, and about the same size 
for T 20 . 

We also found that A is significantly larger when the TME potentials are used. 
Since the previous results in Chapters [| and indicate that the rotational properties 
of the TME wave function are better that for the OME wave function, we interpret the 
increase in A as an indication that extra diagrams need to be included in the current 
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calculation to restore rotational invariance. That is, having a wave function with 
good rotational properties is not sufficient to obtain matrix elements of the current 
operator with good rotational properties. 

Future Extensions 

“I never see what has been done; 1 only see what remains to be done.” 

—Marie Curie 

This work can be extended in many ways. Perhaps the most important addition 
is to add the higher-order graphs to the denteron current that are required for current 
conservation. This may improve the quality of the deuteron form factors when the 
wave function is calculated using the two-meson-exchange potentials. 

Another obvious extension is to include the crossed TME diagrams. This would 
have two effects. Recall that we included the chiral contact diagrams using a pre¬ 
scription to maintain approximate chiral symmetry. When the crossed diagrams are 
included, we simply include the full contact potentials which included chiral symme¬ 
try without approximation. The second effect is that the crossed diagrams are needed 
in principle to obtain an accurate solution for the deuteron. However, these effects 
are not needed to study the breaking of rotational invariance. 

The contribution to the deuteron form factors from the meson form factors is 
another effect that could be investigated. This is possible since our potential is derived 
from a Lagrangian, and we can calculate the mesonic components of the deuteron wave 
function. Obviously, the contribution to the form factor due to the mesons must be 
taken into account when attempting a high precision calculation. However, including 
these components, like including the crossed diagrams, does not affect the overall 
rotational invariance properties of the deuteron form factors. 
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Appendix A 

NOTATION, CONVENTIONS, AND USEFUL 

RELATIONS 

“The best way to be boring is to leave nothing out.” 

— Voltaire 

For a general four-vector a with components (a 0 , a 1 , a 2 , a 3 ) in the equal-time basis, 
we define the light-front variables 

a ± = a° ± a 3 , (A.l) 

aj_ = (a 1 , a 2 ), (A.2) 

so the 4-vector a 11 can be denoted in the light-front basis as 

a — (a + , a~, a±). (A.3) 

Using this, we find that the scalar product is 

a ■ b = = - (a + b~ + a~b + ) — a± ■ 6j_. (A.4) 

This defines g^ u , with g + _ = g_ + = 1/2, g u = g 22 = — 1, and all other elements 
of g vanish. The elements of g ,llJ are obtained from the condition that g^ u is the 
inverse of g^, so g al3 g/ 3 \ = Its elements are the same as those of g MI/ , except that 
g~ + = g + ~ = 2. Thus, 

a ± = 2 a T . (A.5) 

and the partial derivatives are similarly given by 

d ± = 2cU = 2—^—. (A.6) 

T dxT K J 
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To find the physical consequences of this coordinate system, consider the commu¬ 
tation relations \p^,x y ] = ig llL ', which yield 

[p ± ,x T ]=2z, (A.7) 

\P±, X ±}= (A. 8) 


with the other commutators equal to zero. This means that x\ is canonically conju¬ 
gate to p l ± , and x ± is conjugate to p T . Since x + plays the role of time (the light-front 
time) in light-front dynamics, and p~ is canonically conjugate to x + , this means that 
p~ is the light-front energy and that the light-front Hamiltonian is given by P~. 

In any Hamiltonian theory, particles have the an energy defined by the on-shell 
constraint k 2 = m 2 . This implies that the light-front energy of a particle is 

m 2 + k 2 , 


k~ = 


k+ 


(A.9) 


The independent components of the momentum can be written as a light-front three- 
vector fc L f, denoted by 


&lf = ( k + , k ±). 

For dealing with spin, we require the Pauli sigma matrices, which are 


(A.10) 


(J2 ^. 3 \ 

(T , cr , (7 — 



(A.ll) 


The Bjorken and Drell convention for the gamma matrices is used in this 
work. They specify that 

7 1 0 . 

(A.12) 
(A.13) 

y 5 = i7°7 1 7 2 7 3 = . (A. 14) 
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The spin matrices S l then are 


S i = 


S i = 



(A.15) 
(A.16) 


Using E, we can express the hclicity operator as H = p ■ X, which has eigenvalues 
±1. This is useful since the helicity is invariant under rotations, a property that is 
used in Appendix ^|. 

It is useful to define the spinor projection operators A± by 


A± = ^7 T 7 ± = ^7°7 ± = ± « 3 )- 

These satisfy the requirements for projection operators, 


(A-17) 


A+ + A_ = l, (A.18) 

(A±) 2 = A ± , (A.19) 

A ± A t = 0. (A.20) 


We summarize the effect these projection operators have on the gamma matrices: 


A±7° = 7°A t , (A.21) 

A±7 ± = 0 = 7 ± A t , (A.22) 

A±7 =f = 7^ = 7 T A t , (A.23) 

A±7^ = 7" L A±, (A.24) 


and under conjugation, 


7°A±7° = A t . 


(A.25) 
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Appendix B 

CONVERSION TO MATRIX FORM 


“Let us work without theorizing. It’s the only way to make life en¬ 
durable.” 

— Voltaire 


To solve for the bound-state wave function for the Wick-Cutkosky model numeri¬ 
cally, the light-front Schrodinger equation given in Eq. (|2.49|) must be discretized and 
cast in matrix form. The solution of the light-front nucleon-nucleon model proceeds 
along similar lines. The equation is first symmetrized to get 


dpET 


t/2 


dO P V p m (kET, Ok] Pet, # p )AmApet, 0 p ) — AvA^et, 0 k ), (B.l) 


where 


Vp t m(k etj Ok] Pet, 0 p ) —-BAetMAet, Ok)Vp pn (kET, Ok] Pet, O p )A(p ET , O p )B(p ET ), (B.2) 


ET, Ok) — £> Aet) l A(kET, Ok)i> P ,m{kET, Ok), 
A{kET, Ok) = 


2ik^ k<2 kp'rp sin Ok 


(. k c 


-B(/cet) — 


E 2 - 4(fc°) 


Before discretizing the integrals, note that we may write |37 

f( a fu 


dp f ip) — a \ du(f(au) + 


u 


(B.3) 

(B.4) 

(B.5) 


(B.6) 


Using this trick, the integral in Eq. £3 over pet can be written as an integral over a 


finite range. Since we are concerned with a bound state, the wave function is damped 
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for large momenta, and the second term of Eq. (B. (;) converges as u approaches zero. 
Another choice that can be used is [|5§1 ,|59| 


dpfip ) = 


an 


dx 


'-1 cos 2 






/ ( a tan 


Tl(x + 1 ) 


(B.7) 


After using Eq. 


or Eq. (El), the integrals in Eq. |B.1| are over a finite range 


and can be discretized using Gauss-Legendre quadrature. The specific routines for 
the quadrature are given by Numerical Recipes in C ||124||. This conversion gives a 


matrix equation that approximates the original Eq. |B.l 


-V r %(g(E),E)^ m = ^ m , (B.8) 

where the explicit dependence of V^ m on the binding energy E and the coupling 
constant g is shown. This equation must be solved self-consistently for the spectrum 

9(E). 

The approach we use is to first solve for the spectrum for the OBE potential. The 
eigenvalue equation 


V s 

v p,m, OBE 


( E )^p,m = Ot1pp im , 


where 


a = 


(B.9) 


(B.10) 


90Be(E) 2 

The ground-state wave function is the eigenvector that corresponds to the smallest 
eigenvalue a. We calculate the wave function and the smallest coupling constant 
using EISPACK ||125| , |126|| routines for a range of energies to map out the spectrum. 

Using the coupling constant for the OBE potential as a starting point, we can 
use Eq. [B.8| for higher-order potentials that include N meson exchanges. For a given 


energy, the coupling constant g(E) is initially chosen as goBE(E), then we take the 
equation 


N 


E»( B ) 2 ”W,»be(B) 


n=l 


^p,m = 


(B.ll) 
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where 14be is the n-boson-exchange potential, and solve it as an eigenvalue equation 
for (3. The coupling constant g(E) is varied until the the lowest eigenvalue is (3 — —1, 
at which point g(E) is the correct value of the spectrum corresponding to the ground- 
state wave function 
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Appendix C 

AZIMUTHAL-ANGLE AND LOOP INTEGRALS 


We outline the methods used to perform the azimuthal-angle integrations for OBE, 
TME, and MGF potentials, and the loop integrals for the TME potentials. Although 
these methods are specifically for the Wick-Cutkosky potentials (presented in Chap¬ 
ter ^[), they can easily be generalized for the light-front nucleon-nucleon potentials 
(presented in Chapter ||) by using the integrals discussed in Appendix 0. 


C.l Azimuthal-angle Integration of the OBE and MGF Potentials 


“Never express yourself more clearly than you are able to think.” 

—Niels Bohr 

In this section, we evaluate the azimuthal-angle integration of the OBE potential in 
Eq. (|2.61|) and the first term in MGF potential in Eq. ( 2.123 ), using the prescription 
for azimuthal-angle integration given in Eq. (|2.52|) . One of the integrals is easily 
done since since the potential is independent of the azimuthal angle between the two 
perpendicular momenta, so 

dcj) 


V(k + , k ± -,P ,P±) oc 


r-27r 


9(x - y ) 
% - *) 


l 0 A\ + B cos (j) 
r 2n dcj) 

' 0 A 2 + Bcos(j) 


where 


M = (kf - pt)(E - p 1 - k 2 ) - p 2 - p\ - k 2 ± , 

A 2 = (pi - ki)(E - kf -Pz) - P 2 -pi - k 2 ± , 


(C.l) 

(C.2) 

(C.3) 

(C.4) 


B = 2k^p±. 
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The integrals in Eq. E3 are easily done to give 


r»2?r 


d(j) 


-2vr 


(C.5) 


Jo A + B cos (j) V A 2 - B 2 ’ 

which is negative since A < 0. Using this, the azimuthal-angle-averaged OBE poten¬ 
tial is given by 

2 


f M\ 

V 0 BE(k + ,k ± ]p + J p ± ) = - 2vr l — j E 


9(x-y) 9{y-x ) 


yM? - B 2 \Mi - A 


(C.6) 


It is straightforward to rewrite this equation for the potential in terms of the equal¬ 
time coordinates. 

When the other terms in the MGF potential are azimuthal-angle averaged, inte¬ 
grations similar to the one given in Eq. |C.5| are encountered, with the denominator 
squared or cubed. We note that 


-2 7T 


d(j) 


-2n 


A 


lo {A + B cos0) 2 y/A 2 — B 2 A 2 — B 2 ' 
f 2n dcj) -2vr 2A 2 + B 2 


Jo {A + B cos 0) 3 y/A 2 -B 2 {A 2 - B 2 ) 2 ' 

so the azimuthal-angle-averaged MGF potential is given by 


(C.7) 

(C.8) 


1 / mgf(^1) P\i P) 

n „(M\ 2 \ 9{x — y) ^ , N P)2l + W,i 2 , N Pt21 N k)12 \ 

“ \e) 2 £) 1>2 + 2D 1i3 J 

0(?/ — x) ^ Ac,21 + ^Vp,12 + ^Vp,12^Vfc,21 ^ 

y/ A 2 — B 2 V 2/} 2 ,2 2/12,3 / 

where 


Pi,2 = 
A,3 = 


A 2 -B 2 

I ’ 

_ {A 2 - 5 2 ) 2 
2 A 2 + B 2 1 


(C.9) 


(C.10) 

(C.ll) 


and i — 1,2. 
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C.2 Azimuthal-angle and Loop Integration of the TBE Potentials 

“A child of five would understand this. Send someone to fetch a child 
of five.” 

—Groucho Marx 


As in Section |C.1| , we want the azimuthal-angle integrals of the TBE potentials 
given in Eqs. For these potentials, there is also a loop integral that has 

to be done. We start by analyzing the equations schematically. Each of the terms in 
the TBE potentials can be written in the following form, 


V T BE(k + ,k ± - 1 p + ,p ± ) 

q±dq± f 1 


JO 2(27r,r Jo 

I(k + ,k ± ,q + ,q±,p + ,p _ L ) 

r 2n 


3 , dzJ(k + ,q + J p + )I(k + J k ± ,q + ,q ± pp + ,p ± ), 


dfpqdfpp 


>0 


Ai + Bi cos 4> q 


(C.12) 


x 


X 


A 2 + B 2 COS 4> q + C 2 COS <f)p + D 2 COS (4> p - 4> q ) 

1 

A 3 + B 3 COS (j) q + C 3 COS 4>p + D 3 COS (4>p - (j)q) ’ 


(C.13) 


where the A’s, B’s, C’s, and DA may have dependence on k + , k ±, p + , pj_, q + = zE, 
and q ±; they are independent of the azimuthal angles. These functions can be easily 
determined for each potential by examining the forms of the original equations. The 
rotational invariance of the potential about the z-axis allows the </>£ integration to be 
done trivially. 
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In the integrand of /, only the last two terms depend on (f> p . To emphasize this, 
we write 


I(k + ,k ± ,q + ,q ± ,p + ,p ± ) 


-2 IT 


d(j)q 


I 2 (k + ,k ± ,q,p + ,p ± ) 


l 0 “ q (Ai + Bi COS (j) q )(A 2 + B 2 COS 4> q )(A 3 + B 3 COS(j) q ) ’ 
h (k + ,k ± ,q,p + ,p ± ) 

= f 2 * __ 

J 0 (1 + a 2 cos + b 2 sin (p p ) (1 + a 3 cos <p p + b 3 sin <p p ) ’ 
where, for i = 2, 3, 

Ci + D, cos (j) q 


(C.14) 


(C.15) 


d{ 


bi = 


Ai + Bi COS (p q 1 
Di sin 6 n 


Ai + Bi COS 4> q ' 
The integral in J 2 is evaluated to obtain 


(C.16) 

(C.17) 


l2(k + ,k ±1 q,p + ,p ± ) 


2vr 


X 


(a 2 - a 3 ) 2 + ( b 2 - b 3 ) 2 - ( a 2 b 3 - a 3 b 2 ) 2 
a 2 (a 2 - a 3 ) + b 2 {b 2 - b 3 )_ a 3 (a 3 - a 2 ) + & 3 (& 3 - 6 2 ) 


(C.18) 


y/1 - a| - 6| v/ 1 - °l - b i 

The remaining three-dimensional loop integral in Vtbe on q is done using numeric 
techniques. The trick introduced in Appendix [B] is used to convert the semi-infinite 
q± integration into an integration on a compact range. Before doing the z integral, the 
range of integration is limited by using the step functions. Gauss-Legendre quadra¬ 
ture, given by Numerical Recipes in C ||124|| , is used to evaluate all the integrals. 

Since each of the parts of the full TBE potential (TBE:SB, TBE:SX, ... ) given 
in Chapter [2] should be Hermitian and invariant under interchange of particle 1 and 
2, these invariances can be used as a self-consistency check. Each matrix element is 
calculated twice, first by using the straightforward approach, then particle labels 1 and 
2 are interchanged and it is calculated again. The results are compared, and if they 
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differ by an unacceptable amount, the number of quadrature points is increased and 
the element is recalculated. In order to get the numerical accuracy of the potentials 
correct to within 1%, we start with ten points for the q± integral, six points for the (f> q 
integral, and three points for the z integral, resulting in a three-dimensional integral 
using 180 points. 
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Appendix D 


COSINE INTEGRALS 


“According to my point of view, logic and aesthetics cannot be in 
conflict with one another. Perhaps there is something lacking in my logical 
reasoning.” 

— M.C. Escher 

We wish to solve 



(D.l) 


where m is an integer, rij are a positive integers, and a,j and bj real numbers. Note 
that if one of the bj ’s is zero, then we can reduce the integral to one with N — 1 sets 
of parameters. We therefore assume that bj ^ 0. 

Some things to note about this integral: 

1. Since the range of integration is 27T, and all the functions in the integrand are 
periodic in that range, the range of integration can be shifted by an arbitrary 
amount. 

2. When the variable of integration is changed from 0 to —0, the limits of integra¬ 
tion can be rewritten as 




Here we have used item number [[]. Since cos(— z) = cos(z), the only thing that 
changes in the integrand is the sign in the exponential. Hence, I (in) = I(—m), 
or the integral is only a function of the magnitude of m. 
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3. The integral is real, since I(m)* = I(—m ) = I(m). 


In addition, 


d m 1 
da m a + b cos (p 


(—1 ) m m \ 


1 

(a + b cos 0) m+1 ' 


(D.3) 


so that 


r^c^fr 1 

J 0 2i r ( a . + cos0) n J 


A (-!)("'-1) <*"‘-1 1 

(D.4) 

( n t - 1 ) ! ddp l ~ l Jo 27T J-jL dj + bj cos cp' 

Since bj ^ 0, we can factor out a TT^ i -w, and we relabel ^ —s 

7 ! ’ il J=l j, i ’ bj 

• dj. Then, the 

non-trivial integral is 


d 2 y im % j+ L, 


- r fr 2e ‘" 

7o 27 t 11 e 2 *'/’ _|_ 2aje 1 ^ + 1 

(D.5) 


Now we let z = so dcPe 1 ^ = dz/i , which allows the integral to be converted to 
a contour integral over the contour (7, a unit circle in the complex plane. Then, 


I({cLi},rn,N) 



d z Jm \-1 

2m 


N 


n 


2^ 


2 ^ + 2a jZ + 1 



d Z Jm| —1 

27 ri 


N 


n 




2z 

z j+ )(z 



(D.6) 


where 


Zj± = 


—a,j ± sign(aj) 



if \a,j\ > 1 


a j —(— % 




if \a,j\ < 1 


(D.7) 


Here (pi = arccos 
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We will only consider the case where a 3 > 1. In this case, only the only poles 
that contribute are the ones at z — Zj + . The contour integral is done by inspection 
to obtain 


I({ai},m,N ) f 

k =1 Jz ‘ 
N 

= £ 


r dz 

2 ~|m| 

* 9r 

TT 

L k+ 2vr* 

- 2fc+)(2 - 2fc_ 

) z 2 + 2a,-z + 1 


2z 


k+ 


N 

z k+ - z k _. n 

k= 1 ]=l,l 


2 z k + 


z k _|_ + 2 cijZ k + + 1 


Note that 

z k ± T 2dj^-i- T 1 "P 2nfc^A:± T 1 T 2(cq a k ')z k ± 2(cij o,k)^kz 
— z k - = 2 sign(a fc ) - 1. 

We define c*, and d k by 

Ct = 2 ’ 

** = 4+'- 


With the appropriate replacements, we obtain 

jv , ^ 

/({a;},m,IV) = ]J 


, i ^-k . , . ,, 

fc=i j=ijyfe 

Putting the fe/s back in using the a* —> we get 

c fc = sign(afc)^/a| - b 2 k , 


d k 


Cfc Ufc 


r*27T 


dcj) 

In 


tom 


N 1 
IT—-— 

a,- + 6, cos ( 
j=i J J 


n n i 

TT 1 

77 II &7 - 

fc=l Lfc j=l,j^k a j a kb k 


This expression is used in Eq. ( p. 4| ) get 


r*2?r 


d(j> 

—( 

2tt 


,* i(})m 


N i 

TT_1_ 

- L J j- (a 7 - + bj COS 4>) n i 


N 

n 


(_l)(m-i) d"i-i d 


N 


\ Ufc TT r 

L _! fni - ! ) ! 7 c * ,77 % - "7 


(D.8) 

(D.9) 

(D.10) 

(D.ll) 

(D.12) 

(D-13) 

(D.14) 

(D.15) 

(D.16) 

(D.17) 


(D.18) 
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Specificly, for IV = 1, we can drop the subscripts. For n = 1, 


-2tt 




d 


2tt a + b cos 4> c ’ 


icf)m 


(D.19) 


and for n — 2, 


r»27r 


d<f> 


1 d d j, d 

I g 27 t (a + 6cos0) 2 da c c 


, i(f)m 


f'2 = 


a + m c 


<r 


(D.20) 

(D.21) 


Since I(a,b,m ) is like an eigenfunction of (but not exactly, since the “eigenvalue 
is a function of a), we an use recursion to find that for general n > 2, 

r 2n dt 




1 


/ 0 27t (a + bcoscj)) r 


= f n I(a,b,m), 


fn 1 I / 2 /n—1 , fn—1 ) j 

n — 1 \ da 


(D.22) 

(D.23) 


Using this we get 


/i = l, 

a + \m\ c 
J 2 ^5 5 


h — 


3 a 2 + c 2 (m 2 — 1) + 3a|m|c 
2^ ' 


(D.24) 

(D.25) 

(D.26) 


The higher-order /(’s can easily be found from analytic iteration using a program like 
Mathematica, but apparently they do not have a nice analytic form. We note that in 
general, /, is multiplied by a factor of c -2 h -1 ). 

We now check the equation for situations that may be numerically unstable. First 
note that as c —» 0, the /’s are singular, which is due to the form of the original 
integral. Since the singular part is multiplicative, it can easily be treated. 

Instabilities appear when b is very small, in which case we have to take care of d 
carefully as it appears to be singular. Let b = 2 aS, so we can write 

Vl-4(S 2 -l x H 


dk — 


2h 


-S(\ + d 2 (l + 2 + <5 2 (5 + 5\U + 5 2 ...))))); 


|m| 
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This expression demonstrates that d k is not go as b l m l, but as & + H a s b —> 0. 

Now we analyze what happens where there is more than one term in the denomi¬ 
nator. We note first that 


dx 


f(x)g(x) 


E 


n\ 


—' m\(n — m)\ 

m— 0 ' 


d 

dx 


f(x) 


d 

dx 


g(x) 


(D.27) 


This means that Eq. ( P.18 ) can be written as 




,i(pm 


n 


27T 11 (dj + bj COS 4>) n i 


' N nj—1 

IIE 

k= 1 Lz=l i t = 0 . 


N 

E 


(—!)(”*-!-**) d ni - x ~ il f d, 


(ni - 1 - if)! da"* 1 \ c k 


x 


(—1)4) 

d i; I 

fr ill 

(h)' 

da^ ] 

11 6, ( 

J=l,j^k a j ak b k ) on _ 


(D.28) 


The product can be broken into two parts, one where l = k and another where l ^ k. 


f'2ir 


dcf) 


,icf)m 


N 

n 


2vr ^ (a,j + bj cos 0) n ^ 

N r N ni -1 

n e 


X 


E 

k =1 

X 

nk- r 

E 

-*fc=0. 

X 


,l=l,l^=k ii =0 

-!)(*«) d*‘ 


rfni-l-ii j' 4 
Ck 


{ni - 1 — if)! 


(b)! da, 1 


N 

n 


J ak b k 


^!—ifc) (jn k -l-ik ( d 


_(n k - 1 - i k )\ da? k 1 


1 

T 

d lk J 

1 - 

T —1 

(4)! 

da l k | 

11 bj ( 

J=1 ,&k a 0 ak b k ) on _ 


(D.29) 









































165 


Now we note that -f-— = 0, and so on, to get 


dai Ck 


2 vr 


n 


L _1 [flj + bj COS 4>) n J 

(-l) {r " 0 d n '- 1 

k =1 1=1,l^k 


(n, - 1 )! da? 1 


-i 


n —^ 

<j=i,j¥=k a i ak b k 


j ri k -l 

x z_^t f n 


Ck 


i-k —o 



1 

T 

R-' 

cf fc 

1 - 

1 — 1 

^k 

^ k • 

da l £ 

1 11 b :t ( 

lj=l,j^k a j ak b k ) on _ 


The first derivative term in Eq. (p.30| ) can be written as 


# i 

2 t r 


i(j)m 


TT-_ i _ 

11 (% + cos 4>) n i 


N , n k -1 

d k 


(—l) lfc d %k 
c k ^- nk ~ lk i k \ dcfc 


EjrE/. 

fc=l 4=0 


tv 

n 


k j=l,j¥=k \ a j ak b k 
Now, we consider the situation where all the b t have the same value, b. 
where k = N, the complicated part of Eq. ( |D.31| ) is 

, Ar . 1 d lN T-r 1 

h N (t N , N) = —r— -J- II - r— 

da£ “ (% - %) J 


E /ijv-i(*JV—i) N — 1)—-: 

I 1 IT — 7 


1=0 


X 


j(i N -iN-i) 


(lN — *JV-l)! 

1 


da% N ~ iN ~ l] (ajv-i - 


IJV 


^ (*tv - hv-i + ^w-i - 1 )! 
> hjv-i 1 77 -:-TTT-ttt 

i “ 0 - *v-i)K r hv-i - 1 )- 

1 


X 


(n tv —1 - a N ) nN - 1+iN ~ iN - 1 


N 

n 

1 =3 


(h — b-i + Tty-i — 1)! 


* 1 - 1=0 


(ii - i { _i)!(n I _ 1 - 1 )! (a z _i - oqv)^- 1 ^ *'-1 


x 


(*2 + ™1 - 1 )! 


i 2 !(ni-l)! (ai - a zv ) ni+i2 ' 


(D.30) 

(D.31) 
I 11 the case 

(D.32) 

(D.33) 

(D.34) 

(D.35) 

(D.36) 
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From this, we can cyclically permute (or shuffle in any order) the labels to get a 
similar expression for k ^ N. 

For definiteness, we consider a few special cases. If rij = 1 for all j except N, we 
find that 


h 2 (i, 1) 
h 2 (i, 2) 


Si, o 

(02 ~ dlY 2 ’ 

1 

(oi - a 2 ) i+1 ’ 


and 


, i'll— ^ ,0 

j "(a 2 -a 1 )(a 3 ^a 1 )-’ 

* 1 

^ 3 ( Z ’ 3 ' ) (ai - a 3 ) i+1 ( a 2 - a 3 y~i +1 ’ 

Using these we can write 


9k 


f 2n y± e i^m TT 1 

Jo 2 tC Ai + t )j cos fyrij 


n k -l 

fn k -i k hN (ik,k), 

ik =0 

JV , 

E «fc 

—9k- 

c k 


k =1 


(D.37) 

(D.38) 

(D.39) 

(D.40) 

(D.41) 

(D.42) 

(D.43) 


Approximations 


For N = 2, consider what happens when cp ~ a 2 . We can write a 2 = a, cp = a — <5, 
so that 


I(a 1 ,a 2 ,b,n 1 ,n 2 ,m ) 

= f 2n 1 _ 1 

Jo 27r (o 2 + 6 cos 0 — <5) ni (a 2 + b cos 0) n2 

= f 27T ^ c i<t> m _1_y^ (ni~l+j)! 

Jo 2i r (a 2 + 6cos0) ni+n2 (rci — l)!j! (a + bcoscpy 


^(a 2 - ai) J 
.i=o 


( n i ~ 1 + j)! 

(m-i)!j! /ni+n2+J 


d 2 

c 2 ' 


(D.44) 

(D.45) 

(D.46) 
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For N = 3, we could have a 2 ~ « 3 (or equivalently a\ ~ 03 ). In this case, 


I(ai,a 2 ,a 3 ,b, 1,1 ,n,m) 

OO /»27T 7 / -j -j 

= ^ l < D47 > 

OO 

= E (° 3 _ a 2 ) J -T(ai,a 3 ,&, l,n + j + l,m). (D.48) 

i=o 

When ai ~ 02 but ai 7 ^ a 3 , we get 

OO 

/(ai, a 2 , a 3 , 6 , 1 , l,n, m) = ^(a 2 — oi) J '/(ai, a 3 , 6 , 2 + j, n, m). (D.49) 

3 =0 


Summary 

We define 


/({a*}, 6, = 

f 2n d $ c i<t>m TT 1 

(D.50) 

7 o 27T -I-A (a* + bcoscj)) ni 


CM 

I 

rO | 

1 

t-H 

C3 

(D.51) 


V V«fc/ 

dk = 

( Ck - dk \ |m| 

V b k I ■ 

(D.52) 


We assume that rq > 0 and b, 7 ^ (). 
For N — 1, the integral is 


I(a,b,n,m ) = 

c 


(D.53) 


and for IV = 2, the integral is 


^ m—1 


/(a lt a 2 ,6, ni ,n 2 ,m) = (-1)’“ (E ijf Z f 


^n2 —1 


+ (-i) ni E /• 


(j + ni - 1 )! 


d? 


n 2 -j 


3=0 


jK n l - !) ! («2 - Oi)- ?+ni / C 2 


. (D.54) 
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For N = 2, and with ri\ = 1 the integral simplifies to 


I (a i, 02 , b, 1, n, m) 


1 di 
(o 2 - Oi) n Ci 


'n— 1 

E/. 

o=o 


n-j 


(o 2 - o,)J - 1 


d 2 

c 2 ’ 


When N = 3 and ni = n 2 = 1 the integral is 


I(ai,a 2 ,a 3 ,b, 1,1 ,n,m) 

1 d, 


+ 


do 


(o 2 - ai)(a 3 - ai) n ci ( 0 'i - a 2 )(a 3 - a 2 ) n c 2 

"n—1 i 1 


i =0 


y"~0 _ a 3)' ?+1 ( a 2 — a 3 )* • 7+1 


d 3 

C 3 " 


(D.55) 


(D.56) 
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Appendix E 

CHECK OF THE UNCROSSED APPROXIMATION 


“For every problem, there is one solution which is simple, neat and 
wrong.” 

— H. L. Mencken 


In this section, we want to check that how well the approximation 


A cross ~ A ladder Gq A ladder; (E.l) 

works. Since we are using a Hamiltonian theory and are interested in the potentials, 
we compare the potentials defined by 


Etbe:X — "^.9o ^ Go Across G 0 }g 0 \ 

EtBE:UX = -^9o 1 (G() A'l adder GcK \adder G q ) (J 0 ' . 


The notation used here is dehned in sections |2.3.3| and [2.3.4 . 
The TBE crossed potential (TBE:X) can be written as 


(E.2) 

(E.3) 


ErBE:X — ErBE:SX + ErBE:TX + ItBE:WX + ErBE:ZX, (E.4) 

where the potentials on the right-hand side are dehned in section [2.2.2| . Calculation of 
the TBE approximate uncrossed potential (TBE:UX) is straightforward, but tedious. 
We find that 


Vtbe:Ux(A; fei,^) 


Vtbe:Uxi(A; k\, p x ) + Vtbe:Uxi(A; p 1; fci) 

+ Itbe:UX2(A; ki,p{) + Vtbe:UX 2(A; p l5 fci) , 


(E.5) 
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where 


hTBEiUXl^; fei,Pi) — 


M\ 4 f d 2 q A 


E J J 2 ( 2 tt ) 3 


f 1 p j( x - z ) 6 ( z - y) 

lo "(x-z)z 2 (z -y) 

1 


x 


E-p 1 -k 2 -u (q l -p 1 )-uj (fei-qj) 


x 


Vtbe:UX 2(77; fei, Pi) — — 


My f d 2 q± 


E J J 2 (2 tt) 3 


E-Qi ~ k 2 ( fc i - Qi) 

+ {!<-;> 2 }, 


’ f 1 , J(x - z)0(y - z) 

Jo Z ( x ~ z)z 2 (y - z) 

1 

x- 33 ---- 

e ~Qi -P2 - ^”(Pi - Qi) 


(E. 6 ) 


x 


E-q i -k 2 -u (fei - g x ) 

+ {!■«—>' 2 }. 


(E.7) 


Now consider the azimuthal-angle and loop integrals for these potentials. The 
approach used is similar to that of Section |C.2| . Analyzing the potentials reveals that 
each of the terms can be written in the following schematic form, 


VTBE:UX,i(& + , k EP + ,P±) 

f°° q±dq± rl 


Jo 2 ( 27 t ) 3 ju 

Ii(k + ,k ± ,q + ,q ± ,p + ,p ± ), 

r 2ir 


dz J ( k + , q + , p + ) Ii ( k + , k ±, q + , q ±, p + , p ±) 
d<p q d(f) 


Iq Ai ti + B\ i cos 4> q + C\ j, cos 4> \A 2 + B 2 cos (f) q J 


(E. 8 ) 


(E.9) 


where 0 = —<j) q + (j) p , and the A’s, B’ s, and C’s may have dependence on k + , k±, p + , 
p l, q + = and q±] they are independent of the azimuthal angles. These factors 
can be easily determined for each potential by examining the forms of the original 
equations. The rotational invariance of the potential about the three-axis allows the 
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c f>k integration to be done trivially. The </> integral is easily done to obtain 


»2t r 


I(k + , k ±1 q + ,q ± ,p + ,p ± ) = - 2 tt / 


Vh 2 - of, 


/ 1 2 T B2 cos < 


(E.10) 


where 


at = A hi + B u cos 1 


(E.11) 


Further simplification is possible for Vtbe:UX 2 , since for that potential B 12 = 0, 

(27T ) 2 A 2 


I{k + ,k ± ,q + ,q ± ,p + ,p ± ) = 


spZ? - c?, 2V /y - w 2 (ai - Bi) 


(E. 12) 


The techniques discussed in the Section |C.2| are used to do the remaining loop inte¬ 
grals. 

The spectra for the OBE+TBE:SB+TBE:UX potential can be calculated and 
compared to the OBE+TBE:SB+TBE:X potential (recall that this is the same as the 
OBE+TBE potential), the modified-Green’s-function potential, and the ladder plus 
crossed Bethe-Salpeter equation. The spectra are plotted in Fig. |E.1| . The spectra 
for the TBE:UX, TBE:X and BSE all lie close to each other, which shows that that 
the uncrossed approximation is valid. 

This shows that the important approximation in the modified-Green’s-function 


approach is not the uncrossed approximation, but the addition of the extra interaction 
terms in Eq. (|2.102|) which serve to mimic the higher order interactions. 
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Figure E.l: The spectra for the potential derived from the TBE truncation of the un¬ 
crossed approximation (OBE+TBE:SB+TBE:UX, denoted in the figure by TBE:UX), 
the TBE potential (OBE+TBE, denoted by TBE:X), the ladder plus crossed Bethe- 
Salpeter equation (BSE(L+X)), and the modified-Green’s-function potential (MGF). 
Except for the MGF spectrum, all the curves lie close to each other. 
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Appendix F 

SPINORS 

“Everything should be made as simple as possible, but not simpler.” 

—Albert Einstein 

F.l Definition of the Spinors 

The Dirac equation is 

{ip - M)ip(x) = 0. (F.l) 

We can multiply the Dirac equation by {ip + M) to get 

0 =(ip + M)(ip - M)ijj(x) = (-f + M 2 )i)(x) = (- d 2 + M 2 )^(x). (F.2) 

This means that for an eigenstate of momentum, i/) p (x) oc e ±iPttXl1 , with p 2 = M 2 . The 
energy component of p /t is chosen to be positive. 

Looking at the positive energy part of (the part where the sign in the 

exponent is negative), we let 

^ p>+ (x) = «(p)e-^* M , (F.3) 

and for the negative energy part, 

if> p ,-(x) =v[p)e Hp ^. (F.4) 

The u(p) and v(p) are column vectors. They depend on only the spatial part of p, 
since the momentum is required to be on shell. The column vectors satisfy 

(j> — M)u(p) — (jt> + M)v(p) = 0, 
u{p)(J> — M) = v(p)(j + M) = 0. 


(F.5) 

(F.6) 
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From this point on, we will only be concerned with these momentum space represen¬ 
tations. 

To get a Hamiltonian version of the Dirac equation, we multiply by 7 0 to get 


(E — a ■ p — / 3M)u(p ) = 0, 

(F.7) 

(E — a ■ p + j3M)v{p ) = 0. 

(F. 8 ) 


Using the Bjorken and Drcll representation, we can write the Hamiltonian version 
of the Dirac equation as 


E-M -a p W xip) 
-cr ■ p E + M J y u B (p) 

We find that u B (p = wmxip), so 


0 

0 


(F.9) 


u(p) = c 
v(p) = c 


crp 

E+M 

crp 

E+M 

1 


x{p), 

va(p), 


where c is the normalization constant. 

The normalization we choose is u(p, A ')u(p, A) = 5\>\, so 


Sy\ = u(p, A ')u(p, A) 


= c 


crp 

E+M 


1 0 
0 -1 


1 

crp 

E+M 


A ')x(p, A) 


= icr 1- 


p 


(E + M ) 2 


x\p, A')x(p, A) 


2 M 2 + p 2 + E 2 + 2EM + M 2 — p 2 x/w 

= icl 2 ---— TTT^-— xHp+)x(p+) 


= c 


, 2 M 
E + M 


(/-; + m) 2 
x\p, A')x(p, A). 


(F.10) 

(F.ll) 


(F.12) 

(F.13) 

(F.14) 

(F.15) 


(F.16) 
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If we choose also xHPi A')x(p, A) 
so 


1, then we get that c 



where W 


E + M , 


w(p) = 



v(p) = 




crp 

w 

1 


va(p)- 


(F .17) 


(F .18) 


We can choose any two linearly independent vectors for x(p) which correspond to the 
two spins. These different choices lead to the different definitions of the spinors. 


Assorted Spinor Definitions 

“One day Alice came to a fork in the road and saw a Cheshire cat in 
a tree. ‘Which road do I take?’ she asked. His response was a question: 
‘Where do you want to go?’ ‘I don’t know,’ Alice answered. ‘Then,’ said 
the cat, ‘it doesn’t matter.’ ” 

—Lewis Carroll 


The spinors used by Bjorken and Drell [ p9| are polarized in the £ direction, and 
so the x’s are chosen to be 

1 
0 


A+l/2 — 


A—1/2 


0 

1 


The light-front spinors are defined to be |Tb 

1 


wlf (k, A) = 


y/Mk+ 

1 

y/Mk+ 


[M A_ + ( k + + a 1 ■ fc ± )A + ] x LFjA 

[A_(M + oA ■ k L ) + A + k + ] 


(F.19) 

(F.20) 

(F.21) 

(F.22) 


Xlf,a = 


Aa 

0 


(F.23) 
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where xx is the usual Pauli spinor, and the A± are the spinor projection operators 
defined in Appendix |A|. We find that 

u LF (k, A) = -^=== Xlf,a [ A + M + A -( k+ + a± ' k± )] 

= ^== x[f,a [( M - ' k± ) A + + k+A ~] 

Note that these spinors are normalized to satisfy A')wlf(&, A) = Sy\. 

For hclicity spinors, we choose the eigenvectors of the helicity operator (cr • p) as 
the y’s. In particular, (p ■ S)n(p, A) = hu(p, A), where h = 2A. This choice allows us 
to write 


(F.24) 

(F.25) 


u(p, A) 



(F.26) 


where (cr • p)x\{p) = h\x{p)- 

One useful feature of the hclicity spinors is the simplification of the cr • p factor 
living in the prefactor of y in the definition of u. This is an important feature: all of 
the directional dependence lives in the y alone. (For the normal spinors, the y are 
independent of the momentum, but the prefactor has the dependence.) 

Now, we proceed to calculate Xx{p)- This will be done by taking two steps: first 
calculate yA(?), then rotate that to an arbitrary angle. 

The first step is easy. When p = z, the hclicity eigenvalue equation is 


(cr • z)x\(z) = a z x\(z ) 


(F.27) 


= hxx(z). 


(F.28) 


This means that yA(z) are just the usual vectors for the spin-1/2 basis, 


X+l /2 (Z) 
X—1/2 (z) 



7 


(F.29) 

(F.30) 
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Now for the second step. To rotate, use the Euler angle representation of the 


rotation operator [127 


= e ~ iaJy e~ i/3Jy e~ i7j! . 


(F .31) 


This operator allows us to write 


u(p, A) = R^ 6il u(p x z, A) 


= R, 


' 0 , 0,7 


w 


2 M \ hp 
w 




(F.32) 

(F.33) 


Note that we are free to choose 7 , since all it does is give a phase to 7 . The traditional 
approach is to set it equal to —(ft. This is what Brown and Jackson do [ 0 ], and in 
general it is the best thing to do, since it gives the correct rotations starting from 
a state pointing in an arbitrary direction. However, since our state just points in 
the 2 -direction, this angle just gives a phase. It will turn out to be useful to just 
set 7 = 0. This choice simplifies the equations we will use later, but be warned 
that these kets transform differently from those that Machleidt |5?| and Brown and 


Jackson use [||, but they are closer to what Rice and Kim [|128|1 use. Note that for 
totally symmetric potentials, both (ft and (ft' can be set to zero, so this distinction is a 
moot point. However, here (and in Rice and Kim’s work), we cannot set them both 
to zero, and we find that our convention is much easier to work with. 

So we have 


u(p, A) — R<j>,e, 0 


W 
2M 


hp 

W 


X\{z). 


(F.34) 


The operator will not affect the prefactor. Since the spin and the momentum rotate 
together, the factor of er ■ p remains unchanged. The rotation operator only hits 7 , 
so that 


Xx(p) = Rf.e. oXa(S)- 


(F .35) 
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The x is an SU (2) spinor, and the SU (2) representation of the angular momentum 
operator J is cr/2. Note that of = 1, and that in general if A 2 = 1, 


e i9A = 




i2n 

^ (2n)! ^^W2n+1)! 

n =0 V ’ n =0 v ’ 

= cos 9 + iAsmO. 


(-1 ) r 


0 2n+l iA 


(F.36) 

(F.37) 


Thus, 


. , , e~ ia ! 2 0 

'IOlJz _ | 

0 e +iQ / 2 


P ■ ■ P 

= cos — — iGy sm — = 


-iOJ x ® ■ 0 

e = cos- to x sm - = 

2 2 


cos f — sin | 
sin f cos | 


cos | — i sin | 
—i sin | cos f 


This means that R is 


R(f>,e, o — 


e -# /2 C2 _ e -#/ 2 S2 
e +i<t>/ 2g 2 e +i0/2 C2 


where c 2 = cos | and s 2 = sin This gives 


c 2 e 


Xa(^) = ^,0,oXa(2) = Xa (p) = { 


-i<f>/2 

+i(j>/2 


s 2 e 
-s 2 e -^ /2 


c 2 e 


+i(p/2 


if h = +i 


if h — — 1 


(F.38) 

(F.39) 

(F.40) 


(F-41) 


(F.42) 


F.2 Two Helicity Spinors in the Center of Momentum Frame 


We now address the problem of two particles with spin. We work in the center-of- 
momentum frame, and consider the relative momentum between the particles. 

To start with, consider having the relative momentum pointing in the z direction. 
The first particle has momentum pz and helicity Ai. The helicity spinor for this 
particle is just what was discussed in the previous section. 
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The second particle will have momentum — pz and helicity A 2 . We can relate this 
to a spinor with the momentum pointing in the z direction via a rotation, 


u(-pz, A 2 ) = R^nflUipZ, A 2 ), 


(F.43) 


where 0 is a free parameter because it has no meaning when pointing along the z-axis. 
As argued in the previous section, this can be recast as 


Xa 2 ( z) R<f>, it, oX,\2 (^0 • 


(F.44) 


Using Eq. (F.41), we can write 


R(f>,ir, 0 


If we choose (j) = tt, this becomes 


0 - e~^ /2 




0 


1 | i 


0 i 
i 0 


so that 


(F.45) 


(F.46) 


XXoA-z) = b\-A 2 (2j 


(F.47) 


When the relative momentum points in an arbitrary direction, we find that for 
particle 2, 


Xa 2 ( P ) R(f),o,oX\2 ( z) 

iR<f>,6,oX— A 2 (•2<) 

= ix-x 2 (p)- 


This is rather simple. 


(F.48) 

(F.49) 

(F.50) 
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F.3 Helicity Spinor Summary 


We have 


u(p, A) = 


W 
2 M 


f 


Xa (p), 


and 


c 2 e 


Xa(p) = < 


—i<j>/2 
-\-i(f >/2 


5 2 e 
—s 2 e -1 ^/ 2 


c 2 e 


-\-i(j >/2 


if h = +1 


if h = -1 


(F-51) 


(F .52) 


where c 2 = cos |, s 2 = sin |, / = and h = 2A. 

When there are two fermions in the center-of-momentnm frame, we can define 
0 = 0! and 6 = 6i and for particle two 0 2 = 7r + 0 and 9 2 = it — 9. This means that 


“(Pi» A 0 = \/ 7777 . | Xi,Ai(p), 


2 M 


fi 


where * = 1,2 and 


(F.53) 


Xl,Ai(p) = Xa iip), 
X2,A 2 (p) = *X-A 2 (P)- 


(F.54) 

(F.55) 


F.4 Breakup of Dirac Spinor Matrix Elements 


We are interested in matrix elements with the form 


u (p', X')Tu{p, X) = 


Vww t 


2 M 


Xv(p') 




r aa Fab 
Tba Fbb 


vWtf 

2 M 


Xa'^aaXa + Xy^ABXxf 


-xlrBAXA/' - Xy^BBXxf'f) ■ 


1 

/ 


Xa(p) 


(F.56) 
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The elements of the T a fe matrices will be either the identity matrix or Pauli sigma 
matrices. 


For example, we have 


«(P , A )u(p, A) - 2M (1 / /), 

(F .57) 

w(p', A')7 °m(p, A) - xIxa 2M (1 + f'f ), 

(F.58) 

«(?', A')7 5 «(P, A) - 2M (/ / ), 

(F.59) 


and some more complicated ones, 

u{p', A')t *u(p, A) = (f + f) ’ ( F - 60 ) 

u(p', X')a^ u i(p' — p) v u{p, A) = 2 Mu(p\ A')yA) 

- ip' + p)^u(p', \')u(p, A). (F.61) 


We used the Gordon decomposition (Section |F.5| ) to simplify the last equation. This 
simplification of the tensor allows us to only worry about two tricky matrix elements 
for vector mesons. 

With the Gordon decomposition, we can further simplify the vector meson cou¬ 
pling by using the Dirac equation. Note that the expression for T is only valid in this 
context. 


r < = 0 + £) y ‘- f 2 F (p:+pi) '‘- 

Uy = u(p[, Ai)Ti u(p lt Ai ){-g^)u{p' 2 , \' 2 )T 2 u{p 2 , A 2 ) 
= (l + ^ u(Pi, Ai)7 A *w(P 1 , Ai)(-^)n(p2, A(>)7 


'i + 'W 1 


(F.62) 

(F.63) 


'm(p 2 > A 2 ) 

+ ( 1 + ^ A 'i)(/ 2 +^ 2 Mfh> Ai)w(P 2> A' 2 )u(p 2 , A 2 ) 

ff A 2)(/l + A) W (P2, A 2) 

! l 

(Pi + PiT(P 2 + P2) f ij^u{p[, Ai)w(p 1; Ai)u(pa, A 2 )w(p 2 , A 2 ). (F.64) 
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The Dirac equation can be used to write ^ 2 w(pi, A) = (2E'y° — M)u(pi, A), and also 
note that 


(. p'i + PiY l (P2 + P 2 V = 2 (E' + E) 2 - (p' + piY{p[ + Pi) M 

= 2 (E' + E) 2 - 2 M 2 - 2 E'E + 2 p' ■ p. 


(F.65) 

(F.66) 


This means we can write 


Uy= ^1 + 0 «(pi, Ai) 7 '*«(p 1 , Ai )(-g^)u{p' 2 , \' 2 )Yu{p 2 , A 2 ) 

-2(1 + -) ~u{p[, Ai)u(pi, Ai)m(p' 2 , A' 2 )u(p 2 , A 2 ) 

/A fE' + E^, ^ 


+ 1 + - - 


gjg M 


u{p'i, KhHpi, Ai )u(p' 2 , A i)«(p 2 , a 2 ) 


+ ( 1 + -^-^^«(pi, A / 1 )u(p 1 , Ai)«(pi, A' 2 )7°u(p 2 , A 2 ) 


9 ) 9 m 

f\ 2 (E' + £) 2 - M 2 - E'E + p' p 
~g ) 2 M 2 

x «(pi, A , 1 )u(p 1 , Ai)m(p' 2 , \ 2 )u(p 2 , A 2 ). 


(F.67) 


We are interested in the matrix elements for OME and TBE exchange potentials. 
Thus, we want to look at combinations of two uu matrix elements. To save space, we 
will write 


(Ai, Ai|A]., A 2 ) = (x 1 A 'X-a 2 ) , 


(Ai, Ai| o - ! • <t 2 |Ai, A 2 ) = (Xx'yXxi) [XLx’X^-^J • 
We can write the matrix element for the scalar exchange as 
Us = «(pi, Ai)«(p!, Ai) x u(p' 2 , a ' 2 )u(p 2 , A 2 ) 


W'W 
4 M 2 


1 - 


h' 2 hip'p 

W'W 




and the matrix element for the pseudoscalar exchange as 
Ups = u(p[, A / 1 )i 7 5 u(p 1 , Ai) x u(p 2 , \ 2 )i~/ 5 u(p 2 , A 2 ) 


W'W fhip h\p' 
AM 2 \W ~ ~W 


h 2 p h' 2 p' 

w - w 


(A' 1; A' 2 |Ai, A 2 ). 


(F.68) 

(F.69) 


(F.70) 


(F.71) 
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Now for the vector meson. Note that one part is the same as for the scalar meson. 
The other part is more complicated, the vector-vector part: 


f/ vv = u(p[, AibMfh, Ai )(-g l u,)u(p' 2 , X' 2 )Yu{p 2 , A 2 ) 


W'W 
4 M 2 L 


1 + 


h\ hxp’p 
W'W 


h\p h'pp' 

W + ~W 


^ + ^)(A;,A' 2 k 1 .^|A 1 .A 2 ) 


(F.72) 


This means that the full uu matrix element for a vector meson with a non-vanishing 
tensor coupling is 

2 


U v = ( 1 + - ) U vv + - 
9/ 9 

= ( 1 + -) f/vv + - 

9/ 9 


' / /3 -E'E + p'-p 

g\ 2 2 M 2 

• / /3 -E'E + p' 3 p 3 

g\2 2 M 2 


f/s 

C/s 


/ \ 2 P'mPm 


g ) 4 M 2 


£4 . 


(F.73) 


(F.74) 


We can write this out as vector-tensor, and tensor-tensor parts. The vector-vector 
is just t/ vv . 


C/vt. = 


— [C/vv + U s ] 

9 

2/ VF'VF 
~ ~g 4 M 2 


U tt = 


7\ 2 

■(;)' 


ijiwp h 2 h 2 pp\ , , 

r+ L v WW + W'W ) ' 11 i! 27 
/hip /^A /7 2 p fo 2 p / N 
\lh F; 


y (Aj, A ' 2 \(Ti ■ 0-21 Al, A 2 ) 

1 / E'E p' ■ p \ 

U ™ + -2{ 3 -JP + ^) U - 

W'W 


W'W r , , 

^[{A'nA^nA,) 


1 + 


h\hd 


P 


— . K x VFW 

3M 2 -E'E + p'-p / 

l 1 


i + 


h' 2 h 2 p‘ 

W'W 


p 


+ 


2 M 2 
/ h\p h[p' 

+ w 7 


ii'ii 


h' 2 h 2 p‘ 


p 


W'W 


(F.75) 

(F.76) 

(F.77) 


if + lw) ^A^I^-^lAnA 2 ) • (F.78) 
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F.5 The Gordon Decomposition 

“Nothing shocks me. I’m a scientist.” 

—Indiana Jones 

Now, let’s look at the Gordon decomposition. We start by looking at 

ia^ = ^ ( 7 V - YY) = 9*™ ~ YY- (F.79) 

Then, 

u{p', A ')a^i{p' - p) v u(jp, A) = - u(p', X^a^ip^ip, A) - u(p', A ')ipY Ufl u(p, A) 

= u(p r , A ')Yi>u{jp, A) + u(p f , A ')i>Y u iPi A) 

— (jp r + pYu(jp', A ')u(p, A). (F.80) 

Using the Dirac equation, we see that 

u(p f , A ')<j^i(p' — p) v u{p, A) = 2 Mu(p', \')Y u (Pi A) — (p + pYu{p' , X')u(p, A), (F.81) 


which is the Gordon identity. 
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Appendix G 

ROTATION MATRICES 


“There comes a time when for every addition of knowledge you forget 
something that you knew before. It is of the highest importance, therefore, 
not to have useless facts elbowing out the useful ones.” 

—Arthur Conan Doyle 


G.l Relations Between Rotation Matrices 


Brown and Jackson give a good introduction to this topic ||. We start out with the 
rotation operator, parameterized in terms of the Euler angles, 


R a01 = e - iaJ *e- if}J ye- ilJ \ 


We can take matrix elements of this, to get 


(G.l) 


(J'M'\R af3y \JM) = 7 ) 


(G.2) 

(G.3) 


Note that the lower indices in d are in the same order as in T>. This is different 
that the convention in |9|,^], but agrees with [127-129]. However, we also define the 
matrix elements slightly differently, and the combination of those two effects makes 
the calculated Jacobi polynomials the same! 

An alternative interpretation of the TTs is that they are the overlap of JMX kets 
with 9(ft\ kets. To see this, generalize Sakurai’s approach for kets without spin ||127 . 


We start out with |zA). This is a ket for a particle pointed upwards, with A units of 
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spin in the ^-direction. This ket can be rotated with R to get 

R^\zX) = 1 0<f>\). (G.4) 

Brown and Jackson choose 7 = — qf>, however our convention is to choose 7 = 0 and 
get 


R<t>e ol^A) — |$<M) (G.5) 

Applying a bra to the previous equation, we obtain 

(J'M'X'\9(f)X) = (J'M'\'\Rw 0 \z\) = (J'M'\'\Rwo\JM\")(JA4\"\z\). (G.6) 

In the last bra-ket, we must have A = A", since the helicity cannot be changed. 
Furthermore, since A is the J z quantum number for the final state, we also have 
A = M. These facts allow us to write 

(JMA"\zA) = \z) = <t>) = ^ , (G.7) 

which means that 

(JM\\6<t>\) = = lv J uk (4M). (G.S) 

Since the (j) dependence of the T>'s is essentially trivial, we can replace the 0 with a 7 
to get 

= \J (G.9) 

and 


(G.10) 
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Using these relations, we find that 

'D'm'M { a 0l) = 7 ~ ft ~ (G.ll) 

^•m'm(P) = (G.12) 

(<t>, 0, X) = JIKIyjmM = MzJW py (cos 0)e lM * (G.13) 

<4,(« = l/ (j+Af|l ^( cos ' 3 )- (G.14) 

do 0 (/3) = Pj(cos/3), (G.15) 

0) = ^M'M = (G.16) 

Now, let’s use this. Our first thing to do is look at is the addition theorem 

A) 7i) 

M" 

= ^(JM|P a2A , 72 |JM")(JM"|P QlA)Tl |JM') (G.17) 

M" 

= (JM|P a2)/ j 2)72 P ai!i a li71 | JM') (G.18) 

= (JM|P aA7 |JM') (G.19) 

= ®MM'( a )A 7)- (G.20) 


We used the group property of the rotations that the product of two elements gives 
another element. The expressions for (a,/?,7) in terms of the other angles is rather 
complex in general. 
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We can also look at the orthogonality relation, considering the Vs as overlaps, 


2J+ 1 

87T 2 


r»27T /*7T 


da / sin (3 d/3 / d~V':, (afJ'y)V^,.(a/3j) 


2 -W^^ 2jT + l2jT + l l da L SmPd0 l “ 7 

((J 1 M(M 1 |a/3M 1 )e- i7Ml )*(J 2 M'M 2 |a/3M 2 )e- i7M2 


= (2J+ 1)< 


2Ji + 1 2J 2 + 1 


dl c +i'y(M 1 -M 2 ) 


J <Kl{J 2 M' 2 M 2 \SIM 2 ) (OMjI J V M[M X ) 

= (2J+7 V^yTTsjrT' 5 " 1 " 2 

J dn( J 2 M' 2 M x \ rtMi) (OMj I J t M[Mi) 

= (2J+ 1 ) \ZW+~ 12 J 2 1 + 


(G.21) 


So the conclusion is 


2J+ 1 
8n 2 


da / sin (3dfi / d^V^, *(a(3^)V J ^, M (a(3^) 

) Jo Jo 1 2 

= < 5 miM 2 < 5 m(M'< 5 jiJ 2 - 


(G.22) 


Using this, we find 


/” 2.7 + 1 

' 4?f -</’) = S Ml M 2 Sj lJ2 , 

9 / i i 

X —t—Vmw «+'. +)+,„+,», -0) = <5(f2 - sp, 


(G.23) 

(G.24) 


r 1 2 7 + 1 

j dcos9 — dj^, M (6)dj^, M (9) = 

<W 2 > 

9/^-1 

-+- d J MM'( 0 ') d MM'( 0 ) = S(cOs9~ COS 0') 


(G.25) 

(G.26) 


Now, for two particle state in the center-of-momentum frame, we can write the 
ket as |p6 l 0AiA 2 ). The total helicity of this state is A = A x — A 2 , which is easily seen 
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if we choose the momentum to be in the z-direct ion. Since this ket has total helicity 
A, we can transform it using 


\pJMX\X 2 ) — J dO |p$0AiA 2 )(p$(/>AiA 2 |pJMAiA 2 ) 
= I dfl\p9cj)XiX 2 )(9(j)X\JMX) 


2J+1 


47T 


dnV J MX *{9<f) 0 ) \p9cP ArAa), 


(G.27) 

(G.28) 

(G.29) 


where A = Ai — A2. 

We can break apart Eq. (|G.29|) as we did for the orthogonality relations, 


|pJMA 1 A 2 ) = [ dnV J MX (cj ) ,9,0)\p9^X 1 X 2 ) 


47T 


2J + 1 


f27T 


dcos9 d J MX {9)\p9MX\X 2 ), 




|p0MAiA 2 ) = 


#e^ M |p£<MiA 2 ), 


V2n Jo 

and the inverse relations (obtained using the orthogonality relations) 

1 


\p9<j)X±X 2 ) = 




E' 

M 


i(j)M 


\p9MX\X 2 ) 


J,M 


E 0)|pJMA,A 2 >, 


47T 


\p9MXiX 2 ) = Y J \r J ^ d M^ 6 )\P JM XiX2). 


(G.30) 

(G.31) 

(G.32) 

(G.33) 

(G.34) 

(G.35) 


G.2 Jacobi Polynomials 


“Logic, like whiskey, loses its beneficial effect when taken in too large 
quantities.” 

—Lord Dunsany 

The Jacobi polynomials are neat guys. They are defined by 


(J, M \Rqqq\J, M ) — d J M i M {9). 


(G.36) 
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Note that since the Jacobi polynomials are real, 


d J M , M (e) = ((J,M'\R oeo \J,M))* = (J, M\R 0 _q 0 \J, M') = d J MM ,(-0). 

There are other symmetries that the Jacobi polynomials obey. To see 
turn to Schwinger’s model of angular momentum. In this model, we have 


I J,M) = 


, » J+M / , \ J—M 

al_) ( a_ 


VJ + M y/J-M 


|0,0>, 


where a± essentially adds in one half unit of spin, and changes M by ±|. 
What does this state do under rotation about the y- axis? Well, 


R-ooo | -h M) — 


j. | \ J+M / | | \ J—M 

Robo^+Rqoo ) ( Roood—Rooo 


-Roeo\0, 0). 


y/J + M y/J-M 

The vacuum is unaffected by rotations. Now, Sakurai shows that ||127 


9 9 

Rqoq^^Rqqq c, . cos — i sin —. 

This can be used to get an explicit expression for d. 

Consider now what happens to the state under a rotation by 7r. We get a\ 
Thus, 


RonO | J? M) — 


,t 


J+M 


,t 


J-M 


y/J + M y/J-M 

J-M / , \ J+M 


= (-l) 


J-M 




| 0 , 0 ) 

t Y 


v - 


Jj-^M +J + M 


| 0 , 0 ) 


= -M), 


which gives a nice relation for Clcbsch-Gordon coefficients 

(J 1 M 1 J 2 M 2 \JM) = (- 1 ) j i+ j 2+ j - m 1 - m 2- m (J 1 - M x J 2 -M 2 \J -M) 
= (-1 - M X J 2 — M 2 1J — M). 


(G.37) 
them, we 

(G.38) 

(G.39) 

(G.40) 

t * 

(G.41) 

(G.42) 

(G.43) 

(G.44) 

(G.45) 
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Using this result for the Jacobi polynomials, 

d J M ' M (n -9) = {J, M'\R 0n -eo\Ji M) (G.46) 

= (J,M'\R 0 _ e0 R 0n0 \J,M) (G.47) 

= (J, M'\Rq_0 O \J, —M) = (-1 (G.48) 

= (J, —M'\Rq_0 O \J, M) = (-1 ) J - M 'd J _ M , M {-9). (G.49) 


We can rotate by a full 2tt, and get 


+ 9) — d J M , M {0) 

— {J, M'\R 0n0 RQQ 0 R 07r0 \J, M) 

— ( (J, M'\Rq0q\J, M) 


(G.50) 

(G.51) 

-1 ) M ’- M d J _ M ,_ M (6). (G.52) 


We can also use a trick from Ref. 1291. Note that 


= e T 'e Tl ^ v e 


(G.53) 


The J x exponential is calculated in Eq. ( h.4Cj), and for an angle of n it i 


is 


e inJx = i 


0 1 
1 0 


(G.54) 


so under this rotation at —> fat, so 


e ~™ J x | j, m) = (f) 


J+M+J-M | 


J,~M) = (-l) J | 


(G.55) 


Thus, 


d J M ' M (9) = (J,M'\e +i0Jv R o _ eo e- i0Jy \J,M) 

(G.56) 

1 

S' 

o 

1 

c? 

% 

1 

sr 

(G.57) 

— d J _ M ,_ M { 9). 

(G.58) 
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Also, we can use the x-rotation to look at 

0) = {JiM'\Rqk-oo\J,M) (G.59) 

— {J, M'\Ro-eoRoiro\J, M) (G.60) 

= (J, M'\R O -0 O \J, -M) = (-1 (G.61) 

(- l) J - M '(J,-M'\R 0 _ eo \J,M) = (-1 (G.62) 

Summarizing what we have so far, 

= (G.63) 

= d J _ M ,_ M (-0 ) (G.64) 

= (-1 (G.65) 

dMM'i^ ~ d) = (—1) J M d J M ,_ M {—6) (G.66) 

= (-1 (G.67) 


What sort of fundamental symmetries do we expect the Jacobi polynomials to 
obey? 


1. Under 6 —> —6. Combining Eqs. ( G.65 ) and (|G.64|) gives 


d J M ' M {0) = (-1 ) M '~ M d J M , M {-e). 


(G.68) 


2 . 


3. 


Under m ml. Combining Eqs. ( |G.63|) , (|G.64| ), and (|G.65|) , we get 

d J M ' M {0) = (-1 ) M '~ M d J MM '{9). 


Under m —> — m and mn! —> —m’. Looking at Eq. (|G.65| ), 

dUM = (-1 ) M '- M d J _ M ,_ M {9). 


(G.69) 


(G.70) 


4. Under 6 —> n±9. Apparently this is not a symmetry that the Jacobi polynomials 
have. This must be accompanied with a m or m' being negated as well. 
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5. Under 9 —> n — 9 and m —> —m. Combining Eqs. ( |G.66| ) and ( |G.68|) gives 

^M'iw( 7r - 9) = ( — 1) J M d J M /_ M (9). (G-71) 


6. Under 9 ^ n — 9 and m' —> —ml . Combining Eqs. ( G.67 ) and (|G.68| ) gives 


Q) — ( 1) J M 


(G.72) 


As a closing note for this section, Sakurai gives an expression for d J M , M calculated 


using Schwinger’s angular momentum model ||127 


d J M , M W ) = + M)\(J - M)\(J + M')\(J - M')\ 


x 


EM) 


fc-M+VT 
k-M+M' b 2 b 2 


k=k„ 


k\ (k — M + M')\ 


x 


„J+M-k 

c 2 


J-M'-k 


(J + M-k)\ {J-M' -k)V 


k m;n = Max(0, M — M'), 
A)max = Min( J + M, J — M 1 ), 

where s 2 = sin(/3/2) and c 2 = cos(/3/2). 


(G.73) 

(G.74) 

(G.75) 


G.3 Helicity Applied to Potentials with Full Rotational Invariance 


Now, we want to look at matrix elements of the potential 
(p'X'A V|pA,A 2 ) = 


2 J 1 1^ j/ *, ,, s .—/2 J A 1 


At: 




J',M' J,M 

(j>'J'M , \' l \' 2 \V\pJM\i\2). 


An 


0) 


(G.76) 


But if the potential has full rotational invariance, that means that it is a scalar. We 
need to use the Wigner-Eckart theorem, which states that 


(c,'J'M'\Tl\aJM) = (Jj- |=A 


(G.77) 
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For a scalar potential, V = Tg, so 

(a'J'M>\V\aJM) = (G.78) 

V 2 J + 1 

Thus, the potential matrix element is independent of the M value, and is diagonal in 
J. 


This allows us to write 


(p / a;a / 2 |f| p a 1 a 2 ) 


E 


2J+ 1 
47T 


. M 


(p / JA / 1 A'|\/|pJA 1 A 2 ). 


(G.79) 


The M dependence of the matrix element has simply been omitted, since there is no 
dependence. 

Recall the addition theorem for the rotation matrices, 

{afa) = E (ai/3i7i)T ) ]( f ,/ M (a 2 /3 2 7 2 ), (G.80) 

M" 

where 


Ra/3j 7 ?, 


ai/?l7l 77a2/?272 


(G.81) 


e ~iaJ ze ~i0J y g-i^M _ e ~i(3iJ ye ~i{'n+cx2)Jzg-ifoJy^—iyiJz 


(G.82) 


Using the addition theorem, we get 


M M 

= ViM,6",-M). (G.84) 


In general, (j)" ^ 0 2 , and 6" ^ 6' — 6. However, we will be able to make some 
simplifications. 

The sum rule allows us to write 


2J + 1 


(p'A'^IVpAOs) = e-‘« y e‘« A V ___4 A (r)(p'JA' 1 AiV|pJA 1 A 2 ). (G.85) 
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All we are interested is the matrix element buried inside the sum, so we invert the 
equation to get 

{p'J\' l \' 2 \V\pJ\ 1 \ 2 )=2'K J d cos 9" dy X (9")e i ^" x 'e~ i ^ x (p , X l 1 X , 2\V\pXiX 2 ). (G.86) 

Note that there are three free parameters in the integrand. We can choose 0 = 0' = 
6 = 0. This means that (j)" = (j)" = 0, and 6" = 6'. Thus, 

(p'JX'^lVlpJX^) 

= 27 t J dcosd'dy X (6')((p'sin 6', 0,p'cos 6') X , 1 X , 2 \V\pzX 1 X 2 ), (G.87) 

which is what is advertised in Machleidt’s Eq. (E.41) in Ref. [ p>B|j . 

G.4 Helicity Applied to Potentials with Cylindrical Invariance 

Now, we want to look at matrix elements of the potential 

(p'a'^MpAjA,) = V Jliillvs.pw,?, 0)]r 

J>,M 1 v J,M 

(p' J' M' X[X' 2 \V\pJ M XiX 2 ). (G.88) 

Now we consider that the potential operator is not a scalar, but it still conserves M. 
In terms of spherical tensor operators, we can write 

r=£ (G.89) 

J 

This does not bode well for an approach like the one that was used in the previous 
section. 

However, we can write 

(p'WV|pA 1 A 2 ) = £ 

(p'J , MA , 1 A' 2 |y|pJMA 1 A 2 >. 


2T + l^j 


47T 


0 ) 


(G.90) 
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Since the M is all that is conserved, it makes sense to ignore J for now. So we change 
basis to 

(p'KK\V\p*M = Y ^e iA * M (j/d'M\ , 1 \ , 2 \V\p6M\ 1 \ 2 ) i (G.91) 

z ' Z7T 
M 

where A <f> — (j)' — (j). As in the previous section, all we are interested in is the matrix 
element buried in the sum. We invert the equation to obtain 

/»27T 

(p , e , MX , 1 X , 2 \V\p9MX 1 X 2 ) = / (G.92) 

Jo 

There is one free angle in the integrand. Setting (f) = 0 for simplicity, to get 

p2n 

{p'e , M\' 1 X l \V\pdM\ 1 \ 2 )= / d(jj (T l<t> ' M {p'\\ X 2 1 V>3, PaX^X-i) ■ (G.93) 

Jo 

How does this relate to what we got before, with the J, M matrix elements? We 
need to change basis once more, to get 

In J' i r 1 

(p / J'MA / 1 A'|H|pJMA 1 A 2 ) = j ^ d cos O' d J MX , (0') 

x \/ 2J 2 1 y dc ° s0d MA(0) 

x (p9 , M\' 1 \ , 2 \V\p9M\ 1 \ 2 ). (G.94) 


G.5 Summary of Transformations 

“Oh, I have now a mania for shortness. Whatever I read - my own or 
other people’s works - it all seems to me not short enough.” 

—Anton Chekhov 

First of all, we would like to be able to transform to | JMLS). In particular, we 
would like 


(JMSX\JMLS) 


dQ (JMS\\6<j>S\)(6<j>S\\JMLS). 


(G.95) 
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Looking at the first term, 

(jMS\\e<j>sx) = 


j 


An 


9-> 0)- 


(G.96) 


The second term takes a little more work. 


(i 9(f)S\\JMLS) = ( mSX\R\ eo \JMLS) 

= (00SA I JM'LS) ( JM'SL | R\ e 0 1 ,JMLS) 


M' 


Y (00S\\LM l SM s )(LM l SM s \JM'LS) 

x ((JMSL\R^ 0 \JM'LS))* 

Y Y LMl { 0, 0)S\,M s (LM l SM s | 9 ,0) 

M l ,M s ,M ' 


E 

2L + 1 
47T 


2L + 1 

47T 




(L0RA|JA>P^ a *(<M,0) 


(G.97) 


Combining these, we get, by virtue of the orthonormality of the P’s, 

2J + 1 


(JMSA|JMLS) = \ j L + | <L05*A| JA) 


2J+ 1 
2L + 1 
2J+ 1 


47T 


dfi 2^,(0, 0 ,O)Z>ma (0,6>,O) 


(L0EA|JA). 


(G.98) 


Using this, we obtain Machleidt’s equation Eq. (C18) in Ref. |58 


I JMLS) = Y, V^7TT <£0;SA I JA>|JMSA> ' 

A 

For a one particle state, we start with \9cftSX). From that, 


(G.99) 


f27T 


|0MSA) = 




d(t)e iM<t> \9(t)S\), 


I2J + 1C 1 

\JMSX) = \/—I dcos9 d J MX (9) \9MSX ), 


mi 


| JMLS) = \/ 2L + 1 {L0SX\JX) \ JMSX). 
2 J I 1 


(G.100) 

(G.101) 

(G.102) 
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Furthermore, for two particle states, we have I^SiAiS^A^. However, that is rewritten 
as 


IdcpSXS^) = ^(00S , iAi5 2 A2|005 , A5 , i52)|005 , iAi5 2 A2) (G.103) 

Ai A2 

= ^(OOS' 1 A 1 5 2 A 2 |i2^ o i^|OOS'AS' 1 5 2 )|005' 1 A 1 5 2 A 2 ) (G.104) 

Ai A2 

= ^(00^ 1 Ai^ 2 A 2 |00^A5i^ 2 )|^ 1 A 1 ^ 2 A 2 ) (G.105) 

Ai A2 

= = Ai S 2 M 2 = -A 2 | SM = XS 1 S2}\9(j)S 1 X 1 S2X2) (G.106) 

A1A2 

= J](5 , iAi5 2 -A 2 |5A)|005 , iAi5 2 A 2 ). (G.107) 

A1A2 


Now, the \9(pSXSiS2) can be used to get 

1 /*27T 

ftMSXS^) = —= / #e iM ^|d05A5i^ 2 ), (G.108) 

v2vr Jo 

\JMSXSiS 2 ) = I dcos9d J MX (9)\9MSXS 1 S 2 ), (G.109) 

\JMLSS 1 S 2 ) = \ 2L + 1 {L0SX\JX) \JMSXSrf 2 ). (G.110) 

V 2iJ 1 


How does this relate to normal spins? Well, for those, 

WSMsSiSz) = ( s iMi S 2 M 2\ SM s)\9(f>S 1 M 1 S2M2). (G.lll) 

M\ M 2 


But the relation of these to helicity is complicated, since in general every M\ is 
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connected to every Ai. In particular, 

{9(j)SiXi S 2 A 2 | d<j)S 1 M 1 S 2 M 2 ) 

= (OOSi\iS 2 \ 2 \R^ 0 R^Q 0 \OOSiMiS 2 M2) (G.112) 

= (00SiXiS 2 X 2 \e iejy e i<l>Jz e~ i(l>Lz e~ ieLy \00SiA / IiS 2 A / I 2 ) (G.113) 

= (00Si\iS 2 X 2 \e i9Jy e~ ieLy e i ^ Sz \00SiM 1 S 2 A / I 2 ) (G.114) 

= (00S' 1 A 1 5 2 A 2 |e ws “e^ s *|005 , iM 1 S , 2 M 2 ) (G.115) 

= (00100)Ai| e~ l ^ Sz e~ ieSy \S\M\)(S 2 - X 2 \e ieSy e i4>s *\S 2 M 2 ) (G.116) 

= (S 1 \ 1 \e ieSy \S 1 M 1 )(S 2 - X 2 \e idSy \S 2 M 2 )e i<t>{Ml+M2) (G.117) 

e^ Ml+M2 )4 iiAi (d)4 2j _ A2 (d). (G.118) 

This is easily generalizable, to get 

WtSiMWSiM,}) = n e> * M ' d M,J 0 )- (G.119) 

i 

Well, at least this makes sense for 1 and 2 particle states. 

This has an implication for parity, II. We get 

{9<f>{S i X i }\9 ( f>{S i M i }) = {9(l){S i X i }\UU\9(j){S i M i }) (G.120) 

= 5(n - 9 tt + 4>{Si - Ai}17r - On + 4>{SiMi}) (G.121) 

= (G.122) 

i 

= i J] e W(- 1 )«.(- 1 )a-M.4 iA (9) (G.123) 

i 

= (G.124) 

i J 2 

= (G.125) 

i 

So we see that 5 = ITi(—1) S S so 

U\9(f)SX) = (-l) s \n-dn + (j)S-X) (G.126) 

n|d05iAi^ 2 A 2 ) = (-1) Si+S2 |tt - 9n + - Ai^ 2 - A 2 ). (G.127) 
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Continuing on, we can derive other kets from the normal spin kets, 

1 f 2n 

\9M l SM s SiS 2 ) = -= / d^e^e^SMsS^), 

V27T Jo 

| LM l SM s S 1 S 2 )= I dcos9Y LML (9,0)\9M L SM s S 1 S 2 ), 
\JMLSS 1 S 2 )= (LM l SMs\JM)\LM l SMsS 1 S 2 ). 

Ml, Ms 


(G.128) 

(G.129) 

(G.130) 
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VITA 

“The life of every individual, if we survey it as a whole and in general, 
and only lay stress upon its most significant features, is really always a 
tragedy, but gone through in detail, it has the character of a comedy. For 
the deeds and vexations of the day, the restless irritation of the moment, 
the desires and fears of the week, the mishaps of every hour, are all through 
chance, which is ever bent upon some jest, scenes of a comedy. 

“But the never-satisfied wishes, the frustrated efforts, the hopes un¬ 
mercifully crushed by fate, the unfortunate errors of the whole life, with 
increasing suffering and death at the end, are always a tragedy. Thus, 
as if fate would add derision to the misery of our existence, our life must 
contain all the woes of tragedy, and yet we cannot even assert the dignity 
of tragic characters, but in the broad detail of life must inevitably be the 
foolish characters of a comedy.” 

—Arthur Schopenhauer 
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